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ABSTRACT 


This paper deals with the evaluation of the reflection coefficient for 
backward scattering of an electromagnetic wave from a column of ioni- 
zation having cylindrical symmetry. The application of the theory to 
meteors thus depends on the assumption that a meteor instantaneously 
forms a narrow column of ionization which subsequently diffuses radially. 

The theory predicts two qualitatively different categories of meteor 
echoes, depending on whether the electron line density, «, is greater or less 
than about 10!2 per cm, the properties of which may be summarized as 


follows : 


atl no 1014 

If the incident electric vector is parallel to the axis of the trail (parallel 
scattering) it is sufficient to assume that the electrons scatter coherently 
and independently, each contributing an amount 47(e?/m c*)? to the back- 
ward scattering cross section. As the column expands radially, the echo 
amplitude decreases, the form of decay being exponential if the electron 
density is a Gaussian function of radius. 

When the incident electric vector is normal to the column (transverse 
scattering) the echo amplitude differs from the above in that the reflection 
coefficient may pass through a resonance, followed by the exponential 
decay. The magnitude of the resonance (measured as the ratio of the 
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resonant reflection coefficient to that for parallel scattering) depen ds mainly 
on the diffuseness of the column (i.e. on the gradient of electron density at 
the radius of critical density). It has been evaluated for an electron 
density varying with radius as n=ny exp [—(/79)°] and is found to range 
from 2-0 when s=2 (Gaussian distribution) to 10 when s=20. The echo. 
amplitude in both cases is proportional to «, while the echo duration is. 
independent of «. 

The short duration echoes which form the bulk of the experimental 
observations appear to have the above forms. 


(2) «>10!" 

In this case the column is expected to reflect similarly to a metallic 
cylinder. Except initially, when the parallel reflection coefficient exceeds. 
the transverse, the reflection coefficients are equal and substantially 
constant. The duration is approximately proportional to «, and the peak 
echo amplitude is a slowly increasing function of « (~a1’*). 

These predictions are in agreement with the experimental data on long 
duration echoes. 

Some particular electron distributions, including a small diffuse sphere, 


which exhibit multipole resonant scattering are discussed in Appendix IT. 
and laboratory experiments which may be of some interest are suggested. 


$1. [INTRODUCTION 


THE passage of a meteor through the atmosphere is accompanied by the 
emission of visible radiation, and by the creation of an ionized column of 
air, both of which processes reach a maximum at heights between 80: and 
120km. The latter effect has been made the basis of an intensive study 
of the properties of meteors and of meteor streams, notably by Lovell and 
co-workers in Great Britain, and by Millman and McKinley in Canada 
(see Lovell (1948) for a review article). 

It is believed that the major part of the ionization arises from collisions 
between evaporated neutral molecules of the meteor material and air 
molecules leading to a narrow column of electrons and positive ions which 
rapidly diffuses outwards. 

Lovell and Clegg (1948) derived the radar equation for such an ionized 
column, assuming it to be narrow compared with the wavelength and to 
have a small enough electron density to enable the incident wave to pene- 
trate the column without significant modification. In spite of these 
limitations, the theory gives an adequate explanation of many of the 
observed phenomena. Herlofson (1951) has shown that when the 
column is somewhat more than critically dense, this formula is applicable 
only when the electric vector of the incident wave is parallel to the axis 
of the trail (parallel scattering) and that in general the reflection coefficient 
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may exhibit a critical dependence on the orientation of the electric vector 
with respect to the column. Theoretical studies by Herlofson (1951) and 
by Feinstein (1951) predict that if the incident electric vector is normal to 
the axis (transverse scattering) the reflection coefficient will exhibit a 
resonance when the electron density is somewhat above the critical value 
provided that the transverse dimensions of the column are small Chmpered 
with the wavelength. 

The magnitude of the resonance (measured as the ratio of the resonant 
reflection coefficient to that for parallel scattering) will depend on the 
radius of the trail and more critically on its diffuseness, i.e. on the gradient 
of electron density near the region where the column is critically dense. 
Herlofson has considered a model in which the electron density is constant 
from the axis to a given radius and then decreases linearly to zero over a 
small radius interval, while Feinstein has derived the reflection coefficient 
for a column in which the electron density decreases linearly to zero from 
the axis. 

The purpose of this paper is to obtain the scattering parameters of a 
meteor trail without restriction on the electron line density or on its radius. 
Part is devoted to the derivation of an approximate method of solution of 
the wave equation for transverse scattering when the electron line density 
is less than 1012 cm~ and the radius of the column is small compared with 
the wavelength (which are the conditions for the existence of a resonance). 
The method is applied to an electron distribution given by 


HN OXDM (Ti ale le werd week Bell yates (2) 


where n=electron density and r—radius. 

The important physical case of s=2 (Gaussian electron distribution) is 
studied in detail while other values of s are considered in the neighbourhood 
of the resonance. 

The method is shown to be applicable to ionized clouds of geometrical 
configurations other than the circular cylinder. An example is the spheri- 
cal cloud considered in Appendix II, where resonance effects are again to be 
expected. A later paper will deal with the problem of scattering from a 
diffuse elliptic cylinder which can be treated on similar lines and which may 
have a physical significance owing to the possibility of anisotropic diffusion 
in the earth’s magnetic field. 

Both the Lovell—Clegg result for parallel scattering and the above 
mentioned theory break down when the electron line density, «, reaches the 
order of 1012cm-1. It is shown that for larger line densities than this the 
expanding meteor column reflects in a similar manner to a metallic cylinder 
of appropriate radius. On this basis, the reflection coefficients have been 
evaluated for these high densities and for both parallel and transverse 
scattering. By joining smoothly the two solutions in the vicinity of 
a—10!2 a complete solution for columns having a Gaussian electron 
distribution has been obtained, which is not limited as regards the electron 
line density or the radius of the column. 

B2 
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§2. GENERAL THEORY FOR A CIRCULAR CYLINDER 


Transverse Scattering 


The existence of a resonant process for the scattering from a column of 
radius small compared with the wavelength may be demonstrated by 
considering the cylinder as one of dielectric constant, 

/ 2. 

() =e ie Wo TOS aR 
immersed in a uniform transverse electrostatic field. The units in (2) are 
Gaussian (as throughout the paper), n(r)=electron density as a function 
of radius, k=27/A, A=free space wavelength, e, m and c are the charge 
and mass of the electron and the velocity of light respectively. This is a 
purely formal treatment and is not in itself, of course, physically realizable. 

For the simple case of a homogeneous cylinder, along the z-axis of 
cylindrical polar coordinates, having «(r)=«, r<a; «(r)=1,r>a, and ina 
uniform field E, parallel to the y coordinate axis, standard text books give 
the following expressions for the electrostatic potentials, V,; and J,, 
interior and exterior to the cylinder respectively, 


yas 

V=—By (r+ F< S) sin oo Tito fase See 
2r sin d 

Vji=—£E, ae + Atala a, el 5) oe ee 


The first term in (3) corresponds to the applied field and the second to 
the distortion introduced by the cylinder. This latter tends to infinity 
when «= — 1 which, in the electrodynamic case, we might thus expect to be 
the condition for resonant scattering. 

Herlofson (1951) shows that the existence of resonances may be demon- 
strated using a similar model to the above. He evaluates the restoring 
forces between the surface charges of a polarized ionized cylinder and finds 
that the system is a resonant oscillator when «=—1. 

We may likewise treat parallel scattering. The internal field will now be 
equal to the applied electrostatic field, Hy, compared with 2#,/(1+-«) in the 
above case. Denoting the dynamic reflection coefficients by g, and J, we 
would expect that on this simple model 


9 [Gy eit). =. ee, a Oe (5) 


The infinity at «=—1 will be removed physically by the small but finite 
column radius, collision damping etc. 


The generalization of the above treatment, including its extension to 
a non-homogeneous column involves the solution of the equation 


div (« grad V)=0, ee eo es re) 
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which, in cylindrical coordinates and on separation of variables, leads to 


a T2) = mV in 
rey , sary . . ‘ ° 7 6 (7) 


where 
V(r, $)=ZV (1) sin meg ; Faso) tigen iat Be bs} 


kK isa function only of radius, and m is a positive integer or zero whose 
values specify the various modes of excitation of the column. 
An exact study of transverse scattering gives the equations 


- ae ice , 2 
1s (; — <) A,,’/+ (ie 7) Ey) so aeawe i (OF 
A(r, $) = 2: Healt) COS TNO .c6 eae eee ete LO} 


where the magnetic vector has only a z-component, H (primes indicate 
differentiation with respect to r). The term k2«H,, takes account of 
phase variation in the field and (7) and (9) can be shown to be equivalent 
if it is negligible, i.e. if we place upon r and « the restriction | (kr)? |<m? 
or, for the dipole mode, 
|x(kr)?| <1. RN ee era ee 8 
With this restriction, using E=—grad V, curl H=jk«E, where E is. 
the electric vector, we obtain 


her Vm, 
Vin=1H ny |(gkem). 


The solution of (7) in the region outside a column (i.e. where «=1) whose 
parameters satisfy (11) is 


(12) 


TONG Anred Dy ian Meee ery en ass -= ©. (Lo) 
i.e. gies) ci Ani Dtawe) sors jared Uy! Vath kop 
where A,, and B,, are constants and m=1, 2, 3,.... Except in special 


cases, A,, and B,, must be determined by numerical integration of (7). 

The exact solution of (9) for the region outside the column involves the 
Bessel functions of the first and second kinds, J,,(kr) and N.,,(kr), which 
reduce to the independent solutions r” and 7~” respectively when (kr)? <1. 

Let us now consider a plane electromagnetic wave of unit amplitude, 
incident in the negative x direction and polarized in the direction of the y- 
axis. Omitting the time varying factor, exp (jwt), the expansion of the 
incident wave as a series of cylindrical waves is 


He=exp (jhe) =Jolkr)-+227"F ,(ker) cosmd, . « (15) 
1 


remembering that the magnetic vector has only a z-component. 
The reflected wave may be written 


Hyon =StplJnlbr)—JNnlbr)] cos mp, « « + + (16) 
(0) 
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which, at large distances from the column, (kr)? > 1, describes an outward 
travelling cylindrical wave whose value on the positive w-axis, (6=0). is 


(Arent) =0)= Baal [2/(akr)}'? exp [—j(kr—m/4)]. —- (17) 


The total external magnetic component is 


Hop 2 Hm COR MDs) oe ROG are nee 
0 
where 
H,=(29™ btn) (kr) —FtmN m(Ar), m=1, 2, 3, a leg | (19) 
Hy=(1+to)J (kr) —JtoN o(kr). j 


When (kr)?<1, J,,(kr) = (kr)™/(2™m !), N;,(kr) = —2™(m—1) !/[a(kr)™], 


m=1,2,3,.... Using these approximations and comparing (19) with 
(14) we obtain, for the transverse scattering reflection coefficients, f,,,, 
24m ! = 1 AeA 
29 ae ge Oe ae Raa 20 
7c =~14) an Bet ilk ete cascade) 


We have not considered m=0, however it is not difficult to show that to 
the degree of approximation of (20), tj=0. 

With (20), the reflection coefficients, t,,, may be evaluated for any 
electron distribution, n(r), by integrating (7) using the appropriate 
function «(r). Note that (20) is subject to the two conditions | «(kr)?|<1 
within the column, and (kr,)?<1 where 7, is the radius of the column. 

The total reflection coefficient for transverse backward scattering, 
denoted by g,, is, from (17), 


2 . 
2) al 2 
0 


g,= i ope cba isk eiiealt i 
When 4,, is real (which occurs only for a non-diffuse trail), the mth mode 
exhibits a resonance (i.e. | t,, | is a maximum) when 4,,=0. In general the 
condition for resonance is Re (A,,)=0. 

The scattering diameter of the column is defined as the width of the 
equivalent strip of material which reradiates, as an azimuthally isotropic 
cylindrical wave, all of the energy incident upon its surface, and which 
produces a backward scattered energy flux equal to that produced by the 
actual column. Denoting it by d, and using (17), it is given by 


* 
2 


kd=4 


co 
am 
*) tin 


mxdg 30 Gore guie \oattee nn (0) 
The total cross-section for backward scattering, denoted by a, is 
o= 89 ,*r/k=2dr. See ital = Sa Hier mest #45) 


if we have « electrons per cm, then the reflected amplitude at a large 
distance, r, from the column is equivalent to that produced by «4/(rd) 


electrons radiating coherently. The contribution of each electron to the 
scattering cross-section is thus 


Ua 4g ,7/(702), e4 fe 


Seen eh sheen (24) 
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We now consider a plane wave of unit amplitude incident in the direction 
of the negative a-axis and polarized parallel to the z-axis. 
The appropriate wave equation is now 


VE+kxeH=0, 
and hence 
d . dE, 2 2 
De ral le +[k(kr)?—m?]#,,=0, 
where E=ZE,, cos md. 
0 ' 


The expansion of the incident wave into a set of cylindrical waves is 
identical with (15) replacing H by EF. The reflected wave, the total 
external field, and the scattering parameters, will likewise be represented 
by (16), (17), (18), (19), (21), (22), (23), and (24), substituting J,, for t,, 
and g, for g, to avoid confusion. 

When |x(kr)?|<1 throughout the column, there will be negligible 
attenuation of the wave as we pass from the exterior to the axis, the only 
significant mode will be m=0, and the Lovell—Clegg scattering formula 
will be valid. The reflection coefficient, in this case, becomes 


ak (ee) 
n=l l= |. r(1—r)dr a 
=7ae"/(mc?), 
and kd=[2na.e2/(me?)|?, 


o,=47(e?/mc?)?. ; 
The latter is just the cross-section of a free electron for backward scatter- 
ing, completing the identity with the Lovell—Clegg formula. 


Relation between the Reflection Coefficient and the Signal received from a 
Meteor Column 

In the above discussion we have, for simplicity, considered a plane wave 
incident normally to the axis of the cylindrical column, leading to the 
reflected amplitude given by (17). 

In meteor experiments we normally have a transmitter and receiver, 
matched through a common T-R network to an aerial, distant & from the 
meteor column, hence the incident wave is spherical. The curvature of 
the incident wave front results in a central Fresnel zone, seen from radius 
R, of effective length \/(AR/2) instead of 1/(AZ). In this case the ratio of 
the power flux in the reflected wave at the aerial to the incident flux is 


De D0 (tick) sees oso ne pce (26) 

and the power delivered to the receiver terminals is 
PGP : i 
R= 35-4739 watts, ree Fut ae Age 9 8 


where P>=transmitted power and G=aerial power gain referred to an 


isotropic source. 
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§3. HomMoGENEOUS CYLINDER 


, The case of a cylinder of constant electron density represents a consider- 
able idealization as far as a physical meteor trail is concerned, however, it 
illustrates certain important features of the scattering process. The more 
general electron distribution of eqn. (1) is studied in §§ 5 and 6. 

We will consider a dielectric cylinder of radius a and of constant 
dielectric constant x. 
The solution of (7) external to the cylinder is now 


Va=Ae | + (a)anr |. Re ee | 


which reduces to (3) when we put m=1 and introduce the angular factor. 
Thus 4,,/B,,=[(1+«)/(l—«)]a-”, and 


24m !(m—1)!4™" /1 
Se sy pg ea (29) 
tn m(ka)2™ 1—x 
We see that |t,,| is a maximum’ when «=—1 and thus, at resonance, 


| tm lres==2 for all modes. 

Except near resonance only the dipole mode is significant, and even then, 
as we shall see, all physical damping processes tend to reduce the resonant 
reflection coefficient for higher modes below that for m=1. In this case 
g, = |t, | where 


5d 4) l+k 
-#=14 To (). eee pee EI 


and 9g) ,..=2 independent of the radius of the column. ‘The resonant 
scattering parameters are thus 


Aros 8A/77, Ores= L6rA/z. oy sk ok ee 
For parallel scattering, from (25), 
Gy (ka?) (L—K)/4.00 (oe 3 st ee ee 


Ignoring higher modes than m=1, the ratio of reflected amplitudes for 
the two polarizations is 


Pes .m(ka)’(L—K) _ —1 
qe penne (l+x)| , Siig Re RR IN ES 


which, except in the neighbourhood of the resonance, remembering’ the 
condition (ka)? <1, reduces to the simple result (5). 
At resonance, 


OR Poet ua clad i ene ee en hey: 
Fig. | shows g/g, as a function of (1—x) for various values of ka, and in fig. 2 
the resonant ratio is plotted against ka. The broken part of this latter curve 
indicates the behaviour over the region in which the assumptions break 
down. This ratio, unlike the resonant reflection coefficient itself, is a 


function of ka. Tt is a measure of the magnitude of the resonance, which 
decreases with increasing column dimensions, becoming unity when ka = 1. 
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ee condition for resonance to be possible, namely that «—=—1 and 
(ka) ‘<1. becomes 20g mec?*)<1 or, substituting numerical values 
n<2 1022. : 


Hig 


of 1-0 10 
z aan (Oe) 


Resonance curves for a homogeneous cylinder of radius a and dielectric 
constant x. 


Fig. 2 


oO-l a2 O-5 (-O 2-0 
—_> ka 


Resonant polarization ratio vs ka for a homogeneous cylinder of radius a. 


(ka)? |>1, the above formulae break down and the 


When «<0 and |« ate 
n this case, 


column tends to behave as a metallic cylinder of radius a. 
when (ka)?<1, the condition is «<—I1 and the reflection coefficients 


become . : 
Gam | ty [= ba)* 2a eee oer ode eee ee) 


and | 
gy = fo [=| 1—§(2/7) In (yha/2) [> (36) 
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where y=1-781. When (ka)?> 1, the scattering diameter is obtained from 
geometrical optics to be d=7a, and thus, from (22), 

9, =G)— (reals Oe 
If we let (ka)? tend to infinity, the cylinder tends to a plane and the 
condition for metallic reflection is «<0 as we should expect. Using (35), 
(36) and (37) the curve of fig. 3 has been plotted by smoothly joining the 
three branches in the vicinity of ka=1. 

Factors, other than the radiation damping of the cylinder oscillations 
connected with the finite trail dimensions, which tend to reduce the 
magnitude of the resonance are radiation damping of individual electron 
oscillations, collision damping, and the diffuseness of the column, so that 
(34) gives only an upper limit. We will consider only the first two in 
relation to the above model, leaving the third, and most important (at 
least for meteor echoes) to §5. The effect of both is to introduce a complex 
term into the expression for the dielectric constant. 


Fig. 3 


Reflection coefficients for transverse (g,) and parallel (g,) scattering from 
a metallic cylinder of radius a. 


If radiation damping of individual electrons was the sole mechanism 
limiting the resonance, then each electron would contribute an amount 
9A?/(477) to the scattering cross-section. Comparison with (34) shows that 
radiation damping could only be significant if AX8/(3«), which is, in fact, 
an impossible condition. Thus this effect will never compete with that of 
finite radius. 

Collisions between electrons and neutral molecules leads to the following 
expression for the complex dielectric constant (when v<w) : 


aT 47re2n 1 ny 
mae) 


- (38) 
=K 1 —Jke 
where y=collision frequency, w=angular radio frequency. Resonance 
occurs when x,;=—1 and the resonant reflection coefficient for the dipole 
mode is . 
2) aay 


~ 7, lta hae eee ene ae (aY) 
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Seam are thus significant for this mode if the radius at resonance 
satisfies 


GRE y S| (hac) Wen een. 5 agit 10\ 


For A=4 metres and using plausible figures for y in the meteor ionization 
region, we find the radius satisfying (40) to be of the order of only a few 
times the molecular mean free path. Since our assumption that the 
column may be regarded as a continuous dielectric will restrict the radius 
to something greater than this, collisions will in general not be important. 

We will not pursue further the discussion on collisions (except to say 
that the effect will be much larger for the higher modes), since the 
discussion in §5 reveals that in an actual meteor column the effect 
of its necessarily diffuse boundary will be predominant. Collision 
damping considerations may, of course, be important in laboratory 
experiments with discharge tubes of the type made by Denno, Prime, and 
Craggs (1950) and by Romell (1951). 


§4. Fine STRUCTURE OF THE PLASMA RESONANCE SPECTRUM 


Although modes other than m=1 are unlikely to be of importance in 
connection with meteor echoes, there is a matter of some theoretical 
interest concerning the splitting of the resonances. As we have seen, for a 
homogeneous cylinder, the resonant modes are degenerate, all occurring 
when «= — 1 (a higher degree of approximation would in fact have shown 
a slight splitting, tending to disappear when ka — 0). - The reason for the 
degeneracy can be seen by writing r”=x when (7) becomes 


a eae Vin 


with the solution V,,=A,,7+B,,x~1 outside the ionized region. In general 
« will be a different function of x for each mode, hence the condition that 
the real part of A,,, be zero will be different and the resonant modes will not 
be degenerate. Only for a homogeneous cylinder is « the same function of 
x for all m, namely, 
«=const., Cty =a”, 
K(x) = ‘5 LIL, 


and we have 4,,/B,=% 2(1-+«)/(1—x), irrespective of m. 
Examples which show non-degenerate resonant modes are given in 
Appendix II. 


§5. Dirruse CoLuMN oF IONIZATION 


We expect the passage of a meteor through the atmosphere to form an 
jonized column with initial radius of the order of the molecular mean free 
path. The column will then expand according to the laws of diffusion, 
however, since its initial temperature may greatly exceed that of the 
atmosphere, the diffusion coefficie:tt, D, may be a function (decreasing) 
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of radius and time. For simplicity we will assume that D is constant, in 
which case the solution of the diffusion equation leads to the electron 
distribution 

a 


ne i seh Ae ee 


= =e oP [—(r/ro)?] 


where 7,2=-4Dt and we have made the simplifying assumption that at time 

{—0 there is a line distribution of ions and electrons. The effect of a radial 

temperature gradient may be to cause the column to have a somewhat 

less diffuse boundary than Gaussian, but without more information of the 

process involved it is difficult to estimate the extent to which this will occur. 
For the electron distribution (42), we obtain 


«(r)=1—f exp [—(7/")"] 
fHda(e'{me?)|(kr)? | 
which is sketched in fig. 4. 


fone Riga het Fapase 


Fig. 4 
1-0 i aw 
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B 
(e) 
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Dielectric constant of a cylindrical column with electron density a Gaussian 
function of radius. 


The differential eqn. (7) now has singularities at r=0 and r=a, where a 
is the radius at which x0. The first is dealt with by requiring V,, to be 
well behaved at the origin, but the second requires a more careful study. 

Consider a thin region around r=a. Putting z=r™ in (7) we can obtain 


ad AV in i d es 
da \“ “dx =—«3,(=") ; 


which, in this region is approximately 
dV 


df dVy Kk dV, 
mee | 
dx dx a den 
i.e. —=— = S cp [—a# = 
e Fea is (const.) exp [—a/a”] = const. when |7—a|<a, 


ava 
K ENS = h, igen et ey Tate nes (44), 


or 


where / is a constant. 
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Vinp— Vesa | ‘ kadar. Me ake str B, os eco ae (45) 


cD we he K by (38), i.e. we consider collisions, then in the region near 
cai we have x, PAP exis 1 (@) and x, remains substantially constant. 
utting K=(r—a)x«,'(a)—jry into (45) we have 


Ving a Ve ake Ln (=) 


/(a) BA 
iwyeh 1 ant a 
= | nfB)+im—0) |, +» (46) 
where Ln denotes the complex logarithmic function and 


0<6=tan4 | a] oe 


KikKip tke” 


155 
Ka 


When x, > 0, 


ae h icant 
V mB y mA— x’ (a) E | a | +in | ? Oe ey Oy (47) 
Fig. 5 
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Herlofson’s model of a diffuse electron distribution. 


i.e. even when the collision frequency vanishes, we are left with a finite 
imaginary part of V,,. The physical explanation of this is that the term 
In | xp/«, | has a maximum at r=a which becomes infinite when «,=0, 
i.e. the electric field intensity becomes infinite and a negligible collision 
frequency can have a finite effect. The electric displacement, «E, remains 
finite and continuous over the region, in fact (44) is just the condition that 
it be so. The integral of the energy density over the region is also finite 
and represents a negligible proportion of the total electromagnetic energy 
in the column. 

By a procedure of numerical integration outlined in Appendix I, and 
using (47) in the vicinity of the singularity, the ratio A,,/B,, may be 
evaluated for a given electron distribution. 


Herlofson’s Model of a Diffuse Column 


Herlofson considers a column which is homogeneous except in a thin 


region near 7=a (fig. 5). 


14 T. R. Kaiser and R. L. Closs on the 


Within the homogeneous region we may put V,,=7" thus, from (44) 
h=mxa™—! and from (47) 


ma”. 
Ue oe tating re k : 2 eas 
aV»R =m«Ka™ 
Now, outside the column, V,,=A,7"+B,,7~™, 


i.e. Vi, pe Age + Bat 
/ m —m 
aVing =m(A,,a = Die ) 


1 


Ip 2 te Rada By 


Equating the eqns. (48) with (49) we obtain 
Any om (L+H) + Line (ax’)](jr—ln | «) on 
(1—«)+[m«](ax’)](gr—In | |) , a 
1) when ax’(a)>1 and kx~—1. 
A (1+«)+jmk/(ax’) 


Mm ~ q—2m 


B 1—k 


m 


Bos 
which reduces to (5 
(51) 


The value of [ax’(a)]~! at resonance measures the degree of diffuseness of 
the column. 
Substituting (51) into (20) gives 
2jm arg .m ! (m—1)! 4" (1+«)+-jm7K/(ax’) 
ea tenis J a(ka)2™ (1—k) : 
which is the result obtained by Herlofson (1951). 
The effect of the diffuseness on the dipole resonance is significant when 


(ka)?ax' (a) <2, 
or ax’ (a)(e?/me?)x< 1, eee A ee 
: , [ 3-5, a~lol2 
Le. ak’ (a) x 


| (350, er ore LO 


(52) 


The electron distribution of fig. 6 corresponds to 


n(r)==1—f exp [a= (7ito) unk eet eC 
which is similar to the Herlofson model when s>1 and is Gaussian when 
s=2. We have (a/r))*=Inf, ax’(a)=s Inf. 
At resonance, when s>1, f=2 and 
aie (a)2=0-698.9 0 eee een on 
Thus, unless the distribution has a significantly sharper boundary than 
Gaussian, the diffuseness will in general be the major limiting factor. The 
diffuseness limited reflection coefficients using (20), (50), and (55), are 
2(ka)” i 
[tn I 05 Gon Haar | Zein nT | eee ee ere O0) 
(ka)? 


Idan epee 5) | 2s In 2—Jjr | 


~a(e?/me?)|2sIn2—jr|. . . . . (57) 
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- Itis clear from (56) that the diffuseness tends to eliminate the higher modes. 
Using (25) and (57), ignoring modes other than m=1 the ratio of the 
resonant reflection coefficients for the two polarizations is 


(91/91) ves | 0-448—J |, iF ettey 1 Preeti 7(58) 
which is plotted in fig. 9, curve (b). 

The resonance phenomenon will result in the column acting as a 
frequency selective reflector and hence a bandwidth can be assigned to it. 
We can show that the reflection coefficient is reduced to 1 /4/2 times its 
resonant value when the angular frequency differs from the resonant 
frequency by dw given by - 


w $/OW—= 2(9,/9,) res* . . . . . . . (59) 


Fig. 6 


Electron density n=, exp [—(r/79)*] for various s. 


§6. GaussiAN ELecTRoN DISTRIBUTION 


Because of the physical importance of this distribution we have con- 
sidered the reflection problem in detail, and over a wide range of « and 7y. 


(a) Dilute Column , , 

The column is termed dilute when (1—x)<1 throughout, in which case 
the incident wave penetrates the column without modification. Sub- 
stituting for « from (43), the condition for diluteness is 


(krp)?> 11-24 x 10-¥a. el Pe oe, ee OU) 

The reflection coefficient is independent of polarization, being given by (25) 
when (kr,)2?<1, and in general by 

Ju= 9, =72( e?/me’) exp [—(k7y)?], Sir Ore OLS 

where the exponential decay results from phase differences within the 


column (Herlofson 1951). 
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(b) Line Density «<2-4< 10% cm™ 
Eqns. (7) and (25) are valid only when (kr))?<1 and if | «(kr)?|<1 
throughout the column, which latter condition requires «<2-4x 10”. 
Subject to these restrictions, from (25) 


gyn (lore)f/40) 0 ae eee 


Transverse scattering requires the solution of (7) for the dipole mode, 
m=1, which we have evaluated numerically for a range of values of f. 
Because of the diffuse character of the column both A, and B, are complex 
quantities and within the limits imposed by (hry)? <1 the diffuseness is the 
main limiting factor, while for (kr))?>1 the resonance disappears. 

Thus 


—] 


bas bi: 4) A’ 
a2 |t 122 [14+ ea 


(63) 


+ a = tp» —1 —— Ply. 
where A,=A’r, 1, Bj =B'P9. 


Polarization ratio vs f for a Gaussian cylinder, when the resonance is 
diffuseness limited. 


From (62) and (63) the polarization ratio when the resonance is diffuse- 
ness limited is - 
ore ne : 
Il =F | AI Sh Wr) 1s a net wie ERs aot) Me (64) 


The values of A’ and B’ are given in Appendix I together with a 
description of the method by which they were obtained. |The resonance 
condition is Re (A’)=0, which occurs when f= 2-4, i.e. when the dielectric 
constant on the axis is —1-4. Using the computed values of A’ and B’ 
with (64), we have plotted g,/g, as a function of f in fig. 7 which shows that 
the resonant ratio is reduced to 2-0 by the diffuse boundary. | | 

From (43), (62), and (63) we can determine g,/g, as a function of (kr)?/« 
which is given in fig. 8 for various values of x. Since 1,2—=4Dt f 8 
represents the form of the time envelope of the echo obthined with ree 
verse polarization from a column of « electrons per cm. | 
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£ te have also evaluated the solution of (7) in the vicinity of the resonance 
rata bare general distribution, (54). The diffuseness limited resonant 


Fig. 8 
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Polarization ratio vs (kr9)?/a=4k?Dt/a for various values of the electron 
line density, a. 
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Cylindrical column with electron distribution =”) exp [—(r/79)$] 
(a) Diffuseness limited resonant polarization ratio, computed as described 


in Appendix I. . . 
(b) Diffuseness limited resonant polarization ratio derived from 
Herlofson’s model. | 
(c) Resonance condition, fres V8 8, (L—f= axial dielectric constant), 
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ratio is plotted in fig. 9 (curve (a)) as a function of s, together with the 
values obtained from the Herlofson model (curve (b)). The value of f at 
resonance, as a function of s, is also plotted in fig. 9 (curve ‘(c)), which | 
indicates that f,,, > 2 when s > ©. 

When «<2:4X10!2 and (kr))2>1, we see that (60) is automatically 
satisfied and the reflection coefficients follow (61), completing the descrip- 
tion of echo amplitude as a function of (kry)?. 


(c) Line Density «>2-4 10" 

We may now have the conditions «<0 and |«(kr)?| >1 over a region 
within the column. For a homogeneous cylinder, as we have seen, this is 
just the condition for metallic reflection. We will assume here that the — 
Gaussian column behaves as a metallic cylinder of radius r, given by the 
larger of the two values satisfying 


x(r,) (or, P= 1 ne oe ee 


which expresses 7, as a function of 7, for given «. 

The justification for this assumption may be seen as follows. When 
(kr,)?>1 it is obvious, since 7, becomes the radius at which «=0 and the - 
problem can be treated by geometrical optics with reflection taking place 
at the surface of critical density. When (kr,)?<1 we still expect the 
column to behave as a metallic cylinder and (65) must at least give the 
order of magnitude of the equivalent radius. In fact the reflection 
coefficients derived using this radius join smoothly in the region of 
“=2x 10!" to those obtained in (b) above. 

Solving (65) for r, and using fig. 3 with r.=a we find the reflection 
coefficients as a function of (kr,y)?. 

When the column has expanded to a radius satisfying (61), exponential 
decay will set in, and we can join the solution as obtained above with that 
of (a). This is not important when «>10!8, since the coefficient as 
determined above remains large compared with (60) out to some radius and 
then drops sharply. 

Approximate expressions for the maximum reflection coefficient and 
duration, 7, are 


Puiag — UL LOSSES ee oe tec aa 
412 Dr 112450108 eal eee Lemay 


(4) No Restriction on Line Density and Radius 


By joining smoothly the solutions for the two ranges of line density to. 
each other and to the dilute column solution we obtain the complete 
result given in fig. 10, which gives both transverse and parallel reflection 
coefficients for a Gaussian electron distribution, without restriction on line 
density or column radius. In fig. 11 we have sketched some of the curves 
of fig. 10 with linear scales. These curves give the predicted form of the 
time envelope of the echoes. ) 
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———— Parallel scattering. 
—----- Transverse scattering. 
Reflection coefficients vs (kr»)?=4k?Dt for a Gaussian cylinder having electron. 
density a, plotted with logarithmic scales. 
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-—__—— Parallel scattering. 
== -- = Transverse scattering. 


As for fig. 10, but with linear scales. These curves give the expected 
amplitude-time envelope of the received echo, assuming the meteor 
instantaneously to form a narrow column of electrons which expands 
radially at a rate determined by the diffusion coefficient, D. 


C2 
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§7. DiscussIon 


The complete scheme outlined in fig. 10 has obvious application to the 
analysis of experimental results, particularly to the derivation of diffusion 
coefficients, electron line densities, etc. from observed echo amplitudes 
and durations. 

Certain important features are immediately apparent. 

As predicted by Herlofson we may expect a resonance phenomenon for 
transverse scattering to be exhibited by columns having line densities of 
less than 10! cm-!. The resonant amplitude will depend on the diffuse- 
ness of the column, being about twice the parallel scattering amplitude for 
a Gaussian electron distribution and increasing for trails with sharper 
boundaries. Such resonances have been detected by Clegg and Closs 
(1951) who find resonant ratios of the expected order of magnitude. 
Larger ratios than 2 may also be associated with an anisotropy in the 
diffusion coefficient caused by the earth’s magnetic field. A solution of the 
scattering problem for a thin diffuse elliptic cylinder will be given in a 
second paper, when it will be shown that quite a small anisotropy may 
lead, under certain conditions, to a significant increase in the resonant 
ratio. 

We mentioned, in connection with collision damping, that we expect the 
column to be formed with a certain minimum radius and it was sufficient 
for the previous argument to consider this as set by the finite molecular 
mean free path. Another limit is set by the finite meteor velocity, since 

_if we consider its formation we find that the column deviates from the ideal 
picture of an infinite cylinder. We cannot, in fact, consider the column as 
an infinite cylinder until it at least fills the first Fresnel zone of length 
V/ (AR/2), situated symmetrically about the normal from the observer to its 
axis. In the time taken by the meteor to traverse this distance, diffusion 
will have taken place and the trail will be in the form of a paraboloid of 
revolution, of average square radius (over the Fresnel zone) 


To min=D/(2AR)/», Buuit 0 wows oS (68) 


- 


where v is the meteor velocity. We may regard 7,,;, a8 an effective lower 
limit to the column radius set by the finite velocity and for a typical ease, 
with A=4m, R=150 km, v=40 km/s, D=3 x 104 ¢.g.s. units (as determined 
from echo durations at a height of 85 km) we get 77 ~28 em, which consider- 
ably exceeds the mean free path at this height. For the above parameters 
(kro)*min~9-2, thus, referring to fig. 10, the smallest meteor exhibiting a 
polarization effect will be one producing a few times 101 electrons per cm. 
‘The apparatus of Clegg and Closs (1951) was sufficiently sensitive to detect 
a reflection coefficient of about 0-1, i.e. a line density of 10" with parallel 
scattering. It is not therefore surprising that a substantial number of the 
short duration echoes appeared to exhibit ratios, 9,/9,» greater than unity. 
Fig. 10 indicates that we would expect initial ratios of less than unity for 
some echoes and a few such examples have been observed. 
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4 It of interest to compare the durations of echoes having different line 
ensities. In fig. 12 we have plotted 4k2Dz (where + is the duration to 
exp ots I) times maximum amplitude) against «. For a given diffusion 
coefficient, the duration is almost independent of line density up to «=10" 
and then increases approximately proportionately to «. This agrees well 
with the experimental observation that at a given height most echoes 
have about the same duration except for a few which persist for a longer 
time (up to several hundred times as long). Wéith the previous a Ne 
we find 70-03 sec for «<10¥ and +r=30 sec for «=105, ; 
In fig. 12 we have also plotted the maximum reflection coefficient against 
a. It increases approximately linearly wp to c<=10!2 and then much more 
slowly (=o). This again is in agreement with observation. 
The ideal picture of a uniform cylindrical column of ionization fails to 
take account of any distortion of the trail which may occur through 
atmospheric turbulence. Such effects have been considered by Greenhow 


Oo 5 10" 2 5 10% 2 5 jo" 2 Ly 1o!4 ra 5 rae 
Parallel scattering. 
----- Transverse scattering. 
Echo amplitude reflection coefficient g, and duration 7, vs electron line density. 


(1950) as an explanation of the amplitude fluctuations which occur with 
almost all long duration echoes. The theory given here supports the view 
of Greenhow against that of other authors (notably Feinstein 1951) who 
suggest that fluctuations may derive from interference between waves 
reflected at different levels in an expanding column. For this purpose it is 
necessary to postulate a fairly sharp boundary to the cylinder, since it is 
clear from (61) that no such effect exists for a Gaussian electron distribu- 
tion. Now the long duration echoes, on the above theory (and irrespective 
of the radial electron distribution), correspond to columns which are so 
dense that they behave like an expanding metallic cylinder and. it is 
clear from fig. 3 that they will not exhibit fluctuations (at least during the- 
main body of the echo), unless disturbed in some manner. Only short. 
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duration echoes (r~0-03 sec with the previous parameters) could exhibit 

fluctuations of the type suggested by Feinstein and the fact that they are- 
not observed favours a Gaussian, or near Gaussian, distribution. Both 

atmospheric turbulence and loss of electrons by attachment (which may be 

important at the lower heights) may reduce the amplitudes and durations 

of very dense columns, particularly on long wavelengths. 

Summarizing the above, we may say that a qualitative change occurs in 
the character of the echoes at line densities in the vicinity of a few times 
10!2 per cm, which result, surprisingly enough, is independent of both 
wavelength and diffusion coefficient. When « is less than 101°, the echo 
amplitude for parallel scattering obeys the Lovell-Clegg law with an 
exponential decay, while that for transverse scattering exhibits a resonance 
phenomenon when the axial electron density corresponds to a dielectric 
constant of —1-4 (for a Gaussian distribution). In addition, the echo’ 
duration is independent of «, being a function only of A and D, while the 
maximum echo amplitude is proportional to «. These are the typical 
short duration’ echoes which form the bulk of the experimental obser- 
- vations. When «>10!2, the meteor trail reflects as a metallic cylinder 
and, while (kr,)? <1, the echo amplitude for parallel scattering exceeds that 
for transverse. The duration is approximately proportional to «, while the 
maximum amplitude is only a slowly increasing function of «. These are 
the columns which give rise to the long duration echoes (up to several 
minutes) studied by Greenhow (1950). An important result of this is that 
it enables us to determine « (when « >10!°) independently of the equipment 
parameters, provided that D is known. 

The method developed for the solution of the wave equation, namely the 
reduction of the problem to a formal electrostatic one, while restricted as 
to the dimensions of the scatterer, nevertheless has a certain generality. 
The solution may be obtained for any given dielectric constant as a function 
of radius and may be extended to geometrical shapes other than the circular 
cylinder, for instance the diffuse sphere of Appendix II. In a later paper 
the same general technique will be applied to the problem of scattering 
from a thin diffuse column of elliptical cross-section. 

The theory is also of importance in connection with laboratory experi- 
ments with discharge tubes such as have been made by Romell (1951) and 
by Denno, Prime and Craggs (1950), both in the analysis of results and in 
suggesting new experiments. For instance Appendix II enables account 
to be taken of the walls of the discharge tube, and also suggests that it may 
be of interest to study a discharge confined between concentric tubes. 
In this latter case we can arrange for the dielectric constant at resonance to 
range between 0 and — oo and in principle have a method of verifying the 
relation (2) over this range. This may make possible, in particular, an 
ees. a oe Eh Se of the Lorentz term for a 
ideo er picamoreiabies : é 4 roviding that collision damping is 
Hehenes eee oval of the mode degeneracy should enable the 

sonant modes to be detected. 
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APPENDIX I 


EVALUATION OF THE REFLECTION COEFFICIENTS FOR AN ELECTRON 
DENSITY PROPORTIONAL TO exp [—(7/7)*] 


This defines a cylindrical column of electrons, the diffuseness of which 
depends on s, and leads to the relation (54) for x(r). Writing p=r/7o, (54) 
becomes 

(py Lk = frexp ("pe tease # +. > (69) 
We are concerned only with the dipole mode, hence the differential 
equation to be solved is, from (7) 
d dV KV ie 
oe ier) a ee ee, 27) 
dp (° “dp ) pa 
where the mode subscript has been dropped. 


External to the column, «=1, and we may write V=A'p+ B’p~', where 
A’ and B’ are the. parameters in the equation (63) for the transverse 
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reflection coefficient ((63), although derived for a Gaussian electron 
distribution, is valid for all s). 5 
From (25), the parallel scattering reflection coefficient is 


pe a(2/s) ! 
iy flbra)®[ exp (att) de= "EP fbr), (T) 
where x=p?. Using (63) and (71), the ratio of reflection coefficients when. 


the resonance is diffuseness limited is 


ok 
m  (2/s)'f 
reducing to (64) when s=2. 

To determine A’ and B’, V must be expressed as a function of p. When 
p<1, x = const. and the solution of (70) in the vicinity of the origin may 
be taken as V = p. Hence, when p=0 we have V=0, V’=1, «V/p=x(0). 
Integration of (70) thus gives 


Tey 


\ 
a (72 
A! \ ) 


? 


. (73) 
Ve | V'dp, | 
0 
where the prime here refers to differentiation with respect to p. 
By commencing at the origin and taking successive small increments in p 
we can, with the help of (73), evaluate V, point by point, using a method 
of successive approximations. Suppose V and V’ are known for p=p,, 


then they may be evaluated when p=p,—p,+6p as follows. Weassumea 
value for (x V/p), whence 


) 
(pe Vo = (pV), + [eV p): +(«V/p)a 
==(pKV’), bole VO Fee ee eee 


(where subscripts 1 and 2 refer to values at p, and p, respectively). 
From (74) we find V’(p,) and hence 


‘ ) 
V(p2) = V(ps)-+ $ [V'(ex) + V'(p2)] 
=V(p +5). 80, ee ee ee ey 


From (75) we obtain a second approximation to («V/p), and we repeat the 
sequence until sufficient accuracy is achieved. By carefully choosing the 
first approximation we never need, in practice, go beyond the second. The 
only error remaining is that involved in the approximate integrations (74) 
and (75), which may be made as small as we please by taking a sufficiently 
small increment, dp. 

The procedure is illustrated in table 1 for the example s=4, f=2. 
Table 1 (a) is restricted to the range p<0-9124; i.e. «<0, since when x0 
both V and V’ go to infinity, while px V’ passes through a broad minimum. 
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In order to bypass the singularity we use (47) which becomes 


rete). 
=VyrtJjV ar- Marr Airs are es 
In the example we put 
p,=0-88, «,=—0-0980, Vy=1-337, (px V’),-9=— 9-715, 
and the real and imaginary parts of V and V’, in the vicinity of the 
singularity become 


Vy 0-734—0-2595 In | «| | i: 
Vi & —0-715/(px) 
V; = +0-815,, 9) Vy Ce ae 


We evaluate V, and V; separately. Using (77) and (78) we obtain the 
values in the first rows of tables 1 (b) and 1 (c) respectively, and then pro- 
ceed as before until « is effectively unity. 

External to the column, Vp=Ap’p+Bp’p7!, and, since «=1, 


2An’=KVylp+Ve',  2By'=p(Vn—peVn'), « - (79) 


and similarly for A,’ and B,’, where A’=A,y’+JjAy;’, B’=B,p’+ )By'- 
From the last rows of tables 1 (b) and I (c), we get Ap’ =0-073, Bp’ =0-866, 
A;’=— 0-376, B,;'=—0-461, hence |B’/A’|=2-56. Using (72) with 
(2) !=/n/2 (for s=4), we obtain g,/g,=2-89. 

The same procedure has been carried through for s=2, and a wide range 
of values of f, leading to the values of Ap’/Bp’, Ay’ /Bp’, By'/Bp’, plotted in 
fig. 13. When f decreases towards unity, both A,’ and B,’ tend to infinity, 
however, it is not difficult to show that A;’/B;’/=A,’/B,’ at f=1. When 
f<1, of course, both A,’ and B;’ are zero. 

For other s, the solutions have been evaluated in the vicinity of the 
resonance and fig. 14 gives the values of A,'/Bp’, By'/B,’, for the resonant 
condition (A,’=0), as a function of s. , 

It is clear that the above method of numerical integration is generally 
applicable to the solution of (7) for all modes and for any electron density 
distribution which is a function only of radius. 


ACPPo WONG 


Some ParricuLarR ExLecrron DistrRipuTIONS EXHIBITING MULTIPOLE 
RESONANCES 


(a) Stepped Cylindrical Column 


We will put «=«,, when 0<r<a and «=x, when a<r<b. 
The solution of (7), external to the column, is Vin=Ani™+B,70-", 
where ' 


An l bP" (K+ )(1 + eg) +0?" — 1 )(1—K) 


Be ~ Bam B?™ (9+ y)(1— Kg) +a2™(Ky>—K)(1 + kg) : (80) 


' 
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Eqn. (20) then gives the reflection coefficient for the various modes, and, 
if both x, and «, are real, resonancé occurs when the numerator of (80) is 
zero. 

This result may be useful in connection with laboratory experiments 
using discharge tubes since it enables account to be taken of the glass walls 


Fig. 13 


-0-5 


(e) 2 4 6 8 10 (2 14 16 18 20 


by making «, the dielectric constant of the glass, a and b eee and gs 
dii of the tube, and x, the dielectric constant within the discharge. 

" f the glass walls has the effect of removing the degeneracy of the 

es a however the splitting would seem too small and in the 

See ne to explain the subsidiary resonances observed by Romell 

wr 10N1, 


(1950). 


30 T. R. Kaiser and R. L. Closs on the 


The case of «,=1, c2=« <1, corresponding to a plasma confined between. 
concentric tubes, is of some interest. We obtain 
A 1 b2"(1-+-«)2—a2"(1—x)? 


ib se 1 
Be h2m (62m —q@2m)(1 — x?) 2 | (8 ) 


and resonance occurs when c=«, and «=, where 


Ka=Kp l= (a™—b™)/(a™-++6"), nay Sn a 
_ ie. there are two resonances for each mode, 0>«x,>—1, and —1>«,>— 0 
and as m-—> o they crowd towards x=—1. The degeneracy is again 


removed and, when b/a approaches unity, considerable separation of the 
first few resonant modes may be expected. The values of « at resonance 
are given, for the first few modes, in fig. 15, for ba=1-5. 


Fig. 15 


{ 2 Se 4p See Cas, 
O—Multipole resonance conditions for an ionized annular cylinder having 
external and internal radii in the ratio 3/2. 


@—Multipole resonance conditions for a diffuse spherical distribution of 
ionization. 


The above result suggests an experiment (see § 7) to determine the 
magnitude of the disputed Lorentz term in the expression for x, namely 
1 4rne? 

k?mce?+-p47ne?" 
The factor p, which we have taken as zero, is given the value 1/3 by some 
authors. From (83), we see that when n > 0, « > 1—p-, ie. if p= 1/3 
k*—2. If pis finite then, as b/a is increased towards unity, the TesOnAanos 
defined by i, in (82) will, at some point, disappear, even for the highest 
electron densities. 


(6) Spherical Electron Cloud 


k= 


(83) 


An approximate solution of the scattering of an electromagnetic wave 
from a small ((kr)?<1) spherical cloud of electrons may be obtained by a 


similar technique to that used above for evli ill gi 
it used ¢ ylinders. We will 
outline of the method. vibes 


ra 
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We consider, first, a sphere of dielectric constant, «(r), centred at ‘the 
origin of spherical coordinates. If the incident wave is propagated along 


the z-axis and polarized parallel to the z-axis, the appropriate potential 
function is 


V=3V,, CONG Ls (COs 0) aa ta an merle tn (84) 

where P,, is an associate Legendre Function and V,, 18 the solution of 
Cal werd. Vi 

alee = Meal )ieViow,  & + a1 (SO) 


These correspond to the eqns. (8) and (7) for the cylinder and (85) may 
be integrated numerically for any given «(7) in a similar manner to (7) (see 
Appendix I). 

External to the sphere, «1, and we find 


=A Be . . . . . . . (86) 


The exact form of the solution external to the sphere is more difficult to 
obtain than before since the phase of the incident wave is not independent 
of any coordinate and we need to express a unit vector plane wave as a set 
of spherical waves. The total electric field components external to the 
sphere are given by Stratton (1941). Replacing the spherical Bessel 
functions in Stratton’s expressions by their approximate values when 
(kr)? <1 and comparing the result with the field components derived from 
(86), we obtain the following expression for the reflection coefficient, a,,,. 


_ m(2m)! (2m+1)! A 


a Se Mele Sth 2: 7 
i Wel 4m(m-- 1)(m Ng tla SB oe (8 ) 
The total reflection coefficient for backward scattering is now 
Gael (Mies) Omi cae eee es *) (88) 
1 


and the backward scattering cross-section is 
Os 0 2N Tae eo RE ea OO) 
For a diffuse spherical cloud, using the Herlofson model, we obtain. 
for a sphere of radius a, 


A m+ m[x«/(ax’)\(gr7—In | « |) +[1+m«/(m+1)] 


m NS Oo 90 
B,, ma mix [(ax)\ gn —ln |x |) FO —*) 20) 
where «’ is the gradient of « at c=0. 
For a homogeneous sphere, «’ > 00, and 
, An  kK+-(m+1)/m (91) 
Be — arA(1— Ke) . . . . . . . 
Resonance occurs when «=—(m-+ 1)/m, which is given, for the first few 


modes, in fig. 15. 
We will define the magnitude of the resonance as the ratio of the resonant 


reflection coefficient to that of the same number of free electrons oscillating 
in the same phase ; it is thus equivalent to g,/g, for a resonant cylinder. 
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Pa 


Denoting it by 7 and proceeding as in §5 we get, for the diffuseness limited 


dipole resonance, 
3 
asi Sie hath, CP? 


~~ L |) 
1=3/90na—jn) 1 |” 


where the electron density, as before, is proportional to exp [—(7/79)*]. 
We have plotted 7 against s in fig. 16 and although the values derived from 
the Herlofson model may be somewhat less than the true ones (as in the 
cylindrical case), it is clear that the diffuseness of the boundary has a 
greater effect for the sphere than for the cylinder. 


Fig. 16 
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Magnitude of the dipole resonance for a diffuse spherical distribution 
of ionization. 
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II. An Equation for a Particle with Two Mass States and Positive 
Charge Density 


By H. J. Buasna, F.BS. 
Tata Institute of Fundamental Research, Bombay, India * 


[Received September 21, 1951] 


SUMMARY 


A relativistic wave equation in the canonical form (2) is given which 
describes a particle having two different rest masses. In one mass state 
the particle has a spin 3/2, in the other a spin 1/2. The free charge density 
is positive definite so that the equation can be quantized in the usual 
manner in accordance with the exclusion principle. The equation is 
irreducible. 

This example establishes the principle that irreducible multi-mass 
equations exist which can be used to form the basis of a physical theory 
along the usual lines, contrary to the general prevalent view. 


§1. IyTRODUCTION 
THE purpose of this paper is to give an irreducible relativistic wave- 
equation in which the particle has states corresponding to two different 
values for its rest mass. Unlike all previously known equations of this 
type, however, the free charge density is positive definite, as for the 
Dirac equation, so that the equation can be quantized in accordance with 
the exclusion principle in the usual way. The equation is derivable by 
the variation principle from a Lagrangian density of the canonical form 


Bt Dap, +xyps OOS BO he IO ani 5 (1) 
wheret p,——i0/dx* and x is a constant of the same dimension as #,. 
The matrices D and «” have certain standard properties which will be 
stated in the next section. The wave equation has therefore the 
canonical form : | 
(opr i Oe eee rahe Nes so, Da (2) 
where the matrices «* are defined in the next section. o 

Equations of the form (2) which have solutions for several different 
values of the invariant length p*p, of the energy momentum vector p, 
have been known for some time (cf. Bhabha 1945). But it has turned out 
that for all these equations neither the total charge nor the total energy is 
positive definite, so that it has not been possible to interpret the equations 


* Communicated by the Author. Px 
hie, oe are the coordinates and the metric tensor g,,, 18 taken in the form 


Jo0o= —I11 = 922 Isa—l, I= 9 for kAl. | 
SER. 7, VOL. 43, NO. 330.—JAN. 1952 D 
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in the usual manner as describing particles with states of different rest 
mass. This has led to the general belief that irreducible equations 
describing a particle with several values of the rest mass cannot exist. 
Pais and Uhlenbeck (1950), basing their investigations on a much more 
limited form of the Lagrange function than (1), have shown that a 
physically sensible theory cannot be made for multi-mass equations due to 
the lack of positive definiteness of the energy and the charge. The 
present paper therefore proves that the generally held belief mentioned 
above is false, and that physically sensible theories can be constructed on 
the basis of multi-mass equations provided one works within the more 
general framework of Lagrange functions of the canonical form (1). 

Whether the equation given in this paper will prove to describe an 
elementary particle actually found in Nature is not a question which can be 
answered now. The main importance of the present investigation in my 
opinion is that it establishes the principle that certainly one, and probably 
several, irreducible relativistic wave equations exist upon which can be 
based a physically sensible theory for particles having states corresponding 
to more than one value of the rest mass. 

The new equation is defined in §2. Its general properties are established 
in §3. Its irreducibility is proved in §4. An explicit representation of the 
matrices «* is given in the Appendix. 


§2. DEFINITION OF THE EQUATION 


The matrix D is by definition non-singular and hermitian. It is such 
that the four matrices Da" are likewise hermitian. 

The wave function % transforms according to a representation # of the 
Lorentz group, so that a transformation t of the coordinates : 


hth a . . . . . . . . . (3) 
causes a transformation of the wave-function given by 
pT, oo Ce ee eS 


where 7’ is the matrix representing tin Z. The invariance of the expres- 
sion (1) then requires the five matrices D and «* to transform according to 


l™DT=D, ie nica ers Ma bolts aes Da 
Tg Tt cl ee ee (4 c) 
The components of %, and in consequence the matrices of the represent- 


ation #, and the matrices D and «* can always be subjected to simul- 
taneous transformation by an arbitrary non-singular matrix s : 


p> bate y Os ee es 
D + D'=s'Ds, Se een ete ee OY 
(ide A te CPi Re Me eRe A ec 
PAT aga Rae i Me meas (BAS 


without changing anything in the theory. All such formulations are 
equivalent. 
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We take a % transforming according to the representation 


where A(z, 4) and &(4, 4) denote two irreducible representations of the 


full Lorentz group in the usual notation. Under restriction to the proper 
Lorentz group the representation A(k+l, k—l), k>1 breaks up into two 
irreducible representations G(k, 1) and F(I, k) of the proper Lorentz group. 
Thus, 2(3, )~B(1, $)+A(k, 1), while, Z(4, 4)~D(4, 0)+H(0, 4). The 


: : 22 2 
representation Y(k, 1) describes the transformation properties of a spinor 
a; Fise symmetric in its 2k undotted indices and likewise symmetric in its 


21 dotted ones. We assume that & is already in its fully reduced form 
with respect to the proper Lorentz group, and that its irreducible parts are 
arranged in the order 


A~H1, 3)+ AO, 3)+H(0, 3) + HE, I+H(S, 0)+A(4, 0). 


Then 

$a 

bi? 
| p@ 
= ; ry 7 dee Siesta te ened OLE 
ames (7) 
| pide | 

UO ea a) 

Moreover, we arrange the six independent components of the first spinor 
in (7) in a column in the order ¢;", /24;™, 4;, $3"", 1/263", $3”. Simi- 
larly, for the spinor ¢;;° we take the components in the order ¢j;', 1/24)". 
é33', ete. The second, third, fifth and sixth spinors have only two 
components each, and we take them in the order 4;", 45", 4i°, 63, 
; f1, po, etc. 

In a recent paper* (Bhabha 1951 a) I have investigated a class of 
Lagrange functions of the type (1) in which } transforms according to the 
representation &(3, $)+(4, 4). We have only to repeat the methods and 
arguments of that paper to arrive at the general form of the matrices 
D and « in the present case. We first divide up the matrices D and ov 
into submatrices corresponding to the subdivision of % into irreducible 
spinors as in (7). We then arrive at the result that the most general form 
of D, which is hermitian and satisfies (4 6), is 


ely ak i) EEA ea EN} 
where 
dl) ; 
D= . dyh (a) A . . PS e . e (8 b) 
Alia) } 


* Referred to here as (A). 
D2 
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1,,, stands for the unit matrix of n rows or columns, and d,, dz, and dz 
stand for either.1 or —1 independently of each other. 

The most general form of the four matrices a” satisfying (4c) and the 
condition that the matrices Da shall all be hermitian is then 


a’? o* 5 
oo— ( 10 : ) 5 i ( ; tke ) > for — Le 2: ay . (9 a) 
a . —a . 


where A 
f ayuk d,cgv’ —-d¢gv" 
ah : d.c,w* Geen CH ake idee slab: ib. 
|. dacqu® =F dgtya¥ i” | atge™ 


Gy, Ao, Az, C1, Cy and Cy are six arbitrary numbers of which a, a2, a; must be 
real, while c. and c, can be taken to be real without loss of generality. 
For if they are not real, they can be made so by a similarity transformation — 
of the type (5). The sets of matrices uv", v*, w* and z* are reduction matrices ° 
whose form is fixed, except for multiplication by arbitrary numerical 
factors, once the form of ¢ is fixed as in (7). 

The general form of the matrices D and «* given by (8) and (9) is easily 
understood. A reflection turns an upper undotted spinor index into a 
lower dotted one, and vice versa. With the arrangement of the spinors 
in >; given by (7) the transformation r which reverses the directions of the 
three space axes must be represented therefore by a matrix # having the 
form (8). By (46) the matrix D must commute with R since R'=R=R!. 
Moreover, taking the complex conjugate of a spinor turns an undotted 
index into a dotted one and vice versa. D must therefore have the form 
(8) with d,, dy, d; any real numbers.’ By a transformation of the type (5) 
we can then make the modules of these numbers 1. This settles the form 
of D given by (8). 

Denoting by (k’, l’| «| k, 1) the submatrix of « which lies in the row 
corresponding to the representation D(k’, l’) and the column corresponding 
to the representation Y(k, l), I have shown (Bhabha 1945, §5) that 


k’=k+4 
U=144 


Further, according to (4c), R commutes with « and anticommutes with 
a*, k=1, 2,3. This settles the form (9) of the «*. Finally, the require- 
ment that the four matrices Da" shall be hermitian requires that a, wu’, 
a" and a;z" shall be hermitian. Without loss of generality we can 
therefore take u* and z* to be hermitian, and a,, a, and a, to be real 
numbers. Finally, we must have 


(k’, l'| «| k, )=0 unless independently. — 20) 


UE We ee ee Ce. CAD 


The actual form of the matrices uw", v* and z* can be obtained in terms of 


the reduction matrices w', v, u and v" introduced by Dirac (1936) and 
developed by Fierz (1939). An irreducible representation Y(k, 1) of the 
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proper Lorentz group is determined by its six infinitesimal transformations 


i; (k, 1). This set can be connected with two spinors K*" and Lt sym- 
metric in their spinor indices through the relations 


K(k) =—K3(k) =I +7122, 
KM= O46), b. .. (12) 
Ki(k)= ! iI?) 1 4( 1 L 43), 


and 


; 
DES Ney Bg) cea Be | 


L3\(l)= (I 113) + (7? iI), b ee ial heeds 


L2()=(! Wis) (1 ED | 


Then the two spinor matrices u(k) having 2k-+1 rows and 2k columns and 
the two spinor matrices v“(k) having 2k rows and 2k+1 columns satisfy, 
and are defined by, the equations* 


—u,(kb-+4)u"(k+4)=0,(k)yul(k)=2k+1 
v,(k)o"(k+4)=u,(k+ 4)u"(k)=0 | 
v'(k+4h)u,(k+4)=K"(k)+(k+1)8," (14) 
ae ul'(k)v,(k)=K “(k)—ks 


ul and v" satisfy similar equations with L;! in place of K,". One sees. 
at once from (9) and (7) that u*, vo” and z* are just the submatrices 


u*=(1, $| a*| A, AN) | 
REL AL ok (0am ot me anaes, bine (15) 
Af, 0] ot} 0,3) | 


Hence it follows from the form for these submatrices given in the paper 
quoted above (Bhabha 1945, formulae (66) and (67)), that we can take 


uk =o,,ku*(1)v"(1) 


| 
. | 
vb =o, Fu* Lyul'(d) [ Peete tat p tor LG) 
| 
J 


where o’ stands for the 2X 2 unit matrix, and o', o”, o® for the three Pauli 
matrices. The dotted index labels the rows and the undotted one the 
columns of these matrices. An explicit representation for the matrices 
u®, v* and z* is given in the Appendix. 
i a ee 

, * cf. Bhabha (1945), §5. 
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We now maintain that the Lagrange function (1), with the matrices D and 
a* given by (8) and (9), and the constants d, a and c defined by 


dd, a= 
Ges, a, >4, as=4—a,—|A|, . . . (176) 
1 
tas 5 See Ag—% Wale Dep OW A aie | lve 
c;=0, C2 =343 ee b) C3 =3 4 Ag ? 7 ( ) 


describes a particle having states of rest mass x and x/|A| and a positive 
definite free charge density. Its spin in the state of rest mass x is 3/2 in the 
state of rest mass x/|A| it is 1/2. The equation is wreducible. 

Provided a,40 we can always put a,=1 without any specialization 
since this merely changes the definition of the mass in terms of xy. More- 
over, the hermitian matrix 

A, Cy 
(*: 2) ; 


which multiplies z" in (9 b) when d,=d,;= +1 can always be brought to the 
diagonal form through a similarity transformation by a unitary matrix 
without changing the corresponding form of D’. Hence putting c,=0 in 
(9 6) also represents no specialization when d,=d;. Finally, we note that 
c, and c; can be assumed to be purely real without any loss of generality. 
The sign of either is irrelevant since it can be changed by a transform- 
ation of the type (5). Thus (17 c) only defines the absolute magnitude of 
Cc, and cs. 

Introducing (16) into (9 6) and using (17) it can be shown without much 
difficulty that equation (2) is equivalent to the set of spinor equations 


; Oo of . C oO Oo oO 
HOM pai +P bsH)— 5 75 (Pa BO” + Dab) 


C. at . é 
— 375 (Bed + PPO) +xG,"=0, (18.4) 
= Fao Fi + ta pepe + xp=0, Be 
c . me - 
— Py Pe bi + asp P+ xp —=0, . . (18 0) 


Io o ¢ oO o 
MPioba + Piobi")— 375 (Debs + Pi bs”) 


Cy 2 o1.(2 

— 525 (Pi + 09b))+rbsi=0, (18 d) 
eins 

a Pe bai + aePieh + xb,=0, —. (18 8) 


Cc . ; 
ig Pe bit + ag Pigh + x$=0. . (18f) 
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§3. PROPERTIES OF THE EQUATION 
By A, we denote in general a matrix having every element on the 
principal diagonal equal to A,, any element on the diagonal just above the 
principal diagonal equal to 0 or 1 independently of the others, and every 
other element zero. It is well known that every matrix, and in particular 
z, can be brought to the pseudo-diagonal form 


AS, 


Zip 
2; 


at eee ie (19) 


aa 
i 
of 


where A,4A,~... AA, are the eigenvalues 
We now note the following trivial : 
Lemma 1. H and A being two matrices such that H and HA are both 
hermitian, all the submatrices (i | H | j) are zero for 14) of A is of the form (19) 
and all tts ergenvalues are real. 
(i | H | 7) denotes the submatrix of H which lies in the row corresponding 
to A; and the column corresponding to A;. For the hermiticity of HA 
requires that 


ot 


| {@| 1] j)Aj}'S(9 | 411) A, 
Since (| H|j)'=(j| H|2), this requires 
ANGLE )=G HOA, - 
from which follows the statement of the lemma because = AFA. 

Lemma 2. If the submatrix A, in (19) be a multiple i, of the unit matrix, 
then the corresponding submatrix (j | H| 7) of H can be brought to the diagonal 
form without changing the form of Aj. / 

For (j| H|j) being hermitian, it can be brought to the diagonal form 
by a unitary matrix wu: uj] A | jyu=u'(9 | | ju is diagonal. But 
u—*Au=A, since A, is a multiple of the unit matrix. 

If ‘A’ be a matrix of » rows or columns having all elements on the 
principal diagonal equal to A, all elements on the diagonal above the 
principal diagonal equal to 1, and all other elements 0, then we can prove : 

Lemma3. Hand A’ being two matrices such that H and HA’ are hermatian 
and Xa real number, H must have the form 


: h, 
beens 
Pe tow ae 
' (20) 
| hy, he : . 0 © hy—2 hy—4 
{ h, hy hs : : . ° hy_4 h, J 


where hy, hg, »- +» hy, are real numbers. 
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All elements on the principal transverse diagonal are h,, all elements om 
the transverse diagonal below this hj, and soon. For from the hermiticity 
of HA’ follows 

(2) Z| Wry=C| A] +0] A] 2), 
whence 
MIA] 2)=(. |Z] +0] Z| 2), 
or (1|H|1)=0. Proceeding in this way we can prove that (1 | A |= 
for r<n. Hence (r|H|1) is likewise zero for r<n. Taking the next 
column we can then prove that (2 | H | r)=(r | A | 2)=0 for r<n—1, while 
(n—1|H|2)=(n|H|1)=h, say. We eventually arrive at (20). The 


hermiticity of H then requires h,, hz, .. . h,, to be real numbers. 
Consider now the matrix equation 
(Ap, 4-00 ee eee eee 
where p, and y are numbers and y=0. % is a column matrix with the 
elements 1, 2, ... %,. This equation is equivalent to the set 


fF — (Apo+x)b1 +Pof2=0 
(Apo+x)bo+ pos3=90 : 
cia ake: (21 b) 


(Apo+x)bn—1t+Potn=9 
L (Apo +x)bn =0 


If Ap) +x <0 the only solution of this set is J=0. If Apo+x=—0 there is 
only one non-vanishing solution, namely 


%, arbitrary, pop. == 0. Slee 
If A=0, then Apy+x cannot be zero, and there is no solution of the type 
(22). 

We now return to eqn. (2). It has been proved generally (cf. 
Bhabha 1949) that the rest masses of the particle are connected with the 
eigenvalues of «° if y40. Indeed, when y0 one can define a value of 
the rest mass as the value of |p, | for which a non-vanishing solution of 


(p+ x)p=0 fh Kee wont 


exists. Our first step is therefore to bring «° to its diagonal or pseudo- 
diagonal form (19), from which all the mass states of the particle can be 
determined. 

In all equations derivable from the Lagrange function (1), the charge- 
current density vector is ¢'Da*. The charge density is therefore 


ib! Dow? eta aah msg cnet mes cae ee 


In the Dirac equation Dz” is a positive definite matrix. This requirement 
is however too stringent, and it is only necessary for a sensible physical 
interpretation that the expression (24) shall be positive for every free wave 
solution of (2). When this is so, we say that the free charge density is 
positive definite. This will be so from relativistic invariance if it is so for . 
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every solution of (23) corresponding to a particle at rest. We have there- 
fore to follow through the transformation of D under a transformation of 
the type (5 b) when «° undergoes a transformation of the type (5c). 

We first note that a transformation of the type (5) with an s given by 


1 —l il ) 
pM So (10) 70) } Ns be eas ee ent 5: 
val Lio) Lao) ee 
immediately changes D and the «* to 

eas)" 4 ce 70 
bal a). #e(E a) 
bo ( ee) k=1, 2,3 266 
Pees \ ial oe ten ar (ea) 


Introducing the representation for the matrices wu, v° and z° given in 
the Appendix into (9 6), and giving the constants the values (17), we get 


/ 
oO 


ESE tees 


I v2 203 * 302 
i) 1 1 
" v2 ve V2 
1 il 1 
5D} 4/2 “3 4/2 2 
1 
V2 t — 303 — $C 
1 
4 tl 
— C3 V2 Cz As 
I 
2 C3 303 oe 
il 
— $C: Fy ae Cy . . C s a3 
bes 5 
il 
73 : dey . | : 4 . a3 j 
| (27) 


This matrix becomes precisely the matrix dealt with in the preceding 
paper ((A) formula (26)) if we omit its last two rows and. columns. We can 
therefore deal with it exactly like the matrix in that paper’. By a simple 


on ee Le oe en te ae 
* cf. also Bhabha 1951b, Report of an International Conference on 
Elementary Particles, Bombay, 1950. 
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re-arrangement of the rows and columns of (27) we can bring it to the 
form 


2 7 23 22 
\ : as 
lio) - s — 72 %8 ere 
ce a” where «” = 1 3 
on” — 303 a2 Ws. = Os 
—te, as ee As 
Wek S 
(28) - 
while D’ is brought simultaneously to the form 
1 , 
lip) ; ; ; ‘ 
D'>} . Dp’ : where D’= | #3 
edt pie ae 
(29) 


It is clear that as in (A), the matrix «’° must have the eigenvalue 1 four 
times, for the component ¢;"' of the first spinor in (7), which is a component 
of the wave-function belonging to the eigenvalue 1 of «’°, transforms 
according to the representation ZY, of the three dimensional rotation group 
which is of degree 4. Indeed, a transformation of D” and «” of the type 
(5b) and (5c) by a matrix 
| V3 Vocal 
vt —val. + 


ens ee ee 
1 


1 
immediately diagonalizes the eigenvalue 1 of «’° without changing the 
form of D”. Hence D’ and «® are transformed simultaneously to 


1¢4) . . Lig) 
p me . | i ? aie ince . C . > 
7 C 
31 
where (30) 
V/3 /3 
1 3 Daan choi e 
/3 
US —1l e ’ C= ha = Cy A> 
—] 
/3 
ayes : as 


(32) 
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The matrix ¢ has the characteristic equation 
[A—€([={A—(@,+a,—3)2=0. 1°. . . (33) 

Hence it has the two roots 
A=a,+a,—} . . . . . . . . (34) 
and zero. Moreover zero is a double root since € does not. satisfy the 


equation {A—(a.+a,;—4)!A=0. Hence the pseudo-diagonal form of ¢ 
is 


oe aa 
ie ee oe eS Wie eM a 3p) 


When ¢ is transformed to the form A by a matrix s: 


Saco 7/1, har ee bee eee ane (OO) 
» is transformed to 
eS 1S Merge se, Alb etrhs 1.(37) 
and it follows from lemma 1 and 3 that 7’ must have the form 
| ™ : : ; 
ey : Ho Pe ere eee (os) 


| No UE 


It follows from the statement made immediately after (22) that the 
equation 
(Apo+x)$b=0 
has only one non-vanishing solution, namely 
7, arbitrary, oe 0. 
and this is possible only if 


Dee ae I est) iy (39) 


7 || 
For this solution the charge density (24) is 
hyn =n Abs Mere oe hte (40) 
It is therefore only necessary to calculate 7,’. 


Denoting by s, the first column of the matrices in (36), it follows from 


(32) and (37) that 4 


Tiga ROE Behe aes et. . (41) 
=| 811? —[s21 —|Ssi| 
where 51, S:; and ss, are the three components of s,. But by (36) s, 
satisfies 
63, = 8A. 
From this it follows that but for an arbitrary multiplying factor 
1 1 f 
2 I Sa= a Se ° eee ee ee ; (42) 


ee : 23) aa ee a 
iis Bae Co(4s—5) C3(4.—3) 
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Substituting this in (41) and inserting the values (17 ¢) for c, and Cs we get 


3 (d@,—3)(da+a3—4)" 
es coli gerin tat Sree, Lia et EN hE fh Mae er a (sy) 
™ =~ 4 a2a,%(ag—3)(ds—3 oy 


a, and dz being real, 7," is always negative. Hence, if a, is so chosen that 
A<0, the charge density (40) for the state of rest mass x/|A| is positive. 

The equation also has a mass state x which arises from the eigenvalue 1 
of «’° in (31). The charge density is clearly seen from (31) to be positive 
for this state. 

We have therefore proved that the equation (2) determined by the 
values (17) of the constants has two mass states x and x/|A|, and a 
positive definite free charge density. 

The spin in the state of mass x is 3/2. This follows from the fact that 
one of the components of % corresponding to this value of the mass is 4;"", 
which transforms according to the representation Z, of the three dimen- 
sional rotation group. Since the eigenvalue 1 only occurs four times in 
(31), there can be no other spin value for this mass state. 

The eigenvalue A occurs only once in ¢, but ¢ occurs twice in «’® as seen 
from (31). Hence in the state of rest mass x/|A| the spin must be 1/2, 
for the representation Y, of the rotation group has just two components. 
(Spin zero is: excluded by the fact that we are dealing exclusively with 
spinors having one or three spinor indices.) 


$4. IRREDUCIBILITY OF THE EQUATION 


We shall now prove that the four matrices form an indecomposable set, 
. that is, they cannot be decomposed into a direct sum of such sets. Suppose 
this is not so, and let the set be decomposable into the sum of a number of 
indecomposable sets 

oF ag DRL G(R Te 5) te eet ie vera 
In the rest of this section we simply drop the index k for conveniénce. 
Now we have to consider cases of two types. Hither the representation Z 
is not correspondingly decomposable, or it is already correspondingly 
decomposed. 

A case of the former type was investigated in (A). There the «-matrices 
decomposed into a direct sum of four sets of Dirac matrices although 
%~ transformed according to representation #(3,4)+A(4, 4). This is 
because a y transforming in this way can always be written as the direct 
product of a Dirac wave-function times a four vector. In the present 
case a % transforming according to the representation (6) can be written 
as the direct product of a Dirac wave-function times a four vector plus a 
scalar. Thus, if a contingency of the first type arises either the «’s 
decompose into a direct sum of five sets of Dirac «-matrices or a direct 
sum of four sets of B-matrices, where by B-matrices I mean the set of four 
matrices given by Kemmer which occurs in the scalar meson wave-equation. 
In the former case the eigenvalues of «° must all be +1, in the latter +1 
and zero. Since «° determined by (17) does not possess these eigenvalues 
we have proved that the first type of decomposition cannot occur. 
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We have therefore to consider the case in which the «’s decompose as in 
(44) and the representation & decomposes correspondingly. Now either 
the part of # corresponding to one of the sets say, «, in (44) only con- 
tains the representation 4&(3,4), or it contains the representation 
A(3, 4)+A(4, 4). It follows from the investigation in (A) that the 
former alternative is excluded since otherwise the eigenvalue of the 
corresponding matrix «~” would be different from those of «° given by (17), 
and the free charge density would also not be positive definite. The second 
possibility needs further consideration, for it includes the case in which our 
equation might decompose into a direct sum of the canonical form of the 
Dirac—Fierz—Pauli equation for a particle of spin 3/2 given in (A), and the 
Dirac equation. Any other possibility is excluded because the free charge 
density would not be positive definite, as proved in (A), whereas it is 
positive definite for our equation. 

It remains to prove, therefore, that our equation does not decompose 
into the direct sum of the canonical form of the Dirac—Fierz—Pauli equation 
for a particle of spin 3/2 plus the Dirac equation. Suppose such a de- 
composition were possible. Then «* must be transformable through a 
transformation of the type (5c) to the form 

k 
|p whee Meas he ee Geis) 
Lo at , 
where y* are the four matrices of the former equation and «” stands for 
the four Dirac matrices. Without changing the form (45) we can bring 
& to the form 


| Sse eee) 


the first two representations corresponding to the y*. Now we can write 
the «’s of our equation (2) in the form 


j 
Us CVE | c, V" 
ay o(DIk 
a —03W* aga”? bee ye Be dasa (47) 
SS Cale me. azaP 


with the representation # in precisely the form (46). Thus the matrix s 
which transforms (47) to (45) must commute with (46). It must therefore 


have the form 


[ sutas : : | 


s= | P Sgol(4) So3 14) |. Neda en ee (8) 
E : Ssoliy,  S3alcay 


. . . + * x 
where the elements s,, are numbers. Certain submatrices in (48) have to 
0 si : iV i sentations. Now 
be zero since #(2, 4) and A(z, 2) are inequiv alent representatio 
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the top right-hand element in s~taMs is 8,;~'(¢3823-+¢s¢33) VV". This must 
vanish since it does so in (45). Hence, multiplying by s,, which is certainly 
not zero since s is non-singular, 
C3893-+Co83g— 0. ° . . . ° e . (49) 
Moreover, the zero in the last column of (47) must remain unchanged by 
this transformation. Hence 
(8~*)y9893-+(8~ *)o3@3833= 0. 

But by definition 

(8~*)o08n3-+(8~ *)o3833=9. 
Since neither a, nor a; are zero, it follows that 

(sn*) 558550 and (Se )oaheg =O: 
But s,, and s;; cannot both be zero since s is non-singular. Nor can 
(s~*)s5 and (s~'),, both be zero. Hence, either : 
S823=0, 833 70 

or 
But this leads to a contradiction with (49) since neither c, nor c, are zero. 
This rules out the case we are considering and proves that the equation 
defined by (17) is indecomposable. 


APPENDIX 


An explicit representation of the matrices u*, vt and z* is provided by 


1 : 2 : | 


1 | 


wie 


1 
1/2 


u= 6 V/ 2v 


<|H 
bo 


lie 


—_ 


/2 2 
1 
: : _—-—_ 1 
ee 4/2 : 
came 1 * V208= ; 
/2 . 
peal me qi ] 
. /2 2 —— 


with Two Mass States and Positive Charge Density AT 


vi : fa. 
. a © 1 
ia ie | yn 
: ieee ; ? A 20 | ee |” 
ii 1 ice 
a 

2 | 

, | 

| 

; | 

2 


: 
V2 
ie 
= ree feel ks 
ae al 
| 26, TE ela ea 
= 


bole 


w= — 
w= — 


ve 3 
ee eg)? Oe) -C). 


REFERENCES 


Buasna, H. J., 1945, Rev. Mod. Phys., 17, 200-216 ; 1949, [bid., 24, 451-462 ; 
1951 a, Proc. Ind. Acad. Sci. A (in course of publication) ; 1951 b, Report 


of an International Conference on Elementary Particles (Bombay, 1950). 
Dirac, P. A. M., 1936, Proc. Roy. Soc. A, 155, 447-459. 
Frerz, M., 1939, Helv. phys. Acta, 12, 3-37. 


Pats, A., and UHLENBECEK, E., 1950, Phys. Rev., 79, 145-165. 


[. 48 ] 


III. Second Sound and the Thermo-Mechanical Effect at very low 
Temperatures 


By J. C. Warp and J. WILKS 
Clarendon Laboratory, Oxford* 


[Received September 5, 1951] 


SUMMARY 
At very low temperatures the phenomena of second sound and the 
thermo-mechanical effect in liquid helium may be described by using a— 
phonon model independently of two fluid concepts. 


$1. [INTRODUCTION 


Srconp sound in liquid helium was first predicted by Tisza and Landau 
on the basis of the two fluid theory. While this has been successful in 
correlating the velocity of second sound and the viscosity measured by 
oscillating disc experiments, it is still only a phenomenological theory 
with no satisfactory theoretical basis. 

Recent work by Kramers (1951) has confirmed what had before been 
only suspected, namely, that at sufficiently low temperatures (<0-5°) 
liquid helium has a specific heat which varies as the cube of the temperature 
and is about the magnitude to be expected in a Debye type continuum 
in which only longitudinal waves are permitted. The present authors 
(1951) have already shown that a consideration of a Debye type model 
leads directly without: any recourse to a two fluid theory to a velocity 
of second sound equal to the velocity of first sound divided by 4/3 as 
predicted by Landau. Experiments by Atkins and Osborne (1950) had 
shown that at very low temperatures this value is approximately correct, 
but as the temperature is raised it decreases considerably. This fall is 
no doubt due to deviations from the behaviour of the ideal model, but 
so far there is no satisfactory quantum theory of the liquid state to 
indicate what form these must take. In the meantime it may be profitable 
to consider a little more closely the phenomena of second sound in the 
very low temperature region. The object of this paper is to present the 
ideas underlying our previous note and also briefly to remark on the 
thermo-mechanical effect in liquid helium. 


§2. THE VeLociry oF SECOND SounD 


Consider a Debye type model in which the elementary excitations are 
quanta of longitudinal sound waves or phonons. Associated with a 
phonon there is an energy # and a momentum p related by E=pc where 


* Communicated by the Authors. 
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¢ is the group velocity (or the velocity of ordinary sound). The sound 
waves are not independent as otherwise thermal equilibrium could never 
be established ; in all cases there must be interactions between the 
waves. These interactions give rise to the possibility of a periodic variation 
in the phonon distribution function being propagated through the system 
with a definite wave velocity. For this to occur we must assume that the 
interactions conserve energy and momentum since it is only under these 
conditions that an undamped thermal wave can be propagated. Collisions 
in which the momentum is not conserved lead to a thermal resistance — 
and to the heat diffusing in the usual way. 

We previously derived the velocity of second sound by postulating a 
density wave in a phonon gas, and by making use of results already 
available for a. photon gas. However, it is possible to show rather more 
generally that second sound should exist with the correct velocity. We 
first write down the Boltzmann equation which describes the process of 
transference of energy and momentum in the liquid. If f(a‘, p’) is the 
distribution function for the phonons, which are treated as localized 
phonons as in the usual approximations of the theory of thermal 
conductivity of solids, we have : 

of . epi a 

qth F=4.f  (w=lr'p. 
All the phonons are assumed to travel with a constant velocity c and 
A,f is the rate of change of f due to collisions. Assuming conservation of 
energy and momentum in collisions, we then find at once the equations 


OE aPi OPi oT 
oy Es ae 
where ie 
B= | pf Pi= | pifé®p, Pik [EP ra. 


In writing these integrals we have made the assumption that many 
interactions (or phonon collisions) occur in a distance equal to the wave- 
length of any disturbance that may be propagated. This is quite essential, 
for to take an analogy from gas kinetics, it is impossible to propagate a 
density wave in a gas at so low a pressure that it is in the Knudsen region. 
Finally, on remarking that 7'’’=}5'*H for a distribution function zs 
disturbed only slightly from spherical symmetry, these equations give 
directly : eyed rT 

OF 3 dx? 
which is the required result. This proof applies equally to the velocity 
of sound in a relativistic gas, and perhaps clarifies our previous argument. 
We should also mention that we have heard recently from Professor H. A. 
Kramers that he too has devised a proof which though differing from this 
in details is essentially similar. 
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§3. THe THERMO-MECHANICAL EFFECT 


Quite apart from second sound there seems to be another consequence 
of this phonon type of representation. Let us consider a liquid with a 
Debye type excitation spectrum and suppose that it has a vanishingly 
small viscosity. The phonons can be regarded as exerting an osmotic 
pressure in the liquid; this should be observable if one could find a 
suitable semi-permeable membrane, which would be permeable to the 
liquid but not to the thermal waves in the liquid. This requirement 
expressed in the language of two fluid theory means that the membrane 
is to be permeable only to helium at absolute zero, that is to the super- 
fluid. It seems that an effective membrane is a slot or capillary whose 
smallest dimension is of the order of the wavelengths of the phonons. 
The most striking example of the way in which heat may be filtered off 
is seen when liquid helium is forced through very fine channels ; heat 
is developed at the entrance to the channels and the liquid emerges from 
them at a lower temperature. 

Consider now two reservoirs A and B containing liquid helium at 
temperatures 7' and 7'+- 47’ and connected by a narrow capillary. There 
will be a difference JP in the osmotic pressure (P) of the phonons in 
A and B, and since the capillary acts as a semi-permeable membrane a 
difference in the levels of liquid in A and B will be set up equal to AP. 
To calculate 4P we express the pressure P in terms of the internal 
energy H per unit volume and the specific heat per unit volume (C) as 
aT? so that 


4P=44E=1CAT=1aT? AT. 
This is in agreement with H. London’s formula 
AP=SieAT} 


where S is the entropy per unit volume. (Of course any satisfactory 
kinetic description of the thermo-mechanical effort must lead to London’s 
formula as this is known to be a general thermo-dynamic result (De Groot 
1950).) 
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IV. Space Charge Smoothing of Microwave Shot Noise in Electron 
Beams 
By F. N. H. Rozsinson 
Clarendon Laboratory, Oxford* 
‘ [Received September 14, 1951] 


SUMMARY 


A theoretical analysis is given which takes account of both space charge 
interaction between electrons and the multi-valued nature of the flow due 
to the Maxwellian distribution of initial velocities. The theory is in 
agreement with the experimental results of Cutler and Quate (1950) and 
makes possible a coherent account of shot noise at all frequencies. 


Recent work on the travelling wave tube as a signal amplifier (Kompfner 
1947, Pierce 1947, 1950, Robinson and Kompfner 1951) makes it desirable 
to have a theoretical account of shot noise in electron beams at frequencies 
in excess of 2000 Mc/s. In order to be rigorous such a theory must take 
account not only of space charge interaction of electrons throughout their 
flow but also of the multi-velocitied nature of the flow, due to its origin at 
a thermionic cathode. 

We are concerned to calculate the fluctuating component at frequency 
w/27 of the convection current,in a cylindrical electron beam from a gun 
similar to that shown in fig. 1. The shape of the electrodes is arranged to 
cause electron flow in trajectories normal to the cathode. Pierce (1940) 
has shown that the steady state conditions, then obtaining within the: 
beam, are similar to those within a planar diode of infinite extent having 
the same current density and anode voltage. In analysing the gun 
system it is therefore, in the first approximation, permissible to treat it as 
a diode whose anode is permeable to electrons. 
~ MacDonald (1949, 1950) has considered the growth of shot noise in an 
initially smooth electron beam as it passes through regions where the 
potential is constant or varies according to certain laws. He considers 
only the noise generated by the random variation of the transit times of 
electrons due to the velocity distribution and neglects the effect of the 
fluctuating space charge interaction, although this is included in a more 
recent paper (MacDonald 1951) on the equilibrium noise in an electron 
beam. When his theory is applied to an electron gun system used in 
travelling wave tubes it predicts values of noise considerably in excess of 
those found experimentally by Robinson and Kompfner (1951). Clearly 
space charge interaction must play a considerable part in determining the 


noise content of such a beam. 
* Communicated by Dr. K. MacDonald. 
E2 
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Pierce (1950) applies Llewellyn’s (1941) analysis of the diode to the gun 
system and following Rack (1938) assumes that the effect of the Max- 
wellian distribution of initial velocities can be represented by an equivalent 
fluctuation in time of the average velocity of all electrons at the potential 
minimum. This procedure takes account of space charge interaction but 
depends for its validity on the assumption that the electron velocity is 
effectively a single valued function of position and time at points beyond 
the potential minimum. Rack (1938) showed that this is a legitimate 
assumption if the spread in transit time from the potential minimum to 
the anode is a negligible fraction of a cycle at the frequency concerned. 
This has generally been held to be the case, but MacDonald (1950) shows 
that in fact it is not so. Considering a diode in which the potential varies 
as the 4/3rd power of the distance, he finds an expression for the transit 
time 7(H#) of an electron with initial kinetic energy # at the minimum, 


1B)=x(0){1-085(=) |, Aes pe 


Fig. 1 


CATHODE 


COLLECTOR 


FOCUSSING ANODE RESONATOR 
ELECTRODE 


Experimental arrangement for microwave noise study. 


where 7(0) is the transit time of an electron with zero initial energy, and ¢ 
is the anode potential. Typical figures for 7(0) and ¢ for an electron gun 
used at 3000 Mc/s might be 4 cycles and 1000 volts. If we put H=kT, the 
mean energy of electrons originating at a cathode at 1000° x, we find 


7(0)—7(H)=4 cycle= 27/3 radians. 


This will be the order of the spread in transit angle due to thermal velocities 
and is clearly not negligible. 

We are therefore faced with the necessity of finding a theory which, 
while including the effect of space charge interaction, does not neglect the 
multi-valued nature of the electron flow. 

The device of replacing the velocity distribution by an equivalent 
fluctuation will be valid in any region provided that the spread in transit 
angle from one boundary of the region to the other is small. If, therefore, 
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proceeding from the cathode towards the anode we can find a plane such 
that in the region bounded by this plane and the anode this condition is 
fulfilled then we may apply a single valued analysis of this type to this 
region. 

That we can in fact find such a plane follows from eqn. (1). If we 
consider four electrons with initial energies 0, kT’, 2kT, 3kT, their transit 
times (in the numerical example quoted above) are T , 0°91579, 0-89875, 
0°88875. The difference in transit angle between the first two electrons is 
0-677, between the second two is 0-137, and between the last two 0-087. 
While 0-677 is certainly not negligible 0-137 and 0-087 may well be, since 
the effect of a spread dw in transit angle will manifest itself as a difference 
in the phases of the various components of the fluctuation due to individual 
electrons and will be of the order of 1—cos Mw7. This suggests that if we 
choose a plane whose potential is «kT'/e (where «~1) so that the slowest 
electrons have kinetic energy «kT at this plane; then the spread in transit 
angle from this plane, which we shall refer to as the ‘«’ plane, will be 
small enough for it to be permissible to use the single-valued analysis in 
the region beyond it. Direct calculation of the mean reduction in transit 
time from the ‘«’ plane due to the thermal velocities and of the root 
mean square deviation of the transit time shows that this is indeed so. 
The calculation is laborious and it will suffice to note that the situation is in 
fact slightly more favourable to the application of the single-valued theory 
than is indicated by the estimate obtained above. 

Consideration of the results obtained by MacDonald (1949) also makes 
it apparent that the single-valued theory will apply to the region of free 
drift which the electrons enter on passing the anode of the gun. _ The 
initial velocities will not cause any appreciable spread in transit angle in 
this region until some hundreds of cycles have elapsed. 

We, therefore, propose to calculate the shot noise content of the beam 
emerging from the gun on the assumption that it is legitimate to apply a 
single-valued analysis to the region beyond the ‘«’ plane. We shall also 

‘assume that at the ‘«’ plane the current and velocity fluctuations are 
those appropriate to a fully random beam emitted from a cathode at 
temperature 7’. This appears to be a reasonable assumption, since there 
is no @ priori reason for supposing that the potential minimum exerts any 
smoothing action on the components of noise at these frequencies. We 
shall find that the result gives agreement with experiment and that it is 
relatively insensitive to the particular choice of a a’. ‘This latter fact 
gives us some confidence that the device of dividing the analysis at the 
‘%?’ plane is sound. In fig. 2 we show the form of the potential distri- 
bution within the diode and the position of the ‘« ’ plane. ; 

We write the mean square value of the current fluctuation at the 
‘a’ plane in the form 

CROC, Bae ey oe eee C2) 


where J, is the d.c. beam current and B is the bandwidth over which the 
noise will finally be measured. 
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We make use of the result derived in the Appendix to calculate v,? 
the mean square value of the equivalent velocity fluctuation at the 
‘»’ plane. If the mean number of electrons crossing any plane in time T 
is NT, where N—J,/e, then the mean square derivation of their average 
velocity observed in time T from the long term average is 


~~ 


(Sot)}2—(NT)—“1[o?— (3)3], 


where v? and v are the mean square and mean of the velocity averaged over 
the velocity distribution. ; 
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Potential distribution and position of the ‘«’ plane in a space charge 
limited diode. 


By a theorem due to MacDonald (1949) we may then obtain the power 
‘spectrum w(f') of the velocity fluctuations. 


; OO a 
w(f)=Anf |” Sar sin 2nfT aT, i getie 7) 
where 
P= Taso" » cy ae 
The mean square equivalent velocity fluctuation is then 
v2=u(f). B. fete ae ene 
At the ‘«’ plane 
~ (2 (2B Oak dE 
A ae > PY 
; Sil (——) exp (—E/kT) 77, (7) 
~  {° (2H 2akT\ 12 dE 
eae a exp (—H/ET) Fr (8) 
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Evaluating these expressions and substituti i 
tuting the results in equations 
(3), (4), (5) and (6) there results | : 


— e 4kTB 

Eager i F(@), é 3 . ° ° = -, (9) 
where ; 
oe f(%)=14-a—[ot 2+ be%7l/2(1 —erf 1/2) 2 oY eon unt LO) 
n 


x 
erf 2-1 | e dt. 
0 


We tabulate f(«) for a range of values of « in table 1. 
for «>20 f(a)~(4a)-3, l 
Fora=0 f(a)=1—n/4, | 


and eqn. (9) reduces to the expression found by Rack (1938) by a less 
direct and more intuitive method. 


(10 a) 


Table 1 
om f(x) a F(x) a f(a) aX f(a) 
0 0-215 0:8 0-109 2°5 0-060 10 0-020 
0-1 0-182 1:0 0-099 3:0 0-053 13 0-016 
0-2 0-162 1:3 0:087 4-0 0-043 16 0-014 
0-4. 0-140 1:6 0-079 5:0 ~=0:036 20 0-011 
0-6 0-123 2-0 0-069 7-0 0-028 


Expressions have been given by Llewellyn (1941) based on a single- 
valued analysis for the fluctuating components of potential V, convection 
current c and velocity v at any plane in the diode in terms of their values 
at an earlier plane (which we shall identify with the ‘«’ plane). V,, ¢, 
and v,. 

V—V,=Ay +A +A 196, 
v= AT, +Aoo¥,+AosCy red ae LL) 
C=A gl +A5qv,+Azsey | 
where J, is the fluctuating component of total (sum of convection and 
displacement) current and the coefficients A are functions of frequency, 
d.c. current density and the coordinates of the two planes. 

At lower frequencies (<100 Mc/s) shot noise is observed in a circuit 
connected between the anode and cathode of the diode and the measured 
quantity is J, or V. In the experiments with which we are concerned the 
noise is observed by measuring the voltage induced in a resonator or helix 
when the beam passes through it after leaving the gun. The observed 
quantity is the fluctuation c in the convection current. Fluctuations in 
anode potential (V) and total current (Z,) are not therefore of interest. 
Moreover it may be shown that if the total transit time across the diode is 
greater than two cycles the error involved in neglecting terms in V and J, 
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is less than 3°. We can therefore simplify equations (11) considerably 
and also do net have to consider the effect of any external impedance 
connected between the anode and cathode of the gun. 

Eqns. (11) are thus reduced to 


=A 90, +A 3Cy : 
C=As.v,+ A g3ly 

It is convenient to express the coefficients A in terms of J, the steady 

diode current, 7' the cathode temperature, w the angular frequency, 7) the 

transit time of an electron from the potential minimum to the anode in the 


absence of any initial velocities and the normalized potentials « and 
n=ed/(kT’) of the ‘«’ plane and the anode. We then have 


(12) 


A,.=Fy(«, 9), A,,=—K F(a, 7), | (13) 
Az.=KF,(«, 7), Az,=F,(«, 7), 
where : 
K=jaq l,(2el i) = eee ee 
Fy =1—2[y](1+y)-!2{[y]— [«]}, - ¢ eo 2) SORE 
F.=(1+9)-"4{14 [olin inl—[e]}, - - -. . . (158) 
F=(1-+n)94(nl [el en ee ee eee 
Fy=1—(1+y)7{-fe]?, ee. (8) 
and 
Real et Css PE 
(16) 
SG Ue Ube 


Now in using eqns. (12) to evaluate the fluctuating currents and 
velocities at the anode we must bear in mind that c, and v, are uncorre- 
lated. Their contributions must therefore be caleulattd separately and 
then summed as squares. 

Denoting by c, and v, those components of ¢ and v derived from ¢c, and 
by c, and v, those derived from v, we have 


CPC, * Gay eh as be ee eee (17) 


with a similar equation for v2. 
From eqns. (2), (9), (12) and (13), 


C= F",(2ef, BR); a Pen mme oh G b 8, 
mee cn aha SUR te nae aice eBoy r 
Co—KF3(e/m)*?(4kT Bf(a)/Ip)1?, . 2 . . (18 e) 
Vg—=F (e/m) 2(4kT Bf(a)/Io)¥?%. 2. in, (18a) 


Expressing c? in the form 


Cias[ 7 261, Be eee teeta a at} 
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and using eqn. (14) there results, after some re-arrangement, an 
expression for the space charge smoothing factor I? at the anode of the 


piode. PSF 2+ (orl flwF 2, Mt AE E00) 
where F,, F', and f(«) are defined in eqns. (15) and (10). 

We may also calculate the value Ip? of the space charge factor which will 
be found after the beam has drifted for a time Tp in a field free space 


beyond the anode of the gun. The current fluctuation Cp is then given in 
terms of the current and velocity fluctuations c and v at the anode by 


Cp=C COS WoT p—J(w/w9)Lq(v/w) sin woTp, Se ve OCS) 
where w is the steady beam velocity and wy is the plasma frequency in the 
beam (Ramo 1939). 

In applying equation (21) we use the relations 
Ae (2b Tn PEL ten NEN ts Sse a 22) 
WoTp9— a/ 2. . . . . . . . . . (23) 
Again evaluating contributions due to c, and v, separately and taking 
the sum of squares there results : 
Dp {fF P+ oF (1+n) >} sin? (wtp +X;1) 
+(w7)*f(a){F5°+ $F 7(1+7)-+} sin? (wotp+x2), (24) 
oo tan x,;=1/2(1+7)1?F,/F | 
tan y¥2=— V/2(1-+n)?P3/F, 
The dependence of [’,? on the arbitrary parameter « is not in eqn. (24) 
very perspicuous ; we may note, however, that for «=0 and 7 large, 
_ eqn. (24) reduces to the equation derived by Pierce (1940) on the basis 
of an entirely single-valued theory, namely, 
DypP=$hkT|(eb)(w79)"(1—m/4) sin? (weTp +x): P20} 
In order to elucidate further the dependence on « we shall use as a 
numerical example the experimental data of Cutler and Quate (1950). 
These authors explored the variation with 7p of J’, using a gun which had 
a cathode temperature 7'=2000° kK, transit angle w7y=32 radians and 
“anode potential d=900 volts. Writing eqn. (24) in the form 
ty 4 sin? (woTp txXi)+8 sin? (wetptxe), - - + (24a) 
we tabulate A, B, yx, and x, for a range of values of « for this gun in 
table 2. In the last row we have tabulated the amplitude of I)? 
neglecting the small difference between y, and xy» which alters |Ip?| by 
less than 1%, and also gives rise to a small constant term of the order of 
10-4. It is immediately apparent, as we predicted, that the result 
obtained is not critically dependent on the particular value of « between 0 


AS é 
ae experimental values found by Cutler and Quate can be written 


in the form 


(25) 


['p2=0-055 sin? (wtp +66°)+0-01. . . . . (27) 
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The small constant term is attributed by these authors to partition noise. 
When we bear in mind that microwave noise measurements are seldom 
more accurate than 5% and that the type of experiment carried out by 
Cutler and Quate would tend to overestimate the magnitude of Ip? the 
agreement with the results in table 2 appears to be satisfactory. 


Table 2 
Ce eee 
Chey + 4 1 2 4 n 
A= 0:0003 ~=0-010 0:0172 0-0283 0-0451 0-0737 i! 
B= 0-061 0-0405 0:0316  0-0231 0-0128 0-0070 0 
X1= 553° 68-5° 69-7° 70° 70-2° 11-2. 90° 
Y,= 0075. 56°7~ 56-9° 57-9° 58-7° 59-7° 0° 
[52 | 250-061 0-051 0-049 0-051 0-058 0-081 1 


However, in carrying out the analysis we have neglected the effect of the 
finite size of the beam. It is possible (Ramo 1939) to introduce a correction 
for this in the field-free drift region and in interpreting their results in the 
light of Pierce’s theory Cutler and Quate did this. In this region the 
finite size of the beam results in a reduction of the effective plasma 
frequency wy). Inspection of eqn. (21) shows that this results in an 
increase in the contribution to [’,)? from velocity fluctuations at the anode 
and a decrease in y. Making this correction Cutler and Quate found good 
agreement with the experimental value of the periodicity wy; but on the 
basis of Pierce’s equation (26) they obtained a thegretical value 

P7?=0-089:sin? (@otp +42); - - se 
which is outside the limits of experimental error. The present theory with 
the most favourable value of « is capable of yielding a value | I’)? |=0-065 
when these corrections are made but the discrepancy in the value of y still 
remains. : 

Cutler and Quate suggest that this is due to the neglect in the analysis 
of the finite beam size within the gun. Although an exact treatment of 
this region has not yet been attempted, a qualitative estimate does indeed 
suggest that the effect of finite beam size will there tend to compensate 
the effect on | I"p? |, x; and x, of the corrections applied in the region of free 
drift. In the absence of a complete theory of this effect the agreement 
obtained is as good as can be expected. 

At very low frequencies eqn. (20) predicts that [2 > kT (ed). The 
discrepancy between this result and the value [?=0-644k7'/ed obtained by 
Thomson, North and Harris (1940) by a quasi-static analysis, arises from 
the neglect of terms in J, in eqn. (11) which is only “permissible for 
frequencies such that wt, > 4. . ; 

At extremely high frequencies Pierce’s (1950) single-valued analysis 
leads to the conclusion that I? can increase indefinitely as w increases and 
reach values in excess of unity. In the present theory this difficulty does 
not occur, since as the frequency increases the region, within which a 
single-valued analysis applies, contracts. That is to say, that in order for 
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ee mean spread in transit angle 4w7 to be less than any pre-assigned value 
oe and larger values of ‘x’ have to be taken until finally the ‘« ’ 
plane coincides with the anode when eqn. (20) predicts that [=1. 


If the assigned value of the issi 
. permissible spread in t t angle i 
from eqn. (1) « is determined by f eh aaah” 


| 0°85.A-twr9 {kT'/(ed) }/4{(a+1)4—at/4t—= 1, 2. (28) 
ae have only qualitative considerations to guide us in the choice of J. 
A <n/2 space charge interaction will already play a considerable part 
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Space charge smoothing factor as a function of frequency. 
A. Pierce’s theory. B,C, D. Present theory with d=7/2, 7/4 and 7/8. 


in determining the noise content of the beam, but the single-valued 
analysis will give only a very rough account of the behaviour in the region 
beyond the ‘«’ plane. If on the other hand 4 <7z/8, the single-valued 
analysis will be a very good approximation to the behaviour beyond the 
‘ ’ plane but we shall have neglected any ordering of the beam which may 
have taken place before the ‘«’ plane is reached. It, therefore, appears 
reasonable to expect that a correct description will be found with 
1/8<A<n/2. In fig. 3 we plot | Ip? | as a function of frequency between 
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1000 and 100 000 Mc/s for 4d=7/8, 7/4 and 7/2. The values of 7, and 7 
are those appropriate to the gun used by Cutler and Quate discussed above. 
For the purpose of comparison we also give the value of |I°p?| from 
Pierce’s equation (26) and Cutler and Quate’s experimental result. 
Throughout this region the difference between x, and xz is less than 12° 
and is neglected. 

It is satisfying that it is possible to show that the existence of electron 
beams showing a considerable degree of space charge smoothing at micro- 
wave frequencies does not imply that the potential minimum exerts any 
smoothing action on such high frequency components of shot noise. It is 
indeed difficult to see how this could occur, at frequencies so high that the 
transit of electrons from cathode to minimum occupies several cycles. We 
have in fact assumed that thestream of electrons emerging from the potential 
minimum is entirely random. In a forthcoming paper Kompfner (1951) 
will show how the interpretation of micro-wave partition noise given by 
Robinson and Kompfner (1951) and depending on the compensating 
action of the potential minimum, can be adapted to the view of the space: 
charge smoothing process given above. 

The theory we have obtained differs from earlier theories (Pierce 1950): 
chiefly in the interpretation it affords and in eliminating the possibility 
of I? exceeding unity at sufficiently high frequencies ; but at least one 
practical consequence is worthy of note. 

The fact that the current and velocity fluctuations consist of two un- ~ 
correlated components whose minima do not coincide (x, y,) limits the 
reduction in noise possible by the method proposed by Field, Tien and 
Watkins (1951). In this method the beam is suddenly accelerated and. 
then at a later point decelerated, the points of acceleration and deceler- 
ation coinciding with the minima of the current and velocity fluctuations: 
respectively. In the simplest case when the amplitude of the two com- 
ponents associated with y, and y, are equal in magnitude the greatest. 
reduction in noise possible will be approximately sin 2(y,— y,). 

In conclusion it may be instructive to give a brief account of the course: 
of fluctuations in a beam of electrons from the time it leaves the cathode. 
Between the cathode and the potential minimum and beyond the potential 
minimum to the region of the ‘«’ plane the behaviour is dominated by 
the Maxwellian velocity distribution and the beam remains effectively 
random. From the ‘«’ plane to the anode the flow becomes effectively , 
single-valued and Llewellyn’s (1941) analysis correctly describes the 
influence of space charge interaction which is here the dominant feature of 
the flow. Beyond the anode the flow is also at first effectively single- 
valued, a treatment such as that of Ramo (1939) is appropriate, and in 
this region the amplitude of the fluctuations varies periodically with 
distance. : Eventually, however, since the flow is not strictly single-valued 
the velocity distribution will again dominate and the periodic variations 
will decay until finally the current and velocity fluctuations reach the 
equilibrium value predicted by MacDonald (1951) on thermodynamic 
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grounds. This final transition to the equilibrium state has not been 
analysed but it appears plausible to assume that it will occur after a 
transit angle such that the spread in transit angle due to the Maxwellian 
distribution becomes comparable with. According to MacDonald (1949) 


this will require a total drift time of some hundreds of cycles in even the 
most favourable case. 


APPEND sx 


We consider electrons crossing a fixed plane at a mean rate NV per second, 


with an average distribution i in velocity n(v) so that the long term average 
velocity is 


(A 1) 
If we observe the electrons over a period T during which time the actual 
distribution is n;+v, say, then the average velocity observed in this 
interval will be 

~~ (Miers); 
iS Tint bate A2 
tens rape tame ps i? ue en re 

i 


Now if electrons in each velocity class v,; are emitted independently and 
at random, 
Clenda Ene ane poate eae he (AND) 
where the average is over many periods T 
Provided T is long enough and the subdivision into velocity classes is 
coarse enough so that T,; is large, v; will be much small than 7; and we can 
expand (A 2) as 


and so 


et part rae vi 8yi} ° . . . : (A 4) 


We now average (Sut)? over many periods T and obtain 


—— Ta a eee 
= NEL p90, W2Lv 99+ OPEZ 
Nod % 4 j 
(A 5) 


Since the fluctuations in each velocity class are independent this reduces ° 
to’ 


(Sohp a2 [ay tod— Iv Po/ HOPE. yA) 
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We now make use of eqn. ) to obtain 
(Sot a= (Vole ie iden etic i)¢2n,} 
—/(NT)-1fo2— (6)"},, ee 


which is the result used in the text. 
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V. The Coulomb Scattering of Charged Particles in Nuclear Emulsions 
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SUMMARY 


Measurements have been made on the multiple coulomb scattering 
of electrons and protons in normal and diluted Ilford G5 emulsions. 
Empirical constants are given relating the mean angle of scattering to the 

product momentum x velocity for particles of unit charge. 


§1. LyrropuctTion 


THE nature and energy of a charged particle traversing a nuclear emulsion 
may be established by making observations on any two of the following 
three characteristics of its track :— 

(a) The residual range. 

(6) The multiple coulomb scattering. 

(c) The ionization (grain density or delta rays). 

The grain density produced by a charged particle is partially dependent 
upon the method of processing of the emulsion, and can vary from one 
batch of plates to another despite apparently rigid conditions of develop- 
ment. For an emulsion of a given composition the residual range and the 
magnitude of the coulomb scattering of particles with the same mass, 
velocity and charge are always the same, however, within the limits 
imposed by statistical deviations. These two properties may, therefore, 
be regarded as the ‘ constants’ of a nuclear emulsion, and their values 
will depend upon its atomic composition. 

The extremely sensitive nuclear emulsions employed in recent years have 
contained a very high silver halide to gelatine ratio. Their sensitivity was 
mainly dependent upon this feature. Recently, however, Dodds and 
‘Waller (1951) have introduced new types of emulsions, which, whilst 
capable of recording the passage of particles at minimum ionization, 
possess a much smaller silver halide content. These are known as diluted 
Ilford G5 emulsions. 

In this paper values will be given for the ‘scattering constants’ of 
normal, <2 and x4 G5 emulsions (that is, emulsions with the normal 
gelatine to silver halide ratio and with twice and four times the normal 
value of this ratio). A later paper from this laboratory will give range— 
energy relationships for the diluted emulsions. 


* Communicated by Professor P. I. Dee, F.R.S. 
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§2. Mernop 


Measurements have been carried out using protons with kinetic energies 
of 58, 77 and 146 Mev, and electrons with energies of 9-5, 13-5 and 
18-5 Mev. Thus the coulomb scattering has been examined for particles 
of unit charge over a wide range of velocities and masses. 

The proton tracks were measured in plates exposed to the external 
uncollimated beam of the Harwell cyclotron, and the two lower energies 
obtained by interposing aluminium absorbers in the path of the beam. 
The electron tracks were obtained by pair production in the x-ray beam 
of the smaller Glasgow University synchrotron. Electrons with the 
required momenta were selected by means of a magnetic channel. 

The proton scattering measurements were carried out over track 
segments 8000, in length. Segments of 1500 were used for the 
electrons. The scattering measurements were performed on a Cooke 
M.4000 microscope with a x 45 objective and x 15 eyepieces, using the 
coordinate technique developed by Fowler (1950). . 

Grain-density determinations were also carried out on the proton tracks 
in both normal and diluted emulsions in order to examine the variation of 
grain-density with change of velocity. For these measurements a 95 
objective was used. 


$3. RESULTS 
(i) Presentation 
A review of the theories of multiple coulomb scattering of charged 
particles has recently been given by the Bristol group (Gottstein et al. 1951. 
We are grateful to the authors for making their results available to us 
before publication). Using the notation adopted by these authors, the 
mean angle of scattering may be given as 


2e?(N Z?)1/241/22 
eee 
ppe 
where (® ).norq 18 the mean angle of scattering. measured by the coordinate 
method. WN, Z, ¢ are the number of scattering atoms per cm, their 
atomic number and the thickness of the medium in microns (in our case the 
cell-size). ze, p, Be are the charge, momentum, and velocity of the 
scattered particle. L is a slowly varying function of NV, Z, ¢ and B. 
Values of L as a function of ¢ and B have been given in graphical form 
by Goldschmidt-Clermont (1950), using the Moliére theory (1948). 
Eqn. (1) may be reformulated as 


1 
{® aneet x bBe x (+ 


(® chord = L, : : . . ° (1) 


~ 


1/2 
) =2¢(NZ?)12L(100)'2—K, . . (2) 


where K is termed the ‘scattering constant’ of the medium. It is 


usually expressed in units of degrees Mev (100 .)-!!2, and we shall employ 
this convention. ’ 
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We shall also use a parameter K,—the scattering constant determined 
from a restricted selection of the experimental values of O,,.4- This 
parameter is of particular importance in the determination of the energies 
of very fast particles traversing the emulsion. In such determinations 
comparatively few angular measurements can be made on a single track, 
and the presence of one large angle single scattering can appreciably alter 
the mean value of the scattering angle. In the present work we have 
employed the convention that angles greater than four times the median 
value are not included in the determination of K,,. 


In our presentation of the results corrections have been made for the 
following effects :— 


(a) Spurious scattering introduced by observational and microscope 
errors (Levi Setti 1951). 
(6) The inclination of the track to the coordinate axes employed for 


measurement. (This correction was necessary only for the electron 
tracks.) 


Table 1 
Kinetic energy Ki 
| (after correction 
for loss) t=25 p 50 100 200 400 
8:7 Mev | 28:-2+0-6 | 27-1+0-8 | 26-5+1-2 | 25-8+1-7 | 25-3+2-6 
13-4 — 27-4+0°8 | 27-341-2 | 28-441-7 | 26-5-+2°5 
Dze2 — 28:3-40°8 | 26:5-+1-0 25-7-41-5 | 24-:9-12-2 
Vie 
8:7 21-4-+40°5 | 21:9+0-7 | 23-5-40-9 | 24:5+1-7 | 25-3+2-6 
13-4 — 23-°3+0-7 | 24-5+1-1 | 23-2+1-6 | 26-5-+2:5 
17-2 — 23:-2+0-7 | 24-9+1-0 | 22-4+1-6 | 24-9+2-2 


(c) Energy losses by the particles in passing through the emulsion 
(Rossi and Greisen 1941, Smith 1947, Montgomery 1949, Bradner et al. 
1950). 
' (d) Distortion of the electron tracks in the x 4 emulsion. ° A frequency 

distribution of the algebraic values of O4,..q was plotted, and it was found 
that the median was appreciably displaced from the zero position (Menon 
et al. 1951). The displacement was assumed to represent the distortion, 
and its value was added algebraically to the individual measurements of 
Oonora: Tests on the proton tracks in the x 4 emulsion, and on electron 
and proton tracks in normal and x 2 emulsions indicated that distortion 
was negligible. 
(ii) Baperimental Results for N: ormal G5 Emulsions | 3 

The experimental values of K and Kk, for electrons are shown for various - 
energies and cell-sizes in table 1; the errors quoted are the statistical 
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‘probable errors’. It can be seen that, within the limits imposed by the 
statistical error, K and K, are independent of energy, in accordance with 
theoretical predictions. The values of (® ) hora X pBe as & function of pc 
are given in fig. 1. 


Fig. | 
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Comparison of experimental results for the electrons with Moliére’s theory 
of multiple scattering. 


In fig. 2 the weighted mean values of K and K, for the three energies 
are plotted as a function of cell-size ; the theoretical curves for K and K ; 
derived from Goldschmidt’s paper and the Moliére distribution function, 
are included in this diagram. It can be seen that the experimental and 
theoretical values of K do not agree for small cell-sizes. This lack of 
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agreement can possibly be explained by a small difference in the theoretical 
and experimental frequency distributions of the values of nora for large 
scattering angles. Better agreement is obtained between the experi- 
mental and theoretical values of K,. : 
The value of K at a cell-size of 100 p is in good agreement with that 


obtained by other workers (Gottstein ef al. 1951, Corson 1951, Voyvodic 
and Pickup 1951). 


Table 2 
Kinetic energy K 
(after correction 
for loss) t=200 p 400 
140 Mev — 30-0 0-8 
66 33°0-40-6 | 31-9-L0-8 
47* 27-1+0-7 | 26-2+1-0 
IG, 
140 Mev — 28-6 0-8 
66 28-3-40-6 | 28-6-+0°8 
47* 24:2+0-7 | 24:-2+0-9 


* This value is based on residual range measurements 
of protons in the emulsion (Bradner et al. 1950). 


Table 3 


Kinetic 


Type energy K | 
of (after 
particle | correction 
for loss) t=50 pu 100 200 400 800 
Electrons | 18-9Mmev | 22-7--0-7 | 21-0+0-9 | 21:71°3 | 22-5 +2-4 eae 
Protons | 141 = ses = 23-4+0-8 | 22-2+1-2 | 
69 — — 24-7-L0-6 | 22-6 -+0-7 | 20-7-+41-0 | 
49 — — 27-9-+40-6 | 25-7 0-9 | 25-0-+1-1 


K 


c 


Electrons | 18-9 16-7 +0-5 | 16-8 0-7 | 20-0 1-2 | 22-5+2-4 = 
Protons 141 — — — 21-5-40°8 | 22-0-11-2 
69 a — 21-9 0-5 | 21-3-+40-7 | 19-8-+1-0 


The value of K and K, for protons with mean kinetic energies of 140, 
66 and 47 mev (after correction for ionization losses in the emulsion) are 
given in table 2. There is good agreement with the theoretical values 
put the expected small increase in K with increasing cell-size is too small 


to observe in this case. 
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(iii) Diluted Emulsions 

The experimental values of K and K, for x 2and x 4 diluted emulsions 
are given in tables 3 and 4. It can be seen that there is little or no 
variation of K with cell-size. 

The values of K and K, for the electrons in the x 4 emulsions cannot be 
considered as reliable as those for electrons in the other emulsions. In 
addition to the distortion mentioned previously great difficulty was 
experienced in following the widely spaced grains which marked the passage 
of the electrons in these emulsions. 


Table 4 
Kinetic 
Type energy Kk 
of (after 
particle | correction 
for loss) t=50 p 100 200 400 800 
Electrons | 19-1 Mev — 19-2-+0-9 | 17-6+1:3 | 20-4+2-3 a 
Protons 142 — — — 23-1+0-8 | 21-7+1-1 
69 _: — 19-9-+0-7 | 19-8-40-7 | 18-4+1-0 
50 — -- 22-6+0-6 | 19-9-+40-9 | 18-8-+1-1 
Ke 
Electrons | 19-1 — 13-6 +0-7 | 13-8-+1-1 | 20-4+2:-3 — 
Protons 142 — — — 17-3-40-7 | 19-4+1-1 
: 69 — — 16-1 -+0-4 | 18-6-+0-6 | 17-5-+0-9 
50 — as 18-4-40-5 | 17-9-+40-7 | 18-6-+1-0 
Table 5 
Protons, t=400 Electrons, t=100 p 
Emulsion eS 
| 141 Mev | 68 | 49 +5 Mev 
Experiment | 30-0 +0-8 | 31-9-L0-8 | 26-2 +1-0 26-7 -+0-6 
Normal Theory 28-6 28-9 29-0 25-5 
Experiment | 23-4-0-8 | 22-6 0-7 | 25-7 40-9 21-0-+40-9 rt 
2 Theory 23-4 23:7 23:8 20-7 
Experiment | 23-1 --0-8 | 19-8-£0-8 | 19-9-40-8 19-2 +0-9 
x4 Theory 19-5 19-6 19:7, Piel 


{iv) Summary of Scattering Results 


In table 5 the values of K for electrons and protons in the three 
| emulsions are compared with those calculated from the Moliére theory ; 
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ees ee overall agreement. The cell-sizes were chosen so that the 
Ae as eee, is high and any residual effects due to spurious scatter- 

Empirical formulae relating (® ehorac» the mean angle of scattering 
with cut-off, to the product of momentum and velocity are given in 
table 6. These take into account only the variation of K,with 8. The 
rapid change in the value of K, for electrons using cell-sizes smaller than 
100 » has been neglected since considerably larger cell-sizes are normally 
used for experiments which apply the scattering technique. In presenting 
these formulae we have taken into account the data of Gottstein et al. for 
the scattering of protons with energies of 10-20 Mev and. 340 Mev in 
normal G5 emulsions. | 


Table 6 


Piped NG 
(® denora =k, ( ) — degrees. | 


100 pBe 

Kmulsion 
aoe A pea ee 
orma. KD enocdie [27 5—3 5f(B)] (a) a Be 
way) =|22-5 — Sa ae 
x =[22-5—4/(6)] a) x oR 
x4 =[17-5—8 Ea 
=17 5-3/1 (a55) sae 


(v) Variation of Grain Density with Velocity 

An investigation was carried out on the variation of grain density in the 
proton tracks as a function of the velocity of the proton for the three 
types of emulsion. In addition to measurements on the grain density of 
the proton tracks, values were obtained for g,,;,, the grain density at 
minimum ionization in the plates. The tracks of cosmic ray particles of 
great energy and identifiable fast electrons were used for the latter 
determinations. 

The results are presented in table 7 as the ratio of grain density at a 
velocity Bc to that at minimum ionization ; figures for the latter quantity 
are given in the final column of the table. It should be noted that the 
values of g,,,, do not enable a comparison to be made of the sensitivity 
of the different emulsions, since it was found necessary to employ a 
stronger developer for the diluted emulsions than for the normal ones. 

The statistical errors only are quoted in the above table, in addition 

ystematic errors can arise as mentioned in the introduction to this paper. 
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It can be seen, however, that for a given value of the velocity the grain 
density ratios are reasonably constant for the three types of emulsion. 
In fig. 3 the weighted mean values of the ratios are plotted as a function of 


Table 7 
Emulsion BY Os4 ve Gaui 
Be=0-495¢ | Bc=0-36c | Be=0-3le dhe 
Normal 2-8+0-1 5-2-0-1 6-9+0-2 14:9-40°3 
x2 S1280:1 4-3-40-1 5840-2 | 10-0-40-2 
det 3-140-1 5-2-10-2 6-2-0-2 5-7 -£0-15 


* gnin—grain density at minimum ionization in grains per 50 p. 


Fig. 3 


GRAIN DENSITY/GRAIN DENSITY AT MINIMUM IONISATION. 


3 4 Si rece Gy he) 
B 


The variation of grain density at a velocity Bc to that at minimum ionization 
as a function of B. 


the velocity on a logarithmic scale ; in the velocity range Be=(0-3 + 0-5)e 
the data may be represented by the relation 
Grain density at velocity Bc 


; ~ p18, 
Ymin 
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SUMMARY 


Evidence is given that cadmium iodide crystals grow from a low degree 
of supersaturation by the dislocation mechanism. Sequences of 
photographs illustrate the stages in the growth of the crystal plates. 
The screw dislocations are not simple but are groups of simple 
dislocations. The kinking of step-lines is explained by the interaction 
of screw dislocation groups of different strengths. The properties of large 
screw dislocation groups in the crystal plates are discussed. 


§1. INTRODUCTION 


THE growth of cadmium iodide crystals from aqueous solution has been 
reported earlier in a preliminary note (Forty 1951). It has been found 
that, at low degrees of supersaturation, the crystal plates develop in 
thickness by the dislocation mechanism of spiral growth, proposed by 
F. C. Frank (1949, 1950; see also N. F. Mott 1949). The successive 
stages of growth of the crystals can be studied under the optical 
microscope. This paper gives a more detailed account of some of the 
observations. 


§2. MprHops OF OBSERVING THE GROWING CRYSTALS 


The crystals are grown between a microscope slide and coverslip from 
a drop of solution, just saturated at 60°c¢ and allowed to cool slowly 
to room temperature. The crystal surfaces are observed in reflection, 
using an ordinary optical metallurgical microscope. The step-heights 
of the growth layers are usually large (up to 2000 A high; a study of 
step-heights on cadmium iodide crystals is described in Paper II) and 
the step-lines have a high visibility. Several methods for observing the 
step-lines have been employed. 
(1) Vertical, Bright Field Illumination 

The growing surfaces are observed in reflection ; the refractive index 
of the crystals is high and hence the reflectivity is considerable. By 
* stopping-down ’ the illuminating beam to a narrow pencil, step 
visibility is improved (Griffin 1951). Figs. 5-10, 14-23, 24, 26 
(Plates I, IT and IIT}) are photographs using this method. 
ee ee ee eee 

* Communicated by Professor N. F. Mott, F.R.S. 


+ For plates see end of issue. 
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(2) Bright Field : Double Transmission 


When the steps are high, some contrast may be sacrificed to obtain 
stronger illumination for short exposures of photographs. This is simply 
achieved by placing a mirror behind the crystal. Figs. 11-13, 25. 
(Plates I, IIT) are microphotographs taken in this way. 


(3) Reflection Phase Contrast Illumination 


This mode of illumination is found especially useful for observing 
small steps on dry crystals with fairly simple surface structure. 
Fig. 27 (Plate IV) is a microphotograph of a dry, unsilvered crystal 
taken with Cooke, Troughton and Simms metallographic reflection 
phase-contrast equipment. 


-/§3. OBSERVATIONS AND Some DeEpucTIONS 
(1) Growth Sequences 
After nucleation, the crystals grow in the form of extremely thin, 
hexagonal or triangular plates. The well-developed (0001) faces, growing 
parallel to the cover-slip are normally observed. The initially rapid 
lateral growth of the plates may well take place without the aid of 
dislocations, on the less closely packed prism faces of the hexagonal 
structure; the initial degree of supersaturation of the solution is about 
50%. The rapid growth of a large number of small crystals lowers the 
degree of supersaturation and the lateral growth of the plates soon 
decreases. Then growth in thickness is observed by the appearance of 
several small hills near the edges of the plate. Eventually these are 
near the centre of the crystal as a result of continuing lateral growth. 
As the growth rate decreases, the hills can be resolved into step-line 
patterns centred on systems of screw dislocations emerging on the 
(0001) faces. 

The structure and behaviour of growth hills agree with the theoretical 
predictions (F. C. Frank 1949, 1950; Burton et al. 1951), but it must 
be emphasized that the growth layers in this case are not ‘ monomolecular ’. 
Sequences of photographs can be taken to illustrate the stages in the 
growth of these crystals. 

(a) Growth based on a Single Screw Dislocation Group 

The step-line terminating on the emergence of the dislocation group 
on the (0001) surface winds itself into a spiral. Figs. 5-10 (Plate I) are 
photographs taken at intervals of ten minutes to illustrate the 
development of the spiral step-line about a screw dislocation group as 
the growth proceeds. ) 

(b) Growth based on Two Screw Dislocation Groups of the same hand 

Figs. 11-13 (Plate I) are photographs at 10-minute intervals 
illustrating this type of growth centre. The co-operation of two screw 
dislocation groups of the same hand produces a step-line pattern of a pair 
of spirals, one following the other. The two cases described by Burton. 
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et al. (1951), namely for the distance between dislocations 1>2p, and 
for 1<2zp,> where p, is the,radius of the critical nucleus for growth to 
occur, have been observed. The difference between these two cases is 
the interchange of step-lines away from the centre. If the screw 
dislocation groups are very close together, the resultant step-line pattern 
appears to be a double spiral with a single centre, as in figs. 11-13 
(Plate I). 


(c) Growth based on Two Screw Dislocation Groups of opposite hand 


Figs. 14-18 (Plate IL) (taken at 12-minute intervals) and figs. 19-23 
(Plate IL) (taken at 60-minute intervals) show the stages in the expansion 
of the terrace between two screw dislocation groups of opposite hand, 
dominating the growth of a crystal face. The outward growth of the 


terrace generates a closed loop and a new terrace. This cycle is repeated _ 


indefinitely and the crystal thickens at the rate of one layer thickness 
per cycle. 
(d) More Complex Growth Centres 
Growth centres of many screw dislocation groups have been observed. 
All show that the thickening of the crystal plates is due solely to the 


dislocation mechanism of growth. Fig. 24 (Plate III) is part of a surface — 


of a crystal containing many screw dislocation groups. At the top of 
the photograph there are at least three independently active groups, not 
all of the same strength (or Burgers vector normal to the face). On 
the right there is growth from a pair of screw dislocation groups pf the 
same hand, and in the centre there is an irregular spiral based on a single 
screw dislocation group. 

The turns of this spiral are more closely spaced and were seen to be 
advancing faster than the others (fast enough to blur the image slightly 
in the 75 second exposure). Many other ‘smaller dislocation groups are 
present, producing kinks in the step-lines (see later, (5)). It is interesting 
to note that the central spiral eventually spread over the whole surface 
and figs. 5-10 (Plate I) are a sequence of its growth twenty-four hours 
later; the layers becoming polyhedral as the growth rate decreased. 


(2) The Shape of Growth Layers 


It is always observed that, when the dislocation type of growth 
commences, the step-lines are curvilinear, the layers becoming polyhedral 
as the growth rate decreases with the fall in degree of supersaturation. 
From many observations of growth, the phenomena can be summarized :— 


Initial fast growth from a high degree 


: of supersaturation Final slow growth 
(a) Steps close together ..............-.- - Steps widely spaced. 
(6) Fast rotation of step-line about screw 
dislocation 2).:.,.. <3 ee Slow rotation of step-line. 


(c) Step-lines are curvilinear.............. Layers polyhedral. 
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The shape of the polygonal layers displays the symmetry of the crystal 
surface. This is trigonal for a basal surface of cadmium iodide. In 
figs. 14-18 (Plate II) the polygonal growth steps have the angles of a 
regular hexagon but the sides. are alternately short and long. In 
figs. 19-23 (Plate II) which show the same crystal at a later and slower 
stage of growth the short sides of the hexagons have undergone a rotation 
of 30 degrees so that alternate angles are close to 150 degrees and 
90 degrees. Though the step-lines are straight and parallel their directions 
do not apparently correspond to low crystallographic indices. 


(3) Cross-Laced Patterns on Cadmium Iodide 


Fig. 27 (Plate IV) shows a phenomenon similar to one observed by 
Verma (1951: fig. 13). Disregarding for the moment the supplementary 
pyramid on the lower slope of the main pyramid, it is seen that the 
terraces are ‘ cross-laced’ along each pyramid edge. Frank (1951 a) 
gave an explanation of this case. The slowest growing monolayer in 
the stack of monolayers which comprise one repetition of the crystal is 
a different one for two different orientations 60° apart. In each 
orientation, succeeding monolayers pile up behind the one with slowest 
advance, and since this is a different one in adjacent pyramid faces, the 
pyramid edges are cross-laced. Conditions for the same theory to apply 
can exist in a cadmium iodide crystal having a suitable ‘ polytypic ’ 
variant of the standard structure. It is difficult to see how this behaviour 
could be accounted for without polytypism. Further evidence for 
polytypism in cadmium iodide will be presented in Paper IT. 


(4) The Origin and Shape of the Big Screw Dislocation Groups 

As already reported (Forty 1951), the step-heights of the growth 
layers in cadmium iodide are often large : those in the photographs are 
-between 500 A and 1000 A. This indicates that the crystals are. growing 
from screw dislocation groups of large Burgers’ vector, and that groups 
of practically any strength are possible. 

A possible explanation of the origin of these large screw dislocation 
groups has been put forward by F. C. Frank (1951 a). The non-uniform 
distribution of impurities should be sufficient to stress the extremely 
thin plates beyond the theoretical yield stress for the ideal structure. 
A ‘buckle’ followed by slip of the sandwiches (close-packed layers of 
cadmium ions between close-packed layers of iodine ions) past one 
another in the [0001] direction, produces the large os CEELO EHO 
groups. The plate immediately begins to thicken by ‘ spiral growth 

rmed. 
wy eaten the growth of pure lead iodide crystals (these have 
a structure similar to that of cadmium iodide) show that only very nee 
‘screw dislocations occur. Screw dislocation groups of larger ee 
vector can be produced by adding a small amount of camium iodi e as 
impurity in the solution. This indicates the importance of impurities 
in the formation of the groups. The results of a study of step-heights 
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on growing cadmium iodide crystals, given in Paper II, provide some 
confirmation of the ‘ buckle, followed by slip ° process. 

It is somewhat surprising that such large dislocation groups are 
possible, since enormous strain energy in the crystal lattice might be 
associated with them. Frank (1951 b) has suggested that since the unit 
Burgers’ vector for carborundum is greater than 10 4, the core of a screw 
dislocation in carborundum is a hollow tube about one or two microns 
in diameter. This effectively relieves the strain energy around the 
dislocation. He expects large dislocation groups in carborundum to 
have a hollow core consisting of many such tubes. It is only occasionally 
(e.g. figs. 11-13 (Plate I) and 25 (Plate III) ) that there are indications of 
a hollow core for dislocations in cadmium iodide. 

Figs. 14-18 (Plate IT), 19-23 (Plate II) and other photographs taken 
throughout the growth of this crystal show a persistent line about 3 microns 
long at the centre of the spirals. If the dislocation group of total strength 
750 a is distributed along this line, making a ribbon of dislocation passing 
through the crystal, the portions of the lattice on either side of this ribbon 
are rotated only 14 degrees with respect to each other. In this arrangement. 
the strain energy density would not be excessive. 


(5) Effect of Dominated Dislocation Groups on Growth Fronts 
(a) The ‘ Cross-linking’ Type of Interaction 
Since the large screw dislocation groups can have practically any size, 
more than one group of exactly the same strength is rarely found on 
a given crystal face. For this reason, the ‘cross-linking’ type of 
interaction of advancing step-lines with a dominated dislocation group, 
whose Burgers’ vector equals the step-height of advancing layer, is 
seldom observed. (This is the normal type of interaction on Beryl where 
all the screw dislocations have unit Burgers’ vector (Griffin 1950, 1951).) 


(b) The Formation of Kinks in Step-lines* 

A step-line of height equal to the normal component of the Burgers’ 
vector extends in the surface from each point of emergence of a screw 
dislocation group. The growth of a particular face is normally dominated 
by one or a few groups. The step-lines from a dominating growth centre 
develop kinks where they pass over other dislocation groups of different 
strength. 

What can be called the fundamental crystal structure of cadmium 
iodide (G. Hagg 1948) consists of a layer lattice in which the unit cell 
contains two sandwiches each consisting of a layer of cadmium ions 
between two layers of iodine ions. Disordered variants of this structure 
have been reported. If the crystal has the fundamental cadmium iodide 
structure of two sandwich layers and the Burgers’ vectors of all 
dislocations present are multiples of the fundamental period, the 


dislocations are perfect and the crystal structure is continuously ‘ good ’ 
SE Ee ee eee 
* I am greatly indebted to Dr. Frank for these explanations of the kinking. 
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except in the dislocations. Some, however, may have Burgers’ vectors 
corresponding to an odd number of minimal sandwich layers. These may 
be called ‘ odd ’ dislocations for brevity. These are imperfect dislocations 
if the dominating dislocation group is even, so that the resulting structure 
is the fundamental one. If the dominating dislocation group is odd, the 
resulting structure is a polytype (see Paper II for the discussion of 
polytypism) and all the other dislocation groups, odd or even, are almost 
certainly imperfect. A surface of lattice discontinuity extends from every 
imperfect dislocation. Step-lines terminating on imperfect dislocations 
therefore mark fault surfaces and cannot move as easily as steps 
extending from perfect dislocations. 

A geometrically possible configuration for the interaction of a smaller 
dislocation group with dominating growth layers is shown in fig. 1 (a, b). 


Lata yl 


Tt can easily be shown that such a pattern can only be recurrent (in the 
sense that it reproduces itself after a cycle of rotation of the dominating 
spiral) if the approximately radial steps advance slower than the main 
spiral step. This can be the case if, and only if, the dominated dislocation 
is an imperfect dislocation of the structure resulting from the dominating 
one, so that the radial step approximately follows the trace of a fault 
surface in the crystal. Depending on the degree of misfit at us fault, 
various degrees of hindrance in the passage of a growth step over it may be 
anticipated. In the extreme case of complete hindrance we would 
expect a figure such as 2 (a). Then any excess of supersaturation on 
either side of the fault line will lead to different rates of advance of the 
main steps on either side, producing configurations like fig. 2 (b). This 
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Fig. 2 


(a) 


in its essentials, is the commonest form of interaction observed, 
(cf. figs. 5-10, 19-23, 26). (Plates I, I, IT.) 

. If the main growth layers dominate a pair of co-operating screw 
dislocation groups, other forms of interaction are possible. Fig. 3 
represents the formation of a possible configuration if a pair of smaller 


Fig. 3 
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co-operating screw dislocations of opposite hand are dominated. If the 
rate of growth of a composite layer is governed ‘by its lowest member, 
then it is likely that the rate of advance of step EF will be different 
from that of AB, CD although these sections have the same step-heights. 
This leads to a ‘fan-shaped’ region, radiating from the dominated 
dislocations, inside which the steps advance at a rate different from those 
over the rest of the hill. If the strength of the dominated dislocation 
groups is small, the steps such as GE might not be observed but the 
‘fan * will be obvious from the two rows of kinks radiating from the pair 
of dislocations. It is thought that this type of interaction occurs in the 
growth hill photographed in fig. 27 (Plate IV) where the two dominated 
dislocation groups are close together. Measurements have shown that 
the thickness of growth layers inside the fan-shaped region is equal to 
that of the layers on the rest of the hill. ' 

The interaction described in fig. 4 can occur if the directions of the 


Fig. 4 
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dominated pair of screw dislocation groups are opposite to those in 
the case described above. A thinner layer formed by the co-operation 
of the pair of smaller dominated dislocation groups is ‘ pushed ° ahead 
of each dominating growth layer. In this case there is not necessarily 
a connecting step between the main growth layers unless the rate of 
advance of a layer such as ADB is greater than the rate of advance of 
a layer ACB. This form of interaction appears on the crystal 
photographed in figs. 14-18 (Plate IT) where the main growth centre 
dominates the pair of small co-operating dislocation groups on its left. 


(6) Kinking of Step Lines due to Growth from Two Unequal Co-operating 


Dislocation Groups | 
Usually, when two co-operating dislocation groups are large, they do 
not have exactly the same strength. Hence, a step equal to the difference 


80 A. J. Forty on the 


of the Burgers’ vectors extends from the larger dislocation group to the 
edge of the crystal. The step-lines of the growth layers propagated from 
the two co-operating dislocation groups, usually become kinked across 
this line in a manner essentially similar to that of the domination of a 
single smaller dislocation group by another distant group. The same 
explanation holds in principle for the formation of rows of kinks from 


growth centres of any number of dislocations, of either hand. 


Figs. 14-18, 19-23 (Plate II) are photographs illustrating this type of 
kinking for growth from two unequal dislocation groups of opposite hand. 


(7) The Penetration of the Screw Dislocation Group through the Crystal 


The ‘ double-transmission ’ method is particularly useful for observing 
both surfaces of the crystal plates at the same time. It is found that the 
screw dislocation groups penetrate the plates, to emerge on both surfaces. 
If growth occurs on only one surface of a plate containing a single screw 
dislocation group a spiral step-line is observed on that surface and a 
radial step-line on the other (see fig. 25 (Plate III)). If growth takes 
place on both surfaces, a spiral step-line is observed on both, centred 
on the same screw dislocation group. Fig. 26 (Plate IIT) illustrates such 
a case ; the rates of growth of the two surfaces are obviously different. 
indicated by the difference in spacing of the step-lines. 

The cadmium iodide structure is polar (see Paper II). Upper and 
lower surfaces of the plate differ from each other and both have trigonal 
symmetry. Hence with similar degrees of supersaturation on the two 
surfaces it would not be surprising to find that the trigonal shapes of 
the polygonal spirals, centred on the same dislocation group, were rotated 
through 60 degrees of arc with respect, to each other. Actually, it is 
always observed that the rotation is (60-+-6) degrees, where @ is a small 
angle of 5 or 10 degrees of arc. This demonstrates that the segments 
of the polygonal growth spirals do not have orientations corresponding 
strictly to simple crystallographic indices. A proper understanding 
of the phenomenon has not yet been reached. 

According to Eshelby and Stroh (1951), the repulsion between like 
screw dislocations ‘in a crystal plate falls off rapidly at distances greater 
than the thickness of the plate. Therefore, if the large dislocation groups 
disintegrate by internal repulsions, the components should not move far 
apart. Usually, the group remains as a unit through the plate (by the 
internal interference method described in Paper IT, the step-lines of 
the two surfaces have been shown to be of equal height). However, in 
a few cases, growth on a single screw dislocation group produces a single 
spiral step-line on one’ surface, and two or more smaller, radial or spiral 
step-lines on the other surface. This could be interpreted as meaning 
that repulsion between parts of the screw dislocation group causes a 
separation which increases through the plate. This separation should 
become more pronounced as the crystal thickens, but this has not yet 
been observed. 


ee ree er 
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ABSTRACT 


The mode of decay of 87Rb is investigated chiefly by means of a shielded 
proportional tube operating as a spectrometer. It is shown to consist 
simply of a -transition between the ground states of **Rb and *’Sr the 
limiting energy of the particles being 275 kev and the average energy 
44 kev. The spectrum, which was investigated using sources of thickness 
1 mg/cm? and 0-128 mg/em?, shows no maximum, differs very 
considerably from that for an allowed transition and is in agreement 
qualitatively with that for a third forbidden transition with parity change, 
AI=3, yes. The period is found to be 7=6-15X10!° yr. No y-rays, 
X-rays or internal conversion electrons ascribable to the decay process 
are emitted in appreciable intensity. It is suggested that the established 
natural radio-elements with Z<80 form a well-defined and isolated group 
with log ft.=18 approximately. 


INTRODUCTORY 


THE radioactivity of rubidium has been the subject of many investi- 
gations since its discovery by J. J. Thomson in 1906 yet the nature of the 
decay was far from clear at the beginning of the present study. Thomson 
established the fact that the particles were negatively charged and later 
Campbell and Wood (1908) showed that the radiation was heterogeneous 
and less penetrating than that of potassium. The results achieved by 
more modern techniques are summarized in table 1 which gives the source 
thickness in mg/cm*, the half-life in 101° yr., the end-point energy Ey, — 
and the gamma ray energy hv in Mev. 

The period z is very important in connection with the age of the earth 
and mineral deposits and separate studies have been concerned with its 
evaluation. Hemmendinger and Smythe (1937) established the identity 
of the radioactive isotope of rubidium, viz. §*’Rb by isotope separation 
while Hahn, Strassman and Walling (1937) showed that strontium was 
the decay product, Mattauch (1937) identifying the particular isotope, 
*"Sr. Strassman and Walling (1938) estimated the period by the uranium- 
lead method, giving the value 6-3 1019 yr. 

Table 1 indicates the considerable uncertainty and disagreement in 
the results with the most consistent data appearing for +. It has to be 


* Communicated by the Authors. 
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noted that Muhlhoff could not allow for the abundance of the radio- 
active isotope. Haxel et al. consider that Eklund’s value of 7 is fort- 
uitously nearly correct. The same authors stress the importance of 
agreement in 7 values as obtained by purely radioactive measurements 
and the uranium-lead method. Anticipating our own findings we can 
eliminate their suggestion that the source emits two electrons (one decay, 
type and one from photo-conversion) per disintegration, so that it is 
desirable that their value of + should be confirmed or amended. 


Table 1 
Author mg/em?2 T Ey hy Remarks 
~S is) 
Muhlhoff, 1930 1-5 to 2} ~1+1 a Comparison of ab- 
3-0 sorption with Rak 
Klemperer 1935 15to | — 0-25 — Examination of 
3-0 Muhlhoff data 

Libby and Lee g — 0-132 — Screen-wall 

1939 : counter 

Ollano 1941 thick | — | ~0-13 0-034 |Gamma rays de- 
0-053 | duced from conver- 
0-082 | sion photoelectrons 
0-102 
0-129 

Saha 1946 300 — 0-144. = Modified Libby 

device 
Pepnd 1946 a 0-1 5:8 — — | Counter with inter- 
nal source 

Fazzini and 

Franchetti 1948 2 == 0:56 — Cloud chamber 

Haxel and A ~0-1 6-9 | ~0-01 ~0:-135 | Coincidences in 

Houtermans, 1948 G.M. tubes 

Haxel, Houtermans}| ~0-1 6:0 ~0-01 ~0:135 | Tentative alloca- 

and Kemmerich tion of decay B rays 

1948 and electrons 


Kemmerich 1949 0-03 6-0 — -+ Study of effect of 
% . source thickness 


From the stand-point of beta decay theory the transition s7Rb>878r+ Bo 
is of special interest. The spins of the grounds states of both nuclei 
are known (8*Rb=3/2, 87Sr=9/2, Jackson 1933, Heyden and Kopfermann 
1938, Millman and Kusch 1940) giving a spin change 4/=3 in the ground 
state to ground state transition and the calculated ft. value (Feingold 


GZ 
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1951) is 16-5 corresponding to a third forbidden transition. This probably 
unique case offers a searching test of theory, as Wu (1950) has already 
pointed out. 

Our experience in the past with large proportional counters operating 
as B- and y-ray spectrographs assured us that the spectra of sources of 
moderate energy could be analysed successfully even when the specific 
intensity was extremely low. The large acceptance angles, 47 for internal 
sources, 27 for wall mounting make the instrument ideally suited to 


investigations of weak radioactivities such as that of *7Rb. A large | 


source area can be utilized, ~10* cm?, so that a mass of 1 g say yields 
a source thickness of 1 mg/cm?. The “ background’ spectrum of a large 
counter at high pressure (5 to 7 atm) is spread out more or less uniformly, 


Fig. 1 


__ ALUMINIUM 
WINDOW 


~—PROPORTIONAL 
COUNTER 


~~~ SOURCE 
SUPPORT 
2" LEAD 


Proportional counter with lead and anti-coincidence counter shielding. 


over an energy range ~0-5 Mev. Hence for a K-capture source, whose 
radiation is concentrated in a range of less than say 5 kev, the sensitivity 
of the tube is ~100 times better than that of a Geiger-Miller. The same 
is true for sources decaying by emission of y-radiation or x-radiation 
and internal conversion electrons. For normal f-decay sources the 
position is somewhat altered but generally very favourable. Frequently 
the radiation from weakly active sources is of low energy and it covers 
therefore a relatively small part of the total energy range of the instrument. 
Within this restricted range the ratio of counting rate with source 
present to rate in its absence is enhanced. For example, if the 8-spectrum 
covers 50 kev the ratio is increased by a factor ~10 and near the maximum 
of the £-spectrum the factor is very much greater. 

The present investigation is part of a general study of the weak natural 
radio-elements of Z< 80. 


ee 
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EXPERIMENTAL ARRANGEMENT 


The schematic diagram of the spectrometer is shown in fig. 1. A large 
proportional tube of diameter 5-5 in., provided with field nities to give a 
fully effective length of 10 in. (Cockroft and Curran 1951), is surrounded 
with a lead shield, thickness 2 in. An array of Geiger tubes, 2 in. dia- 
meter, operates in anti-coincidence with the proportional counter and 
serves to eliminate most of the background due to cosmic rays. For 
photographic recording of the spectrum the Geiger impulses, suitably 
shaped are fed to the ‘ black-out’ electrode of the cathode ray tube 
on which the amplified output pulses of the spectrometer are displayed. 
Alternatively the output pulses of an anti-coincidence circuit are used to 
quench the time-coincident pulses from the spectrometer tube before 
they are applied to an electronic ‘kick-sorter’. Both electronic and 
photographic analysis of the spectrum were employed in the present 
work. Fluorescence x-rays of silver were used to calibrate the spectro- 
meter, the x-radiation entering the counter when necessary through an 
aluminium window. The proportional tube was maintained at a pressure 
of 5 atm of argon+ 20 cm of methane. 

The source consisted of spectroscopically pure RbCl deposited on the 
inner surface of an aluminium cylinder which was inserted into the tube 
and served as cathode. Two thicknesses were employed, namely 
1-5 mg/cm? and 0-128 mg/em?, the former laid down smoothly with a 
fine brush as a hot water solution and the latter with the help of a specially 
prepared glass wiper, Both appeared uniform in thickness and care 
was taken to keep them within the effective length of the tube. 


RESULTS 


(a) Spectrum Analysis. With the thicker source the spectrum was. 
analysed at different gain-settings and the overlapping portions coincided 
smoothly to give the shape shown in fig. 2. The region of low energy, 
- below 40 kev was examined separately using the thin source (0-128 
mg/cm?) with the results shown in fig. 3. The curve thus obtained was 
fitted to the curve for the thicker source and smooth coincidence from 
40 kev to about 30 kev observed, as seen in fig. 2. 

Below about 20 kev there is clearly a serious absorption of 6-particles 
in the thicker source (curve indicated by dashed line). The final form of 
the complete spectrogram, full line in fig. 2, is believed to be accurate 
down to about 10 kev and over the whole range from 10 kev to the 
end-point at 275 kev no maximum appears on the curve. The curve 
near the end-point is shown in more detail in the inset figure. The curve 
falls very steeply at first with increasing energy and it shows an extremely 
long ‘tail’. No doubt the very unusual form of the spectrum explains 
much of the difficulty of earlier workers. 

The Fermi plot is shown in fig. 4. The curve departs very remarkably 
from the straight line form given by an allowed transition. The extra- 
polation of the Fermi curve in the vicinity of Ho, assuming it linear over 
a short range gives H)=275 kev in agreement with direct observation. 
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(b) Gamma Radiation and Conversion. The non-appearance of 
conversion lines on the spectrum seemed remarkable in view of Ollano’s 
results and a sensitive method of searching for internal conversion was 
adopted. The thick source (1:5 mg/cm?) of RbCl was covered with 
sufficient aluminium to stop the B-rays (58 mg/em?), This cover reduces 
the intensity of K x-radiation of Sr to about 50%. The average efficiency 
of the counter, at 2 atm pressure, for such rays was assessed at 84%. 


INTERVAL 


Fig. 2 
l400 800 
1200 600 
IOOO 400 
800 200 
600+ ‘< fe) 
po ISO 200 pif 250. 300 SSS Cine 


400 


200 


NUMBER OF PULSES PER ENERGY 


O 50 lOO ISO 200 220 Be OO RN SOR cay 
B-spectrum of §*Rb. 


The results of our search are shown in fig. 5. The total number of X-rays 
emitted per unit time was determined from the histogram and comparison 
with the B-rays from the source showed that the x-rays were down in 
intensity relative to the B-rays by a factor of more than 500. The 
‘ extremely low intensity of x-radiation made accurate definition of energy 
difficult but as indicated in fig. 5, allocation to Sr seems slightly preferable 
to Rb. On the other hand the intensity is not inconsistent with the view 
that the x-rays are generated by bombardment of the rubidium by the 
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Fig. 3 
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p-rays. In any case it appears certain that internal conversion in the 
K-shell of 8’Sr does not occur in more than 0-2% of the disintegrations. 
This result completely contradicts the conclusions of Ollano and we are | 
forced to the view that his electron groups are spurious and the result 
of unfortunate statistical fluctuations. ; 

A separate search for L x-ray quanta was carried through since there 
was a remote possibility that y-radiation might convert in the L shell 
only. The low energy of L-radiation (~ 1-8 kev) prohibited the use of 
the B-absorbing screen of aluminium. Using x-rays of Ga (9-2 kev) 
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emitted from “Ge, the threshold energy for detection, as defined by a 
discriminator, could be varied over a wide range; above and below the 
value of 1-3 kev. Now if X-rays are emitted in some disintegrations 
in 50% of such cases the x-ray will not be detected in coincidence with 
a f-particle. Hence a relatively sharp rise in counting rate, could be 
expected as the threshold energy is reduced through the value 1-8 kev. 
The results of fig. 6 show that no appreciable discontinuity can be found. 
Even if the slow rise on the curve is distorted to create a discontinuity 
as shown by the dashed line on the figure an upper limit of about 3% is 
found. It is much more likely however that the upper limit is not prentee 
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than that for the K-radiation namely 0-2°/. Thus there is no real evidence 
of conversion of y-radiation in either K or L shells in more than 0-2°% 
of the disintegrations. i 

Finally, a sensitive search for y-radiation was carried out by Mr. 
Bannerman of this Department, with a scintillation counter operating 
as a spectrometer. About 5-6 g of RbCl was packed round the scintil- 
lator (Nal, Tl) and the counter detected an additional 10 counts/min. 
The crystal was shielded from the B-rays with aluminium. There was 
no indication of homogeneity in the pulses produced and they were almost 
certainly excited by the passage of the f-rays through the absorber- 
bremsstrahlung. Whatever the origin we can conclude that y-quanta in 
the range 20 to 100 kev cannot be released in more than one per 5x 103 
disintegrations and in more than one per 10% disintegrations in the range 
100 to 300 kev. 

Fig. 6 
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Search for L x-rays. Counting rate as discriminator bias is varied over range 
including the L-radiation. 


It would appear certain from this search for y-quanta, K and L X-rays 

that the probability of the B-decay leaving the nucleus *’Sr excited is 

negligibly small. The decay process consists in the transition by 6-decay 
between the ground states of ®7Rb and §*Sr. 


THE PERIOD 7 
The conclusion that the B-decay of *7Rb is simple raises some issues 
concerning 7. The most recent accurate work of Haxel, Houtermans and 
Kemmerich (1948, 1949) was based on false assumptions regarding the 
decay scheme since they assumed that two electrons were emitted in each 
disintegration. Clearly the unusually high intensity of the continuous 
spectrum at low energies explains to some extent their observations, but 
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we are forced to suggest that the coincidences observed in their double 
counter probably had their origin in reflexion at the counter walls. 
However, as they themselves pointed out, their method of deducing the 
number of disintegrations gives a result which is independent of their 
assumptions regarding the nature of the coincidences. Indeed in our 
view it amounts to correcting for reflexion by materials near the source. 
We were interested to compare our value of 7 with theirs. 

Accurate measurements of the counting rate at high gain showed that 
the rate was practically independent of the voltage over a range of some 
hundreds of volts. The difference count with and without the source 
(0-128 mg/cm?) was 1477 counts per min for the total mass of 0-10075 g 
of RbCl, measured to a statistical accuracy better than 1%. Allowing 
for the 27 geometry, chemical composition and isotopic abundance this 
gives directly 7=5-95 x 101° yr. 

Two corrections have to be applied for (a) the reflexion of particles in 
the source and support and (b) the self-absorption of the source. To 
estimate (a) we made use of data due to both Yaffe and Burtt (1950). 
These workers are in fairly close agreement but both plot the variation 
of reflexion coefficient as a function of Hy. Due to the remarkable shape 
of the spectrum of ®’Rb it is probably more desirable for us to estimate 
the reflexion in terms of the average energy, H of the f-rays. Thus we 
find # for 87Rb is 44 kev, less than the values 52 and 54 kev for “C and 
3°S respectively although Hy for **Rb considerably exceeds Hy for both 
4C and *58. In the region involved here the reflexion seems to vary 
rapidly with the energy of the £-rays, further complicating the estimation. 
The value which we adopt is 7-594. Regarding (b) we take the source as 
equivalent in thickness to the range of a B-particle of energy 5 kev. 
We deduce from fig. 2, extrapolated portion between 10 kev and zero 
energy, that the self-absorption reduces the counting rate by 45%. 
Hence the overall contribution to the counting rate conists of an increase 
of (7:-5—4-5)%. The final value for 7 is therefore 5-95x1-03 or 
6-15 101° yr. The chief uncertainties in this value arise from the 
estimation of the corrections and together with statistical fluctuations in 
counting we consider an error of +4°% is possible. With these rather 
conservative limits our value for tr becomes t=6-15-+0-3 x 101° yr: 


DiscussioON AND CONCLUSIONS 

The experiments establish that the -decay of 87Rb is a simple process 
involving transition from ground state to ground state with a maximum 
energy Hy=275 key. The Kurie plot lies very markedly below the 
straight line corresponding to an allowed transition and seems consistent 
with 4J—3 with change of parity according to present theoretical views. 
The calculation of the correction factor Cy, given by Greuling (1942) 
and containing an adjustable factor (the ratio of the matrix elements 
Q3(Ba, 7) and Qs(Bo x7, 7) in Greuling’s notation) is exceedingly laborious 
but it is probable that adjustment can be made. We have not attempted 
to complete this task. 
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The degree of agreement of our 7 value with the most recent measure- 
ment and the geological estimate gives considerable confidence. 

The new value of Hy, 275 kev leads to the value 17-6 for log ft. A 
remarkable degree of agreement in ft. value is now apparent among four 
of the best established natural radioactivities as shown in table 2. Now it 
is to be noted that a considerable gap exists between the value 17-6 for 
*"Rb and the next well-established value 13-65 for 1°Be although there 
is a fairly uniform distribution below this latter value (see the data due 
to Feingold 1951). The Table shows °Nd as an obvious exception on 
the assumption that H)=11 kev and r=5 x10! yr. Some preliminary 
work with the equipment already described has shown that these values 
(Libby 1934) are almost certainly erroneous* since the necessary number 
of low energy f£-rays could not be detected with the counter. The period 
may be underestimated but assuming that it is very approximately correct 
and that H is in error, the calculation shows that log ft for °°Nd can be 
made approximately equal to that for the other elements of the group 


Table 2 
Source Ey 1020 vr. log ft. 
MeV 
Rbé? 0:275 ~ 6:15 17-6 
oo 1-36 0-11 18-05 
buys 0-215 2-4 18-02 
0-4 or 18-91 
Re187 0-043 4-0 17-73 
iINd222 0-011 5 13-73 


when H,—250 kev (see Feenberg and Trigg 1950 for Log f as function of 
energy and Z). Further work is contemplated but at present it appears 
that such a value is not entirely ruled out by our results. If this should 
prove to be the case the natural radio-elements with Z7<80 will form a 
narrow group with log ft close to 18. 


We wish to record our gratitude to Professor P. I. Dee, F.R.S., for his 
interest and encouragement in this work. 
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SUMMARY 


Tension, compression and rolling textures have been derived for the 
orthorhombic metal «-uranium. From a consideration of the relative 
importance of the slip and twinning mechanisms it is possible to 
differentiate between the textures developed after slight and severe 
deformations at both room and elevated temperatures. 


§1. INTRODUCTION 


A METHOD has recently been proposed by which the preferred orientation 
developed in metals as a result of various mechanical operations may 
be derived from a knowledge of the deformation modes of their single 
crystals (Calnan and Clews 1950, 1951 a, b). The treatment has been 
applied successfully to the prediction of tension, compression, rolling 
and drawing textures of face-centred cubic, body-centred cubic and 
hexagonal metals, and the purpose of the present paper is to extend it 
to the orthorhombic structure, «-uranium. For this metal both the 
detailed crystal structure (Jacob and Warren 1937) and the deformation 
modes (Cahn 1951 a) have been established. 


§2. FUNDAMENTAL PRINCIPLES 


The slip and twinning systems observed by Cahn are listed in table 1. 
’ He concludes that at room temperature twinning is the predominating 
- mechanism, the relative importance of the twinning systems being in 
the order listed, and that at temperatures above about 300° c (Cahn 1951 b) 
slip plays a larger part. The slip system (a) is by far the more important. 
Since the amount of deformation associated with twinning is limited, 
considerable slip must occur at all temperatures if high deformations 
are attained. : 

Consider now an isolated single crystal of uranium deforming only on 
the slip system (a) under a single tensile stress whose orientation is 
represented by the point, P on the standard stereographic projection, 
fig. 1. The resolved shear stress on the most favourable slip system, 
namely the direction [100] on the plane (010), varies with the orientation 
of the stress direction in the manner illustrated by the contours in this 
figure. Thus the resolved shear stress corresponding to an applied stress 
direction P is 0-3 of the applied stress, 7',, and slip occurs when 0°37’, 


* Communication from the National Physical Laboratory. 
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reaches a critical value. In a polycrystalline aggregate, however, slip 
will not generally occur when this condition is satisfied owing to the 
constraints imposed by the neighbouring grains. This implies that the 
critical value has not been attained on this slip system and that the 
effective stress, 7',, may be represented by a point on a lower contour. 


Table 1. Slip and Twinning Modes in Uranium 


Direction 


a (010) [100] 
b (110) 2 
Twinning nies Ie or Ne 
¢ (130) |. (110) [310] [110] 
d eam (112) [312] | Irrational 
€ (112) Irrational | Irrational [312] 
~(172) 


100 


Standard stereographic projection of a-uranium showing resolved shear 
stress contours for slip on the (010)[100] system. 


T, may continue to move down the contours until it reaches the 
(001)—(010) boundary where the resolved shear stress on either of the 
(0105 <100) systems is zero, after which no further movement of 7’ ‘ 
possible and no slip occurs. The resolved shear stress for stress are 
orientations on the (001)—(100) boundary also is zero. Thus it will ce 
seen that no duplex or multiple slip is possible and that slip on one 
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system only with associated single crystal rotation can occur. This 
rotation moves the stress axis orientation along the great circle towards 
the slip direction [100] and gives a [100] tension texture. 

The re-orientation associated with twinning on a single system, for 
example (c), will now be discussed. Fig. 2 shows the twinning elements 
of this system and the macroscopic reference circle of the original crystal 
which shears on the plane K, in the direction y, to become an ellipse. 
The plane K,, whose trace is the line OB, N2, is undistorted on twinning 
by movement to the position OD. It can be shown that the change of 
orientation due to twinning is crystallographically equivalent to rotation 


Fig. 2 
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TRACE OF K,, (130) 


TRACE OF 
Ky, (110) 


REFERENCE CIRCLE 


inni i 1 lisplacement. 

Geometry of twinning. Plane of diagram is plane of shear disp 

Ky ea K, are the first and second undistorted planes. 1, is the shear 

direction, 7 is the intersection of the second undistorted plane on the 

plane of shear displacement. The position of K, is that within the 
untwinned crystal (Schmid and Boas 1950). 


through 180° about the twin plane normal although this corresponds to 
neither the macroscopic nor the atomic movements. For simple tensile 
stresses only those which give a positive resolved shear stress on 
the plane K, in the direction 7, will produce twinning, that is to say, 
those within the quadrant COE. Thus in the stereographic projection, 
fig. 3, twinning under tension can occur for those grains only which have 
the stress direction within the quadrant (130)(001)[310], corresponding 
to the quadrant COE of fig. 2. The sector (100)(001)[310] may be 
neglected as here the system K ,=(130) is the more favourable. In order 
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to determine the twinned orientation of a grain or part of a grain whose 
stress axis is represented by the point P, fig. 3, a great circle is drawn 
through this point and the shear plane normal, K,. The angular distance 
between P and K, is measured and the new orientation of P will be at 
an equal distance along the great circle on the opposite side of Ay, 
ie. at the point P’. Thus the twinned orientations corresponding to the 
area (130)(100)(001) lie within the region T 39 1499(001). 


§3. TENSION TEXTURES 

The re-orientation with twinning on system (c), A,=(130), has just 
been derived, and that for the system (d), K,=(172), is illustrated in. 
fig. 4. Orientations with the stress axis originally within the area 
(100)(001)QR twin to the areas TyyTpT,ToTp and To;TpT,. The 
re-orientations for the reciprocal system (e) are very similar and may 
be considered as augmenting system (d). For conditions where twinning 


Fig. 3 
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Standard projection showing the twinning elements for the system 
K,=(130), the orientations for which twinning occurs under tension 
(heavy outline), and the orientations after twinning (heavy broken line). 


predominates the texture will be determined principally by twinning 
re-orientations. Superposition of figs. 3 and 4 as in fig. 5(a) shows that 
for orientations with stress axes within the region (100)(001)C twinning 
may occur on either system (c) or system (d); in the region (001)BC on 
the system (d) alone, and in the remaining area of the unit quadrant 
(010)BA no twinning can occur. The twinned orientations may be 
obtained similarly by superposition within the reference unit quadrant 
also shown in fig. 5(a). The operation of both these twinning systems 
removes orientations from the area (100)(001)EFG to the area 
(010)(001)EFG which thus becomes a region of preferred orientation. 
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The density of orientations is not constant throughout this latter area 
but is highest in the region (010)BH where the superposition of twinned 
orientations is greatest. Subsequent retwinning leads progressively into 
the area (010)BA from which no twinning can occur. This last region is 
therefore the final texture which could be attained by tension twinning 
alone on the assumption that such severe deformation could be reached 
without fracture, fig. 5 (b). 

The rotations for the slip system (a) are those associated with single 
crystal slip in the direction indicated by the arrow along the great circle 
in fig. 5(c), and the preferred orientation resulting from twinning and 
some slip is shown by the cross-hatching. As the slip system (b) has 
not been fully established and, moreover, is believed to be of very minor 
importance it has been neglected in this analysis. Further slip rotation 


Fig. 4 


100 
Twinning elements for the system K,=(172) and re-orientations under 
tension. 


does not, however, lead to a (100) texture since in the neighbourhood of 
(100) twinning is likely to be the more favourable deformation mode. 
The position of the boundary between the regions of slip and twinning 
depends upon the relative critical shear stresses for these modes. In the 
absence of an quantitative data on the magnitudes of these stresses the 
boundary (001)Z has been drawn in fig. 5(d) for a ratio of unity. Thusa 
grain of stress axis orientation near (010) rotates by slip towards (100) 
but on crossing the (001)Z boundary twins back to the neighbourhood 
of (010) from whence it may again rotate by slip. This cyclic process 
moves the highest density of orientations to the vicinity of Z, i.e. near 
the (110) pole and maintains a spread of orientations from Z to the (010) 
end of the (010)(001) boundary. In this was the texture corresponding 
to twinning and a considerable amount of slip is derived as shown by the 


cross-hatching in fig. 5 (d). 
SER. 7, VOL. 43, NO. 336.—JAN. 1952 H 
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§ 4. COMPRESSION TEXTURES 


The compression textures obtained by twinning and varying amounts 
of slip are derived in a manner precisely similar to that just described. 
The orientations of stress axis for which twinning may occur are seen 
by reference to fig. 2. In order to produce a positive resolved shear 
stress in the direction 7, on the plane AK, the simple compressive stress 
must lie within the quadrant COA. Thus for the system (c), table 1, 
twinning will occur on the (130) plane in the shear direction [310] for 
stress axis orientations within the quadrant (130)(001)[310], fig. 6. The 
sector (130)(001)(010) may be neglected since in this region twinning on 
the symmetrical system K ,=(130),7,=[810] is more favourable. ‘Twinning 


Fig. 5 


000i 


100 
Sores ae (@) 
Tension textures, cross-hatching shows density of orientations. 
(a) Twinning re-orientations in the unit quadrant for both twinning 
systems derived by superposition of figs. 3 and 4. 
(6b) Texture resulting from twinning alone. 
(c) Texture resulting from twinning and some slip. Slip rotation indicated 
by the arrow. : 
(d) Texture resulting from twinning and a large amount of slip. 


from the remaining region (010)(001)[310] gives re-orientations such that 
the corresponding stress axis positions are within the area [310](001)X. 
In fig. 7 are shown the re-orientations for the system (d), K,= (172). 
Here twinning takes place from the area (010)YZ to the area lei6k ae 
Combination of figs. 6 and 7 shows that the operation of both systems 
removes orientations from the area (010)Y[310] in fig. 8 (a) to the area 
(100)(001)Y[310]. The overlapping of twinned areas is greatest in the 
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apie ABC and thus the compression texture corresponding to deformation 
y twinning alone is that shown by the density of cross-hatching in fig. 8 (b). 


ope 
--- See 


LJ 100 


Twinning elements for the system K,=(130) and re-orientations under 
compression. 


Fig. 7 
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Z 100 
‘Twinning elements for the system K,=(172) and re-orientations under 
compression. 


The direction of rotation resulting from slip on the system (q) is towards. 
the pole of the operative slip plane (010) indicated by the arrows in fig. 8 (c).. 
A small amount of slip in addition to the twinning will produce the 


HZ 
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texture shown in this figure. For twinning together with a considerable 
amount of slip the texture will be as in fig. 8 (d), where again the slip-twin 
boundary (001)X has been drawn for equal critical shear stresses on the 
slip and twinning systems. 

Fig. 8 
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Compression textures. 
(a) Twinning re-orientations for both systems derived from figs. 6 and 7. 
(6) Texture for twinning alone. 
(c) Texture for twinning and some slip. 
(d) Texture for twinning and a large amount of slip. 


Table 2. Tension and Compression Textures 


Plane normals parallel to Plane normals 
Tiatornation perpendicular to 
processes Tension axis Compression axis compression axis 
(tension (compression (perpendicular 
texture) texture) texture) 
Twinning (O10) Zone (100)-(001) (O10) 
Twinning+ (O10) — (130) Zone (100)-(001) (010) spread 
Some slip broadened towards (110) 
Twinning +- (110) ~(170) Zone (001)—(110) 


large amount 
of slip 
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While none of the tension and compression textures can be fully des- 
cribed in terms of ideal textures, their main features may be tabulated as 
in the first three columns of table 2. 


§5. Rotting TExTuRES 

Sheet Rolling. The compressive stresses normal to the plane of the 
sheet lead to textures in which the planes listed in column 3 of table 2 
lie in this plane. The texture for any direction perpendicular to the 
compression axis (perpendicular texture), such as the rolling direction, 
may be obtained by drawing on the standard stereographic projection 
narrow bands, bounded by small circles, whose intensity corresponds to 
that of their poles in the compression texture. Thus to the regions A, 
B and C in the compression texture of fig. 9(a) there correspond the 
bands in fig. 9(b).. It will be seen that with additional regions inter- 
mediate between A, B and C, forming a compression texture similar 
to fig. 8(b), there is a high density of orientations in the perpendicular 
texture in the neighbourhood of (010). This and the other perpendicular 
textures are listed in column 4 of table 2. The minor component of 
rolling, namely tension in a specific direction perpendicular to the 
compression axis, may then be introduced by applying the tension 
re-orientations to these textures. Comparison of columns 2 and 4 of 
table 2 shows immediately that the tension textures are consistent 
with the perpendicular textures, and therefore the re-orientations may 
occur simultaneously producing the compression texture parallel to the 
sheet normal-and the tension texture parallel to the rolling direction. 

Rod Rolling. The textures derived for the direction perpendicular to 
the compression axis, column 4, table 2, do not imply any additional 
spatial relation between these two directions, and are thus equally valid 
for an array of similar radial stresses. The rod may be envisaged as an 
aggregate of sheets in which each sheet normal is aligned parallel to any 
one of the radial stresses, and each rolling direction parallel to the common 
rolling direction. The preferred orientation of the rod axis is thus the 
texture derived for the sheet rolling direction, i.e. the tension texture. 
On this view the texture in any radial direction should be a random 
distribution of all plane normals perpendicular to the tension texture. 
If, as is more likely, the predominating influence is the compressive 
stress parallel to the particular radial direction under consideration, then 
some evidence of the simple compression texture may be found in this 
radial direction. In general, therefore, a longitudinal section may have 
a complex texture varying from point to point on a given section, and 
from one section to another. 


§6. Errect oF TEMPERATURE AND DEGREE OF DEFORMATION 
On the pasis of Cahn’s observations that the slip mechanism is 
relatively more important above 300° c it might be concluded that the 
room temperature deformation textures would be primarily those 


102 E. A. Calnan and C. J. B. Clews on the 


resulting from twinning alone, figs. 5(b) and 8(b). The amount of 
deformation attainable by twinning is, however, not unlimited, the 
limit depending theoretically on the shears associated with the twinning 
systems and the amount of successive retwinning which occurs before 
the final texture is reached. It may be estimated from a knowledge of 
these shears that this limit is likely to be about 30-40° deformation, 


Fig. 9 


o}fe) 


O10 


Nlustrating the derivation of the texture perpendicular to the compression 
axis for a given compression texture. 
(2) Shows the compression texture, and (b) the corresponding bands 
whose superposition gives the perpendicular texture. te 
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and therefore the room temperature textures for deformations appreciably 
greater than this will show some influence of slip, figs. 5(c) and 8(c), 
tending in the limit to the textures shown in figs. 5(d) and 8 (d). 

Above 300° c, as a result of both slip and twinning, textures similar 
to figs. 5(c) and 8 (c) may be expected in the early stages of deformation, 
while after severe deformation the textures of figs. 5 (d) and 8 (d) will be 
produced. 


| Table 3. Rolling Textures 
eae ee a a 
Elevated temperature 


Room temperature 
(> ~300° c) 


Deformation Pl F Plane normals ; Plane normals 
spain parallel to Planes in parallel to 
rolling rolling rolling rolling 
plane direction plane direction 
Low Zone (010) Zone (010) + (130) 
(100)-(001) (100)-(001) 
: broadened 
High Zone (010) — (130) ~(170) (110) 
(100)-(001) 
broadened 


The textures defined are those for sheet rolling. The textures for the axis 
in rod rolling are given by columns 3 and 5. 


These conclusions are summarized in terms of ‘ideal’ textures in 
table 3. The only experimental work with which these predicted 
textures may be compared is that of Harris (1951) who has studied the 
preferred orientation in rolled bars by Geiger Counter methods. He finds 
that light reduction by cold rolling gives a texture in which (010) plane 
normals lie parallel to the rolling direction with some spread of orientation 
towards (130). Light reduction by hot rolling gives a similar texture 
with a slightly greater spread of orientation towards (110), while a heavier 
reduction by hot rolling gives a strong (110) texture with a spread towards 
(010). These observations are in good agreement with the predictions in 
columns 3 and 5 of table 3 as is more clearly demonstrated by comparison 
of figs. 5(b), 5(c) and 5(d) with the figs. 4 (a), 2 (a) and 3 (a) respectively 
of Harris’s paper. Some refinement of the treatment will be possible 
when the relative magnitudes, at various stages of the deformation, of 
the critical shear stresses for the slip and twinning systems have been 
determined. On this ratio depends the position of the slip-twinning 
boundary discussed above, which in turn defines the limit of the texture 
when appreciable slip occurs, e.g. fig. 5(d). It has been assumed that 
this ratio is unity and the agreement between fig, 5 (d) and fig. 4 (a) (Harris) 
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suggests that this is indeed approximately so, at least in the later stages 
of deformation. The spread of orientation beyond (110) towards (100) 
may, however, indicate that the critical stress for twinning is slightly 
greater than that for slip. 
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ABSTRACT 


The stability of 6-brass, with respect to shears across the (110) plane 
in the direction [110], is considered and is shown to be due, mainly, 
to the change of the Fermi energy as a result of the shearing of the 
- Brillouin zone. Calculations are made which show that this effect 
contributes approximately 10" dynes/cm? to the shear constant 1(c,,—c,), 
_and a rather larger amount to ¢,,. The constitution of the elastic shear 
constants, as the sum of three terms, is examined and a satisfactory 
account of the observed values is obtained. The difficulty experienced 
in attempting to retain the B-phases by quenching is briefly discussed. 


$1. INTRODUCTION 


It is well known that the interaction between closed ionic shells, in metals 
or alloys with the body centred cubic structure, leads to a negative 
contribution to the elastic constant corresponding to a pure shear over 
the (110) plane in the direction [110]. In the alkali metals, because of 
the large atomic volumes, the ions are so well separated that this effect 
is very small and the body centred structure is stabilized mechanically 
by electrostatic interactions. In the other pure metals which have body 
centred cubic structures, e.g. W or Cr, the ionic shells are incomplete 
and these metals are rendered stable, according to Zener (1951), by 
exchange spin interactions. Alloys which are body centred cubic are the 
B-phases, occurring at electron-atom ratios of approximately 3/2, of 
which f-brass is typical. The distance between nearest atoms in B-brass 
is very nearly the same as that between nearest atoms in pure copper 
where it is known that the non-Coulomb ionic interactions very largely 
determine the rigidity. The question arises therefore: what is the 
stabilizing influence in £-brass which ensures a positive value for the shear 
constant 4(c,,—C,,) ? Isenberg (1951) in a recent paper, having in mind 
that the disordered B-phases are stable only at high temperatures and 
cannot, in general, be quenched writes ‘the stability of disordered 
B-brass under a (110) [110] shear is due in some way to thermal 
vibrations’. It is the purpose of this paper to examine in detail one 
possible cause, not directly concerned with thermal vibrations, of the 
observed stability of the 6-phases. 
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A calculation of the elastic constants of aluminium, made by Leigh 
(1951), shows that the change of the Fermi energy of the conduction 
electrons under a shear can lead to large contributions to the shear 
constants. In monovalent metals it is known that this effect is, in general, 
small because in these cases the Fermi surface which contains the occupied 
electronic states is not in contact with the Brillouin zone (BZ) 
boundaries. In £-brass of composition CuZn there are 1-5 electrons per 
atom to be placed inside the {110} zone which contains two states per 
atom. It may be assumed, therefore, that in this metal the Fermi 
surface will be in contact with the zone boundaries and, indeed, this was 
made the basis of a theory of the occurrence of this phase at an 
electron-atom ratio of approximately 3/2 (Jones 1937). 

In the following section we calculate the contribution, arising from the 
Fermi energy, to the two shear constants ¢c,, and }(c,,;—C,,) ; for the sake 
of brevity the latter will frequently be denoted by c’. The conclusion 
is reached that this contribution amounts to at least 1014 dynes/em? for 
the constant c’ in B-brass and therefore provides the explanation of the 
mechanical stability of the phase. In § 3 it is shown that the contribution 
from the Fermi energy along with the estimated electrostatic effect and 
the non-Coulomb interactions accounts in a satisfactory manner for the 
observed values of c,, and c’ as determined by Lazarus (1949). Finally, 
a reason is suggested for the difficulty experienced in attempting to 
retain the disordered B-phase by quenching. 


§2. THE Fermi ENERGY EFFECT 


In this section we shall calculate how the Fermi energy changes as the 
crystal is sheared at constant volume. We shall assume throughout that 
no electronic states outside the {110} BZ are occupied. If the energy 
of a state within the BZ depends upon the momentum in the same way 
as for a free electron, but the above assumption is maintained, then 
the contributions to c,, and c’ can easily be worked out exactly, as shown 
at the end of this section. The approximation of regarding the electrons 
as free for states within the BZ is, however, very unsatisfactory for B-brass 
and we shall take a more general dependence of the energy on the 
momentum into account. In order to carry out the calculation in this 
more general case, we have assumed that the effect on the Fermi energy 
of each pair of planes of the BZ can be calculated separately and the 
total effect obtained by summing over all pairs of faces. It is proved, 
at the end of this section, that for nearly free electrons this last 
assumption introduces negligible errors. 

The wave vector k specifying an electronic state is defined so that 
hk is equal to the momentum of the electron. The number of states per 
unit volume of metal associated with a volume element dk is thus 
(1/47*)dk, where account has been taken of the doubling of each state by 
the spin. If we measure energies in Rydberg units (Ry=21-79.10-12 ergs) 
and lengths in units of the first Bohr hydrogen orbit (4y=0-529.10-8 em) 
the energy « of a free electron in terms of k is given simply by «=k?. 
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If a denotes the lattice constant, the planes which form the boundaries 

of the BZ are given by 
Lape = (75)/0): | 0 |e peace eed ab Wael oo 841) 
where n is one of the set of vectors {110} each of magnitude 4/2. 

We shall now calculate the effect on the Fermi energy of a single pair of 
these planes (+n) on the assumption that the energy at any point k is 
given by 

e=kh*—k2+ p,*f(z), en oe ee. ok re (2) 
where p,= \/2n/a is the perpendicular distance from the origin of k-space 
to either plane, and f(z) is a non-dimensional function of z=k,/p,. The 
axes in k-space have been chosen so that k, coincides with the normal 
to the two planes +n. We need not, at this stage, specify f(z) in further 
detail except to notice that for small z it behaves like 2? and that 

S'(I)=0. 

When the energy has the form (2) the density of states function 

N(e) can be calculated immediately and we obtain 

N(e)=(1/42°) [| grad €|“dS=(1/27?)k vax; mene coe (3) 
where k,,,, is the greatest value of k, on the energy surface « (Jones 
1937). For values of e>e,, where <, is the energy of the surface which just 
touches the planes +n, N(e) is constant and equal to (1/27?)p,. 

The number of electrons per unit volume n is given by integrating 
N(c) from zero to the energy at the surface of the filled states ¢,,,,, and 
if €y3x>€, we obtain from (2) and (3) after an elementary reduction 


271?n/Da® = (maxi Pa®)— [i f(e)dz. rear aC) 
In a similar way we obtain for the Fermi energy per unit volume W 
4n*W |p P=(EmaxlPa®)?— | Pteldz. . - . - - (8) 


If we denote the integral in equation (4) by uw, and that in (5) by w,, and 
if we now eliminate (€,,,./Ppn2) between (4) and (5) we obtain 

4? W = (4274n?/Dy) + 477?2Uy Da? — (Ua— Uy") Pn’- ee Pee (6) 
Thus if we keep n constant, i.e. the atomic volume constant, (6) gives 
the change in W under a pure shear which alters the value of pa. 

Let W, denote the energy of free electrons at the same number density 
n, then (W—W,) can be regarded as the depression of the Fermi energy 
brought about by the pair of planes +n. Thus for the body centred 
cubic structure, if we sum (W—W,) over the six pairs of planes which 
form the BZ, the changes in this sum due to a shear should give a good 
approximation to the change in the Fermi energy under this shear. Thus 
we write for the Fermi energy W” 


W?—SW(p,)—5 Wo. aS eee (7) 
For a pure shear in which all the strain components are zero except 
€,, We find prank aay P 
; fe Ee ee yy (ee W's (55) we ee 3) 
oe ( a ie a) i Ge 0 
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where ¢,," denotes the contribution to the elastic constant C4, arising 
from the change in W”, and the subscript to the derivatives means that 
these are to be taken at ¢,,=0. 

It may easily be shown that to second powers of the strain component 


(a/7) 29 2=|n |?— 219. py+3(%1?2-+N9” Je reat = 2 . (9) 
and thus, the sums in (8) may be evaluated and we obtain 
Cy = hp? W" +8 WW. .. 2 ct. i 2 ee 
In a similar way, for the pure shear e,,=e, €,,=€,,,=(1+e) *—l1, we find 
ow Op 0" Dn 
1 nS 1! Sh fons F TH n ef - ’ “ LI 
Hen—ere)=3 (Ga) awe (Fe) +aw'2 (Ge)... 
where 
(a/z7)*p,2=|n |? + (m,2-+-_2—2ng”)e+3n,7e7, . . . . (12) 
and thus 
4(¢443—C1.) 7? = 4 p?W'+ Epw’. . . .. «. (13s) 


In equations (10) and (13) we have written p for the value of p, in the 
unstrained state, and W is given by (6). 

It will be noticed that (10) and (13) have the same form as the 
contributions to the shear constants, arising from central forces between. 
nearest ions, in a face centred cubic lattice (Fuchs 1936). The reason 
for this is evidently that the reciprocal lattice of a body centred cubic 
space lattice is face centred cubic, and dW/dp, behaves like a central 
force acting along the line joining the points +-n and —n in the reciprocal 
lattice. From equations (6), (10) and (13) we now obtain easily 

(1]t)egg” = —- $v? 4vu— 32 (wp—uy?), - ois a ee 

(1/k)3(Cyy—Cyo)" = — §v2+4vu,— BP (ug—u,”),. . - - (15) 

where x=(p>/27?)(Ry/ay?) and v= 27*n/p*. The number of conduction 

Sey per fies in the alloy is given by zv/1/2, so that for CuZn 

= 3/74/2=0-6752. The lattice constant a=2-945 A.U. and therefore 
ai 4159.10! dynes/cm?. 

No further progress can be made without a knowledge of the form of 
the function /(z). We shall assume a form of the type used by Leigh 
(1951) viz. f(z)=22—Az2’* in which A<1. This satisfies the conditions 
J (2) 2" as z+0 and f’(1)=0, and since f(1)=1—A we see that Ap? is the 
amount by which the energy at a plane of the BZ is less than the 
corresponding energy of a free electron. Using a=2-945 A.U. and 
expressing the energy in electron volts we find that the diminution of 
energy is equal to 8-67A) ev. 


With this form for f(z) the integrals w, and w, can be evaluated and we: 
find from (14) and (15) 


Je, F 2, 4(1—A) 
(1 K Ne — —1 ys : 3A 
tig) t+ gaqy) | 243A 
2(1—A) 
3242) 


pce DE nan 
r 11(1—A)(8-++-38A-++ 2722) 
22+ 3A)¥— SOF AYO SAEPA). 

~ alte (17) 


(16) 


(1/x)3(¢q1—C 19)" = — 82+ 
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Fig. 1 shows the two shear constants, for values of A from 0 to 0-2, calcu 
lated according to these expressions. When A=0-15, for empl the 
lowering of the energy of a state at the zone boundary would be 1-30 ev 
and the shear constants are then ¢,,"=1-8 and ec’? =1-2(10U dynes/cm?). 
It can be seen from the figure that, for values of \ which are at all plausible, 
the shear constant c’” will lie somewhere between 1 and 2 (101! dynes/em?). 
The Fermi energy appears therefore to provide the required stabilizing 
factor against the homogeneous shear (110) [110]. . 

The above calculation of the contribution from the Fermi energy to 
the shear constants depends on the assumption that the effect of each 
pair of zone boundary planes can be calculated separately. This 


Fig. | 


x10" dynes/cm? 


0 0:05 0:10 0-15 rl 0-20 


The Fermi energy contributions to the shear constant according to (16) and (17). 
Upper curve c,,", lower curve $(¢4;—¢12)". 


assumption should not lead to serious inaccuracy when the Fermi surface 
touches the BZ boundaries over only a small fraction of the whole surface. 
Clearly if the zone were almost completely filled the assumption would 
be quite invalid. P ° 
In the special case when A=0, so that the surface of the Fermi 
distribution is a sphere cut off by the planes of the BZ, it is possible to 
calculate the contribution to the shear constants without making the 
above assumption. An elementary calculation, which need not be 
reproduced, shows that in this case 
4? W—= R41 Sp,—tXpy—4h, . . . (18) 
2n2n= R? Sp,—t 2p — PR, . . . . . (19) 
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where R denotes the radius of the sphere which forms part of the surface 
of the electron distribution. The p, are given by equations (9) and (12) for 
the two different shears, from which the sums taken over the six pairs of 
planes forming the BZ can be calculated. In date way we obtain 


(Vje)og, (a — BEANE 8 ee 
and 
(1/x)4(6,,——-Cy5) == — gee 2 oe 
where € is given by the equation 
(10/3) 6Z24+v+2=0. . y . . = - (22) 


The superscript (0) is used to indicate the special case of free electron 
energies within the BZ. 

If we put = sec 6 then @ is the angle which the radius of the circle 
of occupied states on the (110) plane subtends at the centre of the zone. 
For f-brass, using the previously given value of v, we find from 
(22) €=1-0042 so that 0=5°15’. 

If we write £=1+6 equations (20) and (21) give, to second powers 
of 5 


(1/K)eq, (= 26 — 262, and (1/x)$(¢,;—C12) = 5 — 307, . . (23) 

and from (22) we find 
yon 2/34.25-2483,, 8 eee 
so that on substituting this into (16) and (17), taken when A=0, we find 
(1/4 = 25— 582, and (1/k)4(¢,,;—C,) =d —(9/2)87. . (25) 


It appears therefore that the approximate values of the shear constants, 
as given by (16) and (17), differ from the directly calculated values only 
by second powers of 6. The fractional error is the same in both cases viz. 
(3/2)5 which is 0-63%. Thus the approximation involved by 
adding the effects of the different pairs of planes acting separately is 
negligibly small for very small values of A, and it would appear unlikely 
that, over the whole range shown in fig. 1, there can be any serious error. 
Thus, if the energy discontinuities over the BZ boundaries in f-brass are, 
as is generally believed, about 2 or 3 electron volts, we are forced to the 
conclusion that the Fermi energy contributes at least 101 dynes/cm? to 
the shear constant c’ and thus stabilizes the structure with respect to 
homogeneous shears over the {110} planes. 


§3. THE CONSTITUTION OF THE SHEAR CONSTANTS AND THE STABILITY 
OF THE B-PHASE 


Besides the Fermi energy contribution there are two other important 
terms in the shear constants: an electrostatic term which contributes 
171-2 a-4Z¢? and 23-0 a~4Z,_? in units of 1011 dynes/em?2 to the two 
constants ¢,, and c’ respectively (a in A.U.), and the non-Coulomb 
interaction between the ionic shells. If a*w(r)/8 denotes the potential 
energy of a pair of ions at a distance r apart, arising from the non-Coulomb 
forces, the contributions which these forces make to ¢,, and c’ are 
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respectively (4/9)r2w”+-(8/9)rw’ and (4/3)rw’. It is known that w(r) varies 
very rapidly with r and it is usual to assume that the variation may be 
represented by exp (—r/p), where p is a length approximately equal 
to 0-3 A.U. Let rw’=X, then r2w’= —(r/p)X and the two contributions 
to cy, and c’ become (4/9)(2—r/p)X and (4/3)X respectively. 

To determine the electrostatic contribution we need to know the 
value of Z.~¢ which is very difficult to estimate accurately. It must, 
of course, lie between 1 and 2, and for the sake of obtaining approximate 
numerical results we shall simply assume Z,,2—=2-0. For the same reason 
we shall assume A=0-15 which would correspond to what is believed to 
be a likely energy gap over the BZ of about 3-0 ev, so that the Fermi 
energy contributions according to fig. 1 are c,,=1-8 and c’=1-2 
(1011 dynes/em?). ~ 

The actual values of the shear constants of B-brass have been measured 
by Lazarus (1949) who gives 


C4g—= 8-24, c’=0-97 (104 dynes/cm?). 
Adding the three contributions together and equating to the observed 
values we get the equations 
1-2+-0-61-+ (4/3) X =0-97, 
1-8-+4-4-55+ (4/9)(2—r/p)X=8-24. 

From which, since r=2-549 A.U., we find X=—0-63.10" ergs/em? and 
p=0-29 A.U. This value of X corresponds to a positive energy of 
1-33 keals per mole arising from the repulsive forces between the ions, 
which appears to be quite reasonable, as does also the value of p. The 
contribution to the bulk modulus from the non-Coulomb ionic interaction 
is readily seen to be —(4/9)(2+7/p)X which, according to the above 
values, amounts to 3-01.10! dynes/em?. The observed value is 
11-62.10% dynes/em? whilst the Fermi energy contribution, assuming 
free electrons, is 6-54.10'! dynes/em? (Jones 1949). Since a change of 
volume affects the energy of the lowest state of the conduction electrons, 
and therefore a term occurs in the bulk modulus which is absent in the 
case of pure shears, the agreement is as good as can be expected. 

According to this theory of the constitution of the elastic constants, 
therefore, the B-phase is rendered stable with respect to homogeneous 
shears across the (110) planes in the direction [110] by the effect of the 
zone boundaries on the Fermi energy. We are forced, therefore, to the 
conclusion that even the disordered B-phase would be mechanically 
stable at low temperatures with respect to any homogeneous shear. The 
reason that this phase does not occur at low temperatures is evidently 
because it becomes thermodynamically unstable relative to neighbouring 
phases, e.g. the ordered f’-phase or the ¢ and y types (Hume Rothery, 
Reynolds and Raynor 1940). 

If, as Zener (1947) suggests, the characteristic temperature of the 
B-phase is less than that of neighbouring phases the free energy of the 
B-phase decreases, as the temperature rises, relatively more rapidly than 
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that of the other phases. This would explain the occurrence of the 
B-phase only at high temperatures and the general V shape of the phase 
boundaries. 

The question now arises why it is so rarely, if ever, possible to quench 
the disordered B-phase and so obtain this structure in the metastable 
state. We suggest that the reason for this lies in the peculiar nature of 
the stabilizing influence of the Fermi energy. This operates only with 
respect to homogeneous shears over regions containing a good many 
atoms. It is, therefore, quite different from the stabilizing effect of 
forces which operate mainly between nearest neighbouring atoms. At 
any temperature, however low, local distortions should occur very easily 
in the B-phase and hence make the existence of any metastable state very 
- improbable. 
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X. Quantitative Measurement of Preferred Orientation in 
Rolled Uranium Bars 


By G. B. Harris* 
[Received September 26, 1951] 


SYNOPSIS 


The preferred orientation: of uranium bars hot-rolled and cold-rolled 
in the « range has been determined quantitatively by x-ray diffraction 
using a Geiger counter. It is shown that hot rolling with a slight reduction 
causes the 010 planes to become strongly and the 110 planes weakly 
oriented perpendicular to the rolling direction. Further hot rolling 
causes the 110 planes to become strongly and the 010 planes weakly 
oriented perpendicular to the rolling direction. Cold rolling with light 
reduction causes the 010 planes to become strongly and the 130 planes 
weakly oriented in the planes perpendicular to the rolling direction. 


$1. INTRODUCTION 


‘THovucH the tendency of the crystals to align themselves in particular 
preferred orientations in rolled or otherwise mechanically deformed metal 
bodies has long been recognized, and techniques for qualitatively 
observing preferred orientations by x-ray diffraction have been 
extensively developed, Barrett (1943), comparatively little theoretical 
or experimental work has been done on the COPE determination 
of preferred orientation by x-ray methods. 

The Geiger counter X-ray spectrometer, in which the intensity of 
x-rays diffracted from a specimen is measured directly with a Geiger 
counter, has obvious advantages for this work over the more usual 
photographic method of recording the diffracted x-rays, and has been . 
-successfully applied to it by Decker, Asp and Harker (1948), and by Field 
and Merchant (1949). In both these investigations the Geiger counter 
was set to pick up x-rays diffracted from one set of lattice planes, and 
the specimen, in the form of a flat plate, was rotated about two axes 
perpendicular to the x-ray beam, so that the variation in the Geiger 
counter response enabled a pole figure to be drawn, showing the relative 
amounts of crystalline material with these lattice planes oriented at 
different ‘latitudes and longitudes’ relative to the normal of the plate. 
‘The results had to be corrected for the different geometrical diffraction 
conditions at different settings of the plate ; Decker, Asp and Harker 
_ verified their calculated corrections by measurements on the silver 
grains in an exposed x-ray film, which were assumed to be randomly 


oriented. 
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With «-uranium, the comparatively low symmetry of the crystal 
structure (orthorhombic, a=2-852 4, b=5-865 a4, c=—4:9454; space 
group Cmem) enables an alternative method to be used, which yields a 
pole figure more directly related to the anisotropy of physical properties 
associated with preferred orientation. Because of the low symmetry, 
there are a large number of different sets of planes in the crystal lattice, 
which have slightly different spacings by which they can be distinguished. 
If the Geiger counter and specimen are simultaneously rotated in the 
spectrometer through the Bragg angle settings for these planes, it is 
possible to measure directly the relative amounts of crystalline material 
with each set of planes oriented parallel to the face of the specimen. 
Hence a pole figure can be drawn, showing the relative preponderance 
or scarcity of crystals with different crystallographic directions oriented 
along one axis in the specimen. This type of pole figure is convenient 
for relating the value of any directional physical property (e.g. electrical 
conductivity) measured in this direction to the values along the principal 
axes in a single crystal of the metal. 


Fig. 1 


TRANSVERSE TANGENTIAL RADIAL 
Types of specimen section examined. 


The method outlined above has been applied to the measurement of 
preferred orientation in cast, hot-rolled and cold-rolled uranium bars of 
diameters ranging from 0-7-1-5 in. (1-9-3-8 em). 


§2. PREPARATION OF SPECIMENS 


Specimens for the study of preferred orientation were cut from uranium 
bars with a carborundum slitting-wheel. The three types of specimen 
examined are illustrated in fig. 1. j 

A transverse section was obtained by cutting a bar along two parallel 
planes perpendicular to its axis and examining a plane face of the slice 
removed. 

A tangential section was obtained by cutting a transverse slice of a 
bar along a plane containing its axis and examining the face of this 
cut in the semi-cylindrical specimen obtained. The normal to the plane 
of the section at any point is perpendicular to both the axial direction 
and the radius passing through the point. 
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A radial section was obtained by cutting a tangential specimen along 
a plane parallel to the axis of the bar, but at some distance from it, 
and examining a region of the section half-way between the edges in 
which it meets the cylindrical surface of the bar. In this region the 
normal to the plane of section is along a radius of the bar. 

In practice the above geometrical conditions were fufilled with an 
accuracy of +3° in all the sections cut. 

All sections were hand-polished, then polished electrolytically and 
examined under the microscope between crossed polarizers. The detailed 
procedure was identical with that described by Chipman (1950). Grains 
of different orientation appeared with different brightness, because of the 
anisotropic reflectivity of uranium, so that the grain size of the specimens 
could be estimated, and the presence of preferred orientation inferred 
when the contrast between grains was weak. 


§3. PHoToGRAPHIC X-RAY EXAMINATION 


Before quantitative measurements were made on the Geiger counter 
spectrometer, all the specimens except those from bar 287/10 (see below) 
were examined by photographic x-ray methods. On a photograph it 
is possible to estimate, from the ‘ spottiness’ of the diffraction lines, 
whether enough crystals are covered by the x-ray beam to give a true 
picture of the statistical distribution; with the Geiger counter 
spectrometer anomalies caused by covering too few crystals are detected 
only after measurement of the integrated intensities and calculation of 
the pole figure. 

x-ray diffraction photographs of sections were taken in 11cm 
diameter cylindrical cameras with filtered cobalt Ka radiation (wavelength 
1-790 a). Specimens were placed with the plane of section on the axis 
of the camera at an angle of 20° or 38° to the incident beam. With an 
angle of 20° lines from the following lattice planes were focussed on the 
film: (020), (110), (021), (002), and (111). Lines focussed with an 
angle of 38° were (131), (040), (023), (200), (041), and (113). 

By comparing the relative intensities of the lines in sections of cast 
bars (assumed to be randomly oriented) and rolled bars, the preponderance 
or scarcity of crystals with the corresponding lattice planes parallel to 
the plane of section in the rolled bars could be estimated. 


$4. GEIGER COUNTER SPECTROMETER MEASUREMENTS 


In order to obtain a more detailed picture of their preferred orientation, 
the specimens were examined on a ‘Geiger counter spectrometer. In 
this instrument the x-rays diffracted from a specimen fall on a Geiger 
counter which, through suitable electronic circuits, supplies a voltage 
proportional to the x-ray intensity to an electronic potentiometric 
recorder. The specimen face being at an angle @ to the beam, diffracted 
rays are brought to a focus in front of the Geiger counter on a slit at an 


WZ, 
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angle 26 to the incident beam. Over a period of several hours the X-ray 
tube current slowly fluctuates, sometimes by as much as 10%. To 
compensate for these changes, which are ascribed to frequency changes of 
the A.C. supply, the measured area of each peak is divided by the tube 
current measured as it was recorded. Using filtered copper Ka radiation 
(wavelength 1-542 4) and values of @ up to 45°, lines of lattice spacing 
greater than 1-11 a could be observed. 

It is necessary to observe the orientations of a large number of crystals 
in order to be sure that any departure from randomness is statistically 
significant. The main criterion used in this work for deciding when 
enough crystals had been observed was the absence of spottiness in the 
Debye-Scherrer rings on a photograph. 

Smooth rings could not be obtained from stationary specimens of either 
cast or rolled uranium. Consequently specimen holders were made in 
which the specimen face was translated or rotated in its own plane, so 
that the x-ray beam scanned a fairly large region of the surface. In view 
of the well known fact that in rolled bars and slabs of many metals, 
deformation is greater in regions adjacent to the rolls than in the centre 
of the bar or slab, Orowan (1943), Kasz and Varley (1950), arrangements 
were made to scan narrow bands on the surface at constant distance 
from the axis of the bar. 


§5. INTERPRETATION OF GEIGER COUNTER SPECTROMETER RESULTS 


Preferred orientation in a polycrystalline specimen, being a non-random 
distribution of orientations of the crystal axes relative to directions in 
the specimen, can have a definite meaning only if the part of the specimen 
examined contains a statistically significant volume of crystals in every 
orientation. Furthermore, any systematic variation of crystal orientation 
between different regions of the specimen can be expressed only by 
determining separately the preferred orientation in a number of regions, 
each containing enough crystals to satisfy this condition. 

These considerations inevitably limit the accuracy in practice of the 
quantitative treatment described below, which appear to be implicit 
in most published work on preferred orientation. The accuracy of the 
method increases with the ratio of specimen size to grain size, other things 
being equal. 

In the Geiger counter spectrometer measurements described above, 

_ X-ray reflections contributing to a particular line were received from all 
crystals with the corresponding lattice planes parallel to the specimen 
face, i.e. from all crystals in which the a, 6 and ¢ crystallographic axes 
are at three particular angles, «, 8 and y respectively to the normal to 
the specimen face. By measuring the intensities of a large number of 
lines, with different values of «, 8 and y, it is possible to build up a picture 
of the relationship between these angles and the volume of crystals with 
the corresponding orientations, which can be consistently represented by 
a standard stereographic projection of the crystal axes and related 
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directions in a crystal of «-uranium, with contours superposed to represent 
the preponderance or scarcity of orientations in which the various 
crystallographic directions are perpendicular to the plane of section. 
Such a diagram will be called an inverse pole figure (e.g. figs. 2-4). If 
obtained from a transverse section, so that it shows which crystallographic 
directions are oriented along the rolling direction, it will be called a 
W inverse pole figure (German ‘ Walzrichtung’). If obtained from a 
tangential section, so that it shows which crystallographic directions are 
oriented along a tangential direction, it will be called the Q inverse pole 
figure (German ‘ Querrichtung’). These names are consistent with the 
nomenclature used in the literature to describe orientations in rolled 
sheet, Barrett (1943), if it is assumed that the compression direction is 
radial in bar-rolling. 

The fraction of the total volume of material observed, which consists 
of crystals whose orientations are within a small solid angle dQ about a 
point («, 6, y) on an inverse pole figure, can be expressed as 


pla, B, y) AQ) 427, 


where p(«, 8, y) is a quantity which, in general, depends only on «, 6 and y. 
Since every crystal must have some orientation, it follows that 


Z| | ol 8) ao=1, 


the integral being taken over the whole inverse pole figure. For randomly 
oriented material p is independent of «, 8 and y, so that, in this particular 


case — A 
Pp 
Lf = jhe 
eae ; 


p=l. 


Values of p greater than 1 at any point show that the corresponding 
crystallographic direction is preferentially oriented perpendicular to the 
plane of section; values less than 1 show that such orientations are 
‘ preferentially avoided ’. 

If p¥(«, 8, y) is the p value of a point on the W inverse pole figure and 
pa, B, y) is the p value of a point on the Q inverse pole figure, and if 
all directions in the bar perpendicular to the rolling direction are 
equivalent in respect of orientations, it can be shown that 


hence 


1 ies , , / 
D%(a, B, y= xe | BM, Bi 7’) Ad 


where («’, B’, y’) lies on a great circle perpendicular to («, B, y) and d is 
its angular distance from some arbitrary zero on this great circle. A 
perfect ‘fibre’ texture in a rolled bar can thus be recognized by 
examining two sections. 
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6. CALCULATION OF p VALUES FROM INTENSITY MEASUREMENTS 


Tt can be shown, Field and Merchant (1949), that the intensity of the 
x-ray line diffracted from the surface of a thick specimen, when a beam 
is incident on it at the Bragg angle 6 for the particular line examined, is 


given by 
Lj=AD/2p 


where A is the intensity of the incident beam, 
D is the fraction of the incident intensity coherently diffracted 
by the specimen per unit path-length, 
uw is the x-ray absorption coefficient of the specimen. 
I, is thus independent of the Bragg angle 6, other things being equal. 
The quantity D can be expressed as the product of three factors, viz. :— 
2 
n(a, B, y) ase 


the fraction of the total volume of material oriented within an angle 60 
of the appropriate orientation (86 being the effective ‘ line breadth ’) ; 

F?(hkl), a factor depending on the crystallographic indices of the line, 
having the same value for a particular line in all specimens of the same 
material (the ‘ F? value’ used in crystal structure analysis by X-ray 
diffraction) ; 

B, a factor depending on the imperfections in the structure of the 
crystals, constant for all lines from the same specimen. 

Thus 


‘AAJ us 9 
ORY a, B, y) F(hkl), 


for a given set of reflecting planes of index (hkl). B and « are constants 
and experiment showed that for a given specimen 60 is also constant. The 
quantity J(hkl) measured by the Geiger counter spectrometer is equal 
to KI,/A where K is an instrumental constant and hence 


I(hkl)=C . p(«, B, y)F2(hkl) 


where C is a constant which includes B, x, 66 and K and may vary from 
specimen to specimen. 
For a randomly oriented specimen 


I' (hkl) =C' F2(hkl). 


Hence the ratio of the line intensities from preferentially and randomly 
oriented specimens is given by 

(hkl) C 

T (hkl) See p(x, B, y). 


To eliminate the factor C/O’, it can be assumed that, as the mean value 
of p(x, B, y) for all orientations is 1, the ratio of the sums of the intensities 
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of a large number of reflections for preferentially and randomly oriented 
specimens is equal to C/C’. Hence 


3 (hkl) ZI'(hkt) 
pla, B; y)= DT (hkl) . T(t) * 


Values of I’(hkl)/XT'(hkl) were obtained by measuring the intensities of 
20 lines from transverse, tangential and radial sections of 1 in. cast bar 
GH/9. The value for each line from each section, the mean of the three 
values, and the corresponding value deduced from Jacob and Warren’s 
calculated line intensities are given in table 1. There are variations up 
to a factor of 2 between intensities of the same line from different sections, 
but, as they do not fit into consistent pole figures, these anomalies can 
be attributed to insufficient scanning. The mean values, which were 
used in subsequent calculation of p values, do not differ greatly from 
Jacob and Warren’s calculated intensities ; the differences can probably 
be reconciled by small adjustments in the assumed crystal structure of 
#-uranium. p values calculated from the experimental figures may be 
expected to be correct to +20% 


RESULTS 


The preferred orientations observed are summarized in table 2. 

Bar 24692. This bar, 14 in. (3-8 em) in diameter, was hot-rolled 
in the «-range from a cast billet. x-ray photographs of a transverse 
section, taken with scanning annuli of mean radii 0-15, 0-3 and 0-45 in. 
(0-4, 0:75 and 1:15 cm), showed enhancement of the (020) line and 
reduction of the (002) line, both effects increasing with the radius. These 
changes indicate preferred orientation of the (010) planes perpendicular 
to the rolling direction ; the degree of preferred orientation increases with 
distance from the axis of the bar. _ 

Geiger counter intensity measurements were made on transverse and 
tangential sections at a mean distance of 0-4 in. from the axis. The 
results are summarized in fig. 2, the W and Q pole figures. The 
(010) planes are strongly oriented, and the (110) planes very weakly 
oriented, perpendicular to the rolling direction; the (113) planes are 
oriented radially in the region examined. However, since the (010) and 
(113) planes are not quite perpendicular, the W and Q pole figures are not 
quite consistent with each other, though each is consistent in itself. 
The anomaly can be explained if the preferred orientation is different 
along different tangential directions, but it can only be inferred from the 
results that there are strongly preferred orientation directions 
perpendicular to the axis of the bar, i.e. that the texture is far from a 
perfect ‘ fibre ’ texture. 

Bar 067921. This bar, lin. (2:5 cm) in diameter, was hot-rolled in 
the a-range with a heavier reduction than bar 24692. An X-ray 
photograph taken with a scanning circle of 0-25 in. radius shows great 
enhancement of the (110) line and slight enhancement of the (020) line, 


f 
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Table 1. Values of [(hkl)/T(hkl) for Sections of Cast Bar GE/9 


I(hkl)/XT (hkl) 
hkl Jacob and 
Transverse | Tangential Radial Meat Warren’s 
Section Section Section calculated 
value 

114 0-0203 0-0360 0-0415 0-0326 0:0183 : 

133 0-0271 0-0419 0-04.26 0-0372 ~ 0-0290 

221 

on 0-0758 0-1439 0-1008 0-1068 0-0761 

042 0-0139 0-0157 0-0119 0-0138 0-0171. 

132 0-0109 0-0083 0-0077 0-0090 0-0077 

113 0:0227 0-0225 0:0285 0-0246 0-0239 

041 0-0105 0-0134 0-0245 0-0161 0-0068 

200 | 0-0223 0-0177 0-0101 © 0-0167 0-0162 

023 0-0408 0-0264 0-0217 0-0296 0-0325 

040 0-0114 0-0141 0-0239 0-01 65 0-0136 

131 0-0791 0-0616 0-1236 0-0881 | 0-0735 

130 — 0-0061 0-0063 0-0089 0-0071 0-0077 

112 0-1045 0-1056 0-1364 0-1155 0-0853 

022 0-0061 0-0124 0-0090 0-0092 0-0051- 

lll 0-0997 0-1345 0-0891 0-1078 0-1247 

002 0-1303 0-0649 0-0884 0-0945 0-1041 

021 0-1680 0-1144 | 0-1052 0-1292 0-2140 

110 0-1407 0-1534 0-1128 0-1356 0-1485 

020 0-0098 0-0069 0-0129 0-0099 0-011 1 | 
ee ee eee eee 

024 


29 was not detected in either section 
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Table 2 


Bar No. Treatment 


Preferred Orientation 
Planes perpendicular to 
rolling direction 


24692 | Hot-rolled with light 
reduction 


(010) strong 
(110) weak 


067921 Hot-rolled with heavier 
reduction 


(110) strong 
(010) weak 


-287/10 Cold-rolled with light 
reduction 


O10 ben 


(O10) strong 
(130) weak 


12t 


Pole figures showing distribution of crystallographic plane normals in bar 24692 
after light reduction by hot-rolling, (7) along the rolling direction, (5) along 


a tangential direction. 


¢ 
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with weakening of the (021) line and disappearance of the (002) and 
(112) lines. These changes indicate a strong orientation of the (110) planes, 
and weak orientation of the (010) planes, perpendicular to the rolling 
direction. 

Geiger counter intensity measurements were made on transverse and 
tangential sections at mean distances of 0-25 in. and 0-2 in. respectively 
from the axis of the bar. The results are summarized in fig. 3. The 
W pole figure confirms the photographic results, and the Q pole figure 
shows the (023) planes radially oriented. Since the (023) and (110) planes 
are not quite perpendicular, it can only be stated that there are preferred 
orientation directions perpendicular to the axis of the bar. 


Fig. 3 
1:06 - 


Pole figures showing distribution of crystallographic plane normals in bar 067921 


after heavy reduction by hot-rolling, (a) along the rolling direction, (b) along 
a tangential direction. 


Bar 287/10. This bar was cold-rolled with a light reduction. 
Micrographic examination showed a large grain size without perceptible 
elongation of the grains along the rolling direction. Line intensities were 
measured on the Geiger counter spectrometer with transverse and 
tangential sections, at a mean distance of 0-1 in. from the axis. The 
results are summarized in fig. 4. The W pole figure shows that the 
(010) and (130) planes are oriented perpendicular to the rolling direction 
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and the Q pole: figure shows that the (100) 
perpendicular to ( 
examined, 


planes, and all planes 
130) are preferentially oriented radially in the region 


WY 


1-16 


0-42 


MY 
\'20-3 0:38 0-15 


“a, 
Le 5 
oO 


are 
S 


» Pole figures showing distribution of crystallographic plane normals in bar 287/10 
after light reduction by cold-rolling, (a) along the rolling direction, (6) along 
a tangential direction. 
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XI. CORRESPONDENCE 


Metallic Conduction—The Internal Size Effect. An Addendum 


By D. K. C. MacDonaLp* 
Clarendon Laboratory, Oxford+ 


[Received July 26, 1951] : 


In a recent paper under this title (MacDonald 1951), the writer has 
pointed out that internal ‘ barriers’ or scattering boundaries in a metal 
may, in particular, result in anomalous behaviour under an external 
magnetic field. Thus, specifically, an initial decrease of resistance can 
take place when a magnetic field is applied ; the Hall field characteristics. 
may also deviate significantly from those appropriate to an ideal, 
homogeneous, bulk metal. We wish here to make one or two observations. 
in the latter connection. 

First, let us mention a recent note on ‘ The Hall Potential across an 
Inhomogeneous Conductor’ by J. Volger (1950). Volger is not, concerned, 
as we were, with the situation arising in a region of metal when the 
electron mean free path becomes comparable with the size of the region ; 
in general, however, whenever the length of an element is comparable 
with the breadth (across which the Hall field is generated) deviations from 
the simple theory may arise and Volger also suggests that this may be of 
significance when interpreting experiments on inhomogeneous conductors 
—particularly granular semi-conductors. 

Secondly, we wish to amend somewhat our qualitative discussion of * 
the generation of the Hall field in a ‘ striated’ conductor, where the 
striation-separation is comparable with the electron free path. An 
internal boundary or ‘ striation® might be regarded as having various 
limits of behaviour, as follows. Firstly, if all electrons penetrated the 
barrier without significant scattering, elastic or otherwise, then there will 
of course be no effect either on the conductivity or on the Hall field. If, 
however, electrons were reflected elastically then the conductivity in the 
absence of a magnetic field would be unaffected but the Hall field arising 
when a magnetic field is applied would take the form sketched in our 
fig. 1(b). (loc. cit.)}. Diffuse, inelastic, reflection will again cause the Hall 
field to take the same qualitative form, but now also the resistivity (in the 
absence of a magnetic field) will be increased over that of the homogeneous 


* Communicated by the Author. 

+ Now at the Physics Division, National Research Council, Ottawa, Canada. — 

{ We may remark that, in recent personal discussion, Professor C. J. Gorter 
has questioned our assumption (MacDonald and Sarginson 1950) that the Hall 
field must fall to zero as the external boundary of a metal is approached, 


suggesting the possibility of additional ‘locally bound’ electrons forming a 
true surface charge-layer, 
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metal; application of a magnetic field can then in general produce a 
decrease. Finally, diffuse scattering will again increase the resistivity 
over the homogeneous bulk metal but the Hall field will no longer fall to 
zero at the internal boundaries, but may be expected to follow qualita- 
tively the variation of our fig. 1 (c) ; again we may expect that a magnetic 
field can produce some initial decrease of resistance although this will be 
less pronounced than in the corresponding case of diffuse reflection. 

In general, experimental studies of the Hall effect in metals. at low 
temperatures, especially those exhibiting the anomalous minimum of 
electrical resistance, should be of considerable interest. 
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The Moving Griffiths Crack 


By C. J. TRANTER 
Military College of Science, Shrivenham+ 


[Received November 1, 1951] 


In her recent paper ‘ The Moving Griffiths Crack ’, Dr. E. H. Yoffe (1951) 
reduces the problem of the calculation of the stress field about a straight 
crack moving through an elastic medium to that of finding a function A(s) 
which satisfies the dual integral equations © 


{ sA(s) cos xs ds=constant (—a<u<a), 
0 


ne f 
A(s) cos as ds=0 (|e | =a). | 
J0 

These equations are included in the more general pair considered by 
L. V. King (1935-36) and I have recently given solutions (1950, 1951) 
to still more general pairs. It is the purpose of this note to point out 
that the solution of the above dual equations can be obtained immediately 
and more simply from King’s result. 
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High-energy Components of the B-particle Spectrum of MsTh, 


By C. G. Camppety, W. J. HenpERson* and J. KYLEs 
Department of Natural Philosophy, University of Edinburgh 


[Received October 17, 1951] 


A new investigation of the B-particle spectrum of MsTh, has been made 
using one half of the double semi-circular: focusing spectrometer of ; 
Feather, Kyles and Pringle (1948). 

The method of source preparation followed closely that outlined by 
MacLane and Peterson (M.D.D.C. 1742) and Peterson (M.D.D.C. 1709). 
In order to obtain activities of the order of 1 mc the cerium carrier was. 
not removed from the MsTh,, since this process proved to be one of low 
yield, but the amount of carrier added was kept down to about 0-015 mg. 
The nitrates of cerium and MsTh, were evaporated on 0-2 mg/cm? 
aluminium foil of area 1 cm by 0-2 cm. These sources were used in the 
spectrometer at an emission angle of 45° so that the effective breadth 
was reduced and sharpness of focusing was improved. 

Scintillation counters were employed as detectors (E.M.I. Type 5045 
photomultipliers with anthracene as the phosphor). Tested with 
monoenergetic B-particles these counters gave a plateau in the counting 
rate/multiplier voltage curve, starting at a voltage which increased as the 
energy of the incident B-particles decreased. 

The magnetic field was determined from a Fermi plot of the 
ThC spectrum obtained with a source of: thorium active deposit, the 
end-point for ThC being taken as 2-256 Mev. That the MsTh, sources 
were free from the decay products of thoron, was shown by the fact 
that there was no evidence of the F-line of the ThB spectrum. Moreover, 
for B-particles of 1-80 Mev energy the half-life o the sources used agreed 
with the accepted value for MsTh, to within 1% 

For the investigation of the high- energy romion of the MsTh, apecuatt 
the multiplier role was set so that unnecessary sensitivity to low 
energy (scattered) particles was avoided. This adjustment at the same 
time reduced the effective y-ray background registered by the counter. 
Only as low energy regions were explored was the multiplier voltage 
appropriately increased. 

A Fermi plot of the spectrum was constructed and was found to be 
non-linear except near the end-point. Extrapolation of the final linear 
portion gave a maximum energy of 2:16-0:03 mev. This value is 
considerably greater than the values of Lecoin (1938) and Lee and Libby 
(1939), but is in reasonably good agreement with the absorption 
measurements of Feather (1930, 1938). An attempt was then made to 


deduce the end-points of possible partial spectra by resolving the Fermi 
—ee—__ 
* On leave from the National Research Council of Canada. 
+ Communicated by Professor N. Feather. 
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plot into straight lines, and this showed evidence of partial spectra with 
end-points at 1-6-+0-1 Mev and 1-85-L0-1 Mev. At lower energies the shape 
of the continuous spectrum is masked to a large extent by conversion 
lines, but there seems to be evidence for the existence of at least one 
other partial spectrum with an end-point at about 1 Mev energy. 
Coincidence experiments are now being undertaken in an attempt to 
establish these values more exactly. 
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A Note on the Statistics of the Weiss Field 


By G. M. Brut 
The Clarendon Laboratory, Oxford* 


[Received October 17, 1951] 


CoNSIDER an assembly of NV magnetic dipoles of magnitude p, with two 
orientations relative to an applied field, H. If N, are oriented along the 
field and N, against it, the relative magnetization is given by 
m=(N,—N,)/N. Suppose that the V dipoles occupy all sites of a lattice 
of co-ordination number q, that each parallel pair of nearest neighbours 
has an interaction energy, —w/2, and each parallel pair, w/2, where w>0. 

Let N,, be the number of anti-parallel nearest neighbour pairs. Then, 
leaving out certain constant terms, the nearest neighbour interaction 
energy can be written as V,,.w (Rushbrooke 1949). 

In a random arrangement N,.—qN,N./N=qN(1—m?)/4 and in a 
zero-th order approximation this is taken as the actual value of Nj. The 
nearest neighbours thus contribute —qNwm?/4 to the energy (omitting 
the constant gwN/4). A similar term, with different values of q and w, is 
given by each type of neighbour pair. (It should be, noted that, in the 
limiting case where a system’s interaction with any other system does 
not depend on the distance, the zero-th order method becomes 
accurate.) The interaction energy for the assembly is thus —NkAm?/2 
where k is Boltzmann’s constant and a constant depending on the lattice 
structure and exchange forces. 

Now it is well-known that, using classical (localized) statistics, m now 
has the same value as in an effective field H1, where H'=H-+-kAm/p. 
However, correct values of m and #,, cannot simultaneously be derived 
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from a magnetic partition function which depends on H? as if the latter 
was an external field in the theory of paramagnetism. This can be seen 
from writing 
—pH} 
(BF a= {exp er prt EXP a a} 


= 2 ak exp {( N,—N,) ey 


yi+m.-nW,! Ne! 

All values of NV, and N, satisfying the relation N=N,+ NV, are summed 
over, but all except one pair of such values are incompatible with the 
given value of m or (V,—N,)/N. 

There is the same inconsistency when the Weiss Field is used with 
Fermi-Dirac Statistics (Stoner 1938, 1939). However, a_ partition 
function from which correct expressions for the energy and magnetization 
are derivable can be constructed by summing, for each value of m, over 
only those sets of occupation numbers which are consistent with that 
value of m. Thus, in the case H=0, put 


pH+ a 


—AmEN f Jide 
12. Ie. =D. [ex - pe ex sits | 
iy 2T 2d,=3(1+m)N P ET: 
—2A,€ 
x px exp es i): We to te 2 
eta Pr } (2) 


The d, are the occupation energies of the microscopic energy states 
denoted by ‘s’ and can take only the values 0 or 1. The method of 
steepest descents is used to evaluate the sums inside the square brackets 
which are taken over the states with positive and negative spin 
component, respectively. It is then found that 


—dDd.< 00 e 
] ty ps Aisin me Ce : wipe - 
‘ ae Pane } i pt [exp 7 er} 2 | de—Nny 
(3) 


where n(e)de is the number of levels in the energy interval de and N, Is 
given by 


Jn [ex {er =n} +a] de=N,=}N(1+m). . . . (4) 
Similarly 
Li di€, ise € 
sa cHiemw Perf =], om [exfas— gp} +1] ean, 
% ; agers eats 
where, i (| exp gp —nop +1 | de=N,=}N(1—m). . . . (6) 


Eqns. (3) and (5) are now substituted into (2). Then, using 
(4), (6), (P.F.),, is differentiated logarithmically to obtain its largest 
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term, which by the usual property of thermodynamic functions yields 
its effective value. From the differentiation we have 


Grea AMD Oe Wee ee el Sr eer) 

Defining a new parameter 7 by 7= a(n +72), (4) and (6) become 
lr n(e) [xpd gp —n— amt} 1] de=4N(1+m) . . . (8) 
i; nle)| exp fp <b Am|P) 1] ae PIN any eo) 


and the free energy is given by 


F 
op Nat $NA? |/T— is n(e) In | exp {nbd |? — ert - 1] ile 


val As (c)In == €.. ; 
I, n(e) | exp {n— Am|T— in} +1 |e (10) 


It can be seen that a partition function constructed so as to sum over 
only consistent values of the sets of occupation numbers gives correct 
expressions for the magnetization and free energy by differentiation, 
without any necessity for the addition of extra terms. 
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The Angular Distribution of Protons in the Photo-Disintegration 
of the Deuteron 


By K. PHItures 
High Voltage Research Laboratory, Metropolitan-Vickers Electrical Company* 


[Received November 1, 1951] 


SUMMARY 


Using a thin heavy wax target the angular distribution of photo-protons 
has: been observed in photographic plates, for the energy range 6 to 


15 Mev. 


In the past few years many publionions have appeared on the experimental 
investigation of the angular distribution of photo-protons in the photo- 
disintegration of the deuteron. Most of the investigations have been 


* Communicated by the Author. 
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earried out using monochromatic y-ray sources of 6-1 and 7 Mev, e.g. 
Goldhaber (1948), Phillips, Lawson and Kruger (1950), Hough (1950), 
Gibson, Grotdal, Orlin and Trumpy (1951) and at 17-6 Mev by Hough, 
Waffler and Younis (1949). Another recent paper by Fuller (1950) 
describes an experiment using a 20 Mev betatron as x-ray source. These 
experiments especially at the higher energies are likely to yield important 
evidence concerning the nature of nuclear forces. 

This short report describes some preliminary results on a similar 
experiment using a 20 Mev betatron for the x-ray energy range of 5-15 Mev. 
A collimated beam of x-rays was incident on a 200-micron thick Ilford C.2 
nuclear emulsion over which was placed a 10-micron thick deuterium wax 
target. The plates were inclined at 15 degrees to the horizontal. The 
exposure to the betatron Jasted only for a few minutes so as to prevent 
undue x-ray blackening. The plates were processed in Azol developer 
and 1% potassium bromide in the usual manner. 

With the aid of a Vickers projection microscope several thousand tracks 
were examined and only protons with angles of dip in the emulsion within 
the interval 5-45 degrees were used in the distributions. The upper limit 
was adopted to reduce the error due to inaccuracies in range measurements 
for tracks with large angles of dip. The lower limit was introduced in 
order to eliminate large errors due to the passage of the protons through 
the wax target. A study of the neutron flux in the betatron enclosure 
was carried out. by examining emulsions exposed with ordinary paraffin 
wax targets. The effect of these neutrons was found to be almost negligible 
except at the very low energies, i.e. less than 5 Mev. From the geometry 
of the experiment it can be shown that the angle of the emitted proton 
with respect to the direction of the x-ray beam is given by the expression 


cos 0=sin 8 sin 15°-+-cos 15° cos B cos 4, 


where f is the angle of dip in the emulsion, corrected for shrinkage 
and 6 is the angle of ejection of the photo-proton measured in the plane 
of the plate. 

The protons were divided into four convenient energy groups 5-8, 
8-10, 10-13, 13-15 Mev, and the histograms of the number of tracks 
emitted for unit solid angle and corrected for loss due to the angle of 
dip criteria, are shown in figs. 1, 2, 3 and 4. The smooth curves are the 
sin’ 6(1+ By, cos @) distributions, predicted by the theory of Marshall 
and Guth (1950). 

Assuming a distribution of the form 

N(9)=Ayzaytsin? 6(1+ By,x cos 8) 
the values of the constant A,,;, and B;,, have been determined by 
statistical analysis of the data and are given in the table for the different 
x-ray energies. The values of Aj, indicate an increase of the isotropic 
term with x-ray energy as is characteristic of the parameter due to 
non-central interaction of the nucleons ; whilst the values of By. , agree 
with the theory of Marshall and Guth within the experimental error, 
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Fig. 1 
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The histograms show the angular distribution of the photo-protons between 
0° and 160°, with respect to the direction of the x-ray beam, for four 
different x-ray energies. The dotted curves indicate the expected 

distribution predicted by the theory of Marshall and Guth. 
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Energy 5-8 8-10 10-13 , 18-15 
hy Mev 


Arie 0-063-+0-015 | 0-083-+40-013 | 0-119+0-010 | 0-135+40-017 
Buss 0-22 +0-2 0-24 40-16 | 0-44 40-23 | 0-49 10-25 


The limits in the above table have been taken as plus and minus twice 
the standard error. 


Fig. 3 
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The histograms show the angular distribution of the photo-protons between 
0° and 160°, with respect to the direction of the x-ray beam, for four 
different X-ray energies. The dotted curves indicate the "expected 
distribution predicted by the theory of Marshall and Guth. 
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A Theorem involving Binding Energies in Relation to B-disintegration 


By. N. Fratusr, F.R.S. 
The University, Edinburgh* 


[Received November 12, 1951] 


From time to time many expressions, equivalent so far as they are 
essentially correct, have been put forward to relate the energy released 
in B-disintegration to the masses of the nuclei involved, the total energy 
of binding of the neutrons and protons in such nuclei, or the masses 
of parent and daughter atoms in the neutral state. Sometimes the 
difference between the disintegration energies of successive f-active 
bodies (as of the successive members of a ‘ chain’ of fission products), 
or the difference between the disintegration energies of succeeding 
members of a series of bodies of constant isotopic number, is of interest. 
In such cases the ‘ theoretical ’ expression for the energy difference may 
be neater than the corresponding expression for the energy of a single 


disintegration. In this note a form of expression is deduced for the 
[A van 
difference between 7’; ( Zs and 7’; gy 


B-disintegration of two species ( a and ( Ak which exhibits this 


, the total energies of negative 


neatness in a marked degree. So far as the writer is aware this expression 
has not previously appeared in the literature. 
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A A 
If the exact mass of the nucleus ( a) is represented by NV e and 


if m is the mass of the electron, we have 


re()-[G)-*eh) Jee?) 0 


Z representing the energy derived from the adiabatic 


Z 
the term e 
reorganization of the extranuclear electrons in the process of disintegra- 


A 
tion. Ifthe energy of binding of the last-added proton in the nucleus ( =) 


Z and that of the last-added neutron in the same 


ae a little: consideration will show that (1) becomes 


A 
is written as B, 


nucleus as b,, 


A A A 1 Z+1 
re(2)-m(d) M(B) (CE). 8 


Again, consideration of the process of formation of the nucleus 
A+1)\~ A—l | 

ae) from the nucleus ( Z *) , first by the addition of a neutron 
followed by a proton, then by the addition of a proton followed by a 
neutron, leads to the result 


A Awl A ALIN Se a 
2.) +2672, (Ga eee 


From (2) and (3) we obtain 


eee A+l1 A+l1 eat ve a8 | 
#5 (2) =o(ai) Pear) (CZ) 


and from (4), neglecting the very small difference « ic, ‘) ae ts a _) ; 


Z Z' 
a YAN 3, (Al A+] Atl 
5 (2) te (2) 8 (a1) 8 (zar1) Bo (Ze) 
Ate 
+B, (ean (5) 


Eqn. (5) is the expression referred to above. For positron disintegration 
(or capture transformation) the corresponding expression will be found 
to be 


OY ANe aE A+1 A+ A'4] 

ti (z) ti (a) =") 87") Bao) 
A’+1 
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When (5) and (6) refer to active species ah ( 


) having the same 


isotopic number D, it will be observed that in one case (5) the binding 


energies are those of neutrons and protons in nuclei both of isotopic 
number D—1, and in the other case (6) the binding energies are those of | 
neutrons and protons in two nuclei each of isotopic number D+1. It is 
this simplicity which makes the expressions particularly useful in 
discussing the variation of disintegration energy with mass or charge 
number for a series of bodies having the same value of D(=A—2Z). 
Recently the writer has given a survey of the total energies of 
B-disintegration of species of isotopic number ranging from 40. to 53 
(Feather 1951). For each value of A—2Z a curve is given showing the 
variation of 7; with neutron number A—Z. Certain of these curves 


_show discontinuities in the region A—Z=126, but the curve for 


A—2Z=45 does not. In applying (5) to this smooth curve we are 
concerned with the differences between the binding energies of the 
last-added protons and the last-added neutrons in a series of nuclei for 
which A—2Z—44. This isotopic number has unique’ significance in | 
the range 39 <A —2Z<52 as the difference between the ‘ magic ’ numbers 
126 and 82. The terms entering into.-(5) in this connection are the 
differences between the binding energies of ‘ equivalent’ protons and 
neutrons, in respect of shell closure at Z=82, A—Z—=126, respectively. 
The fact that the curve for 7; against A—Z for A—2Z=45 is perfectly 
smooth then implies that the contribution of specifically nuclear forces 
to the variation of proton binding energy with proton number around 
Z=82 is paralleled precisely by the contribution of specifically nuclear 
forces to the variation of neutron binding energy with neutron number 
around A—Z=126—at least within the series of nuclei for which 
A—2Z=44. This isa significant result to deduce from what is, in a sense, 
purely qualitative information. 

It is hoped to give a more complete discussion of the experimental 
material from this point of view in a later publication. 
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XII. Notices of New Books and Periodicals received 


Ferromagnetism. By R. M. Bozorru. * [Pp. 968.] (London: Macmillan & Co. 
Ltd.) Price £6 10s. (New York: D. Van Nostrand Company Inc.) Price 
$17.50. 


Tuts is an excellent reference book on ferromagnetism. It should be available 
to everyone working in a field related to ferromagnetism. The book is divided 
into four parts: a short introduction, one dealing with the properties of 
magnetic materials, another on magnetic phenomena and theories and a short 
one on the technique of measurements. This makes it of great use to engineers 
and physicists alike. The enormous amount of information is presented in a 
very pleasant way. The number of graphs is about equal to the number of 
pages and there are many tables and an up-to-date list of references covering 
70 pages. 

It is obvious that this work is not meant to be a textbook on ferromagnetism. 
Fundamental theories are only briefly reviewed but great attention is paid to 
the description of basic physical phenomena and to the theory necessary to 
understand the physical concepts. The results of recent ‘research are amply 
used in the presentation of the general picture of ferromagnetism. The author 
has wisely made no attempt to include theories of a more speculative nature 
‘which could be advanced for the explanation of many experimental 
observations. DAP. 


The Interpretation of X-Ray Diffraction Photographs. By N. F. M. Henry, 
H. Lipson and W.' A. Wooster. [Pp. 258.] (London: Macmillan.) 
Price 42s. 


Few books have hitherto been written about x-ray diffraction techniques which 
deal adequately with the indexing and interpretation of diffraction photographs, 
so that this new, authoritative and detailed practical textbook will be widely 
appreciated. In writing it, the authors have drawn jointly upon considerable 
experience gained not only in the practice but also in the teaching of their 
subject, and they have exercised sound judgment in the choice of subject matter. 
This is presented largely in a geometrical way but an elegant treatment is also 
given of the corresponding analytical approach, the style throughout being lucid 
though concise. 

The first two chapters deal in an elementary way with crystal lattices, their 
symmetry and representations, and with the fundamental physical properties of 
x-rays. In the third chapter, this leads naturally to the theory of x-ray diffrac- 
tion, into which the concept of the reciprocal lattice is introduced at an early 
stage. The methods of indexing and interpreting single crystal photographs are 
then discussed, and are illustrated by a number of examples in which the more 
complex problems are approached in a systematic and simple way. Other 
chapters are devoted to the powder method and the accurate determination of 
cell dimensions, the study of orientation in single crystals and twins, preferred 
orientations in polycrystalline aggregates, the calculation and measurement of 
the intensities of x-ray reflections and the measurement of grain size. The 
reproduction of photographs and charts to scale, together with the inclusion of 
worked examples and tables of the physical constants and mathematical 
functions most frequently required by the experimenter contribute greatly to the 
practical value of the book. It can confidently be recommended to all who are in 
any way concerned in the practice of x-ray diffraction. HeDsK 
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Figs. 5-10 
Showing sequence of growth about a single screw dislocation group. 
: Photographs at 10-minute intervals. x 250 
Figs. 11-13 
Growth from two screw dislocation groups of the same hand. 
Photographs at 10-minute intervals. x 250 
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Figs. 14-18 
Growth from two co-operating screw dislocation groups of opposite 
hand. Photographs at 15-minute intervals. x 300 

Figs. 19-23 


Growth from two co-operating screw dislocation groups of opposite 
hand. Photographs at 60-minute intervals, (Same crystal five 
days later.) x 300 
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Fig. 27 
Illustrating cross-lacing phenomenon and domination of a pair of 
screw dislocation groups of opposite hand. 


(Reflection phase contrast illumination.) (Step-height 294 4+ 1-5 A.) 
x 650 


DESCRIPTIONS TO PLATE III 


Fig. 24 


Showing step-lines on the surface of a crystal containing many 
screw dislocation groups. x 250 


Fig. 25 
Crystal plate with spiral step-line on one surface and radial step-line 
on the other, both centred on same screw dislocation group. 500 
Fig. 26 
Crystal plate with spiral step-line on both surfaces, 250 


